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Part  I of  the  report  deals  with  stick-slip  of  selected  rocks  at  different 
t-  r.peratures  and  under  normal  stresses  up  to  20  kg/cm^.  Stick-slip  dominated 
the  sliding  behavior  of  roughly  ground  surfaces  of  limestone,  basalt,  and 
quartzite,  but  tests  of  lightly  loaded  limestone  below  100°C  almost  inevitably 
involved  smooth  slip.  This  was  attributed  to  absorbed  water  indicated  by 
creep  data  und  rate  process  theory.  Smooth  slip  was  less  readily  obtained 
under  higher  normal  loads,  and  in  the  harder  rocks,  suggesting  penetration  of 
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20.  ABSTRACT  (Continued) 

absorbed,  water  films  by  mineral  contact  asperities.  Because  of  the  strong 
temperature  dependence,  stick-slip  is  attributed  to  asp  verity -to- asperity 
bundi n g . 

Part  XI  presents  results  from  a simple  but;  carefully  controlled  laboratory 
experiment  to  study  the  relation  between  dilatance,  sliding  friction,  and 
resistance  to  deformation  of  densely  packed  assemblages  of  steel  or  Teflon  rods. 
Results  did  not  fellow  the  current  theories  af  alia tan t behavior,  but  suggested 
a simpler  model:  (a;  sliding  first  occurs  along  the  path  of  .lowest  cumulative 

sliding  friction,  (b)  such  sliding  brings  new  surfaces  into  contact,  increasing 
friction  ana  causing  arrested  slip,  (c)  through  repetition  of  this  process  a 
- < .lure  cone  grows  and  continues  to  grow  in  thickness  until  the  number  of 
p t'-ntial  shear  paths  is  depleted  and  catastrophic  failure  becomes  imminent . 
f ' ages  of  primary  secondary  and  tertiary  creep  are  preai.  >»-d  if  time  is 
expressed  as  cumulative  probability.  The  arrested  slip  model,  while  yet 
unverified  in  the  field,  appears  preferable  to  the  usual  assumption  of  cohesive 
cr  clay-like  behavior  of  joints  in  rocks. 
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SUMMARY 

This  Final  Report  consists  of  two  sections  which,  while  involving  dis- 

X 

tinctively  different  research  methods  and  analysis,  are  closely  interrelated 
in  their  significance  for  the  stability  of  rock  masses,  e.g.  rock  slides,  under- 
ground blast  caverns,  etc. 

Part  I deals  with  stick-slip  of  selected  rocks  at  different  temperatures  and 
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under  normal  stresses  up  to  20  kg/ cm  . Stick-slip  was  found  to  dominate  the 
sliding  behavior  of  paired  roughly  ground  surfaces  of  limestone,  basalt,  and 
quartzite.  An  alternating  short-and-long  slip  with  intermediate  creep  oc- 
curred uniquely  with  the  basalt,  attributed  to  slow  dissipation  of  higMy  lo- 
calized high  temperatures  developed  from  sliding.  Tests  of  lightly  loaded 
limestone  below  100°C  almost  inevitably  involved  smooth  slip,  attributed  to 
adsorbed  water,  supported  by  creep  data  and  rate  process  theory.  Smooth 
slip  was  less  readily  obtained  in  limestone  under  higher  normal  loads,  and  in 
the  harder  rocks  under  any  loading,  suggesting  penetration  of  adsorbed  water 
films  by  mineral  contact  asperities. 

The  coefficient  of  sliding  friction  nearly  always  was  initially  low  and  in- 
creased with  sliding  to  approach  a stable  value;  higher  temperatures  increased 
friction  and  caused  longer  sticks  and  larger  slips.  Because  of  this  strong  tem- 
perature dependence,  stick-slip  is  attributed  to  asperity-to-asperjty  bonding. 
The  effect  of  rock  debris  added  between  the  sliding  surfaces  was  in  general  to 
simulate  tests  under  lower  normal  pressures,  and  to  speed  up  sliding  equili- 
bration. 

Part  II  of  this  report  presents  results  from  a simple  but  carefully  con- 
trolled laboratory  experiment  to  study  the  relation  between  dilatancy,  sliding 
friction,  and  resistance  to  deformation  of  densely  packed  assemblages  of 
steel  of  Tef!on  rods.  The  experiment  is  intended  to  model  the  behavior  of 
fractured  rock  masses  deforming  under  load;  but  the  results  also  may  apply 
on  a smaller  scale  to  sliding  of  rock  grains  along  incipient  fractures  of  grain 
contacts.  Both  occurrences  ofdii&f»«oy  have  been  confirmed  in  nature. 
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Shearing  tests  were  conducted  with  controlled  stresses  and  monitored 
movements  and  volume  changes.  Results  did  not  follow  the  current  theories 
of  (tilatant  behavior,  but  suggested  a much  simpler  model  based  on  proba- 
bility: (a)  As  the  ratio  of  principal  stresses  increases,  sliding  occurs  first 
along  the  path  of  lowest  cumulative  sliding  friction,  (b)  Since  sliding  along 
the  path  of  least  resistance  brings  new  surfaces  into  contact,  the  new  con- 
tacts almost  inevitably  have  a higher  cumulative  sliding  friction,  causing 
arrested  slip,  (c)  Sliding  then  will  occur  along  the  next  least  resistant  path, 
and  probably  In  turn  will  become  arrested,  (d)  Through  repetition  of  this 
process  a failure  zone  grows  and  continues  to  grow  in  thickness  until  the 
number  of  potential  shear  paths  is  depleted  and  catastrophic  failure  becomes 
imminent. 

It  is  shown  by  theory  and  by  experiments  that  arrested  slip  of  cohesion- 
less assemblages  can  deliver  an  increasing,  decreasing  or  constant  stress 
ratio,  the  latter  being  characteristic  of  rock  cree»*  A classical  time- 
deformation  curve  showing  stages  of  primary,  secondary  and  tertiary  creep 
is  developed  based  on  an  assumed  equal  likelihood  of  occurrence  of  any 
combination  of  1,  2,  3,  ,.n  slip  planes  where  n is  a discrete  number.  A 
distinction  from  viscous  creep  is  the  requirement  for  an  increasing  shear 
zone  which  might  be  monitored  in  field  situations.  The  arrested  slip  model 
explains  rock  ’ noises"  preliminary  to  failure,  and  appears  referable  to  the 
usual  assumption  of  cohesive  or  clay-like  behavior  of  joints  in  rocks.  The 
model  is  sustained  by  the  experiments  in  Part  I of  the  Report,  in  that  stick- 
slip  dominates  in  widely  varying  rock  types,  and  friction  tends  to  increase 
with  slip  distance. 
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INTRODUCTION 
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Statement  of  the  Problem 

In  naturally  occurring  bodies  of  rock,  movement  of  engineering  significance 
does  not  normally  occur  in  the  body  of  intact  rock,  but  along  joints  and  fracture 
planes.  These  features  may  be  part  of  the  geologic  structure  of  the  rock,  or 
may  be  induced  by  the  engineering  construction  effort  in  the  vicinity  of  the  rock. 
Whether  the  joints  are  caused  by  earth  stresses  or  by  blasting,  stress  relief, 
or  some  other  mechanism,  these  planes  of  weakness  control  the  strength  of  a 
rock  body. 

As  a result  of  engineering  construction  work  such  as  cuts  in  rock  or  under- 
ground exca\  ation,  previously  heavily  loaded  rock  undergoes  a stress  relief  in 
the  direction  of  excavation.  This  may  provide  a sufficient  imbalance  of  forces 
to  cause  movement  along  the  joint  planes.  This  movement  will  be,  initially  at 
least,  at  a very  slow  rate.  The  mode  of  sliding  of  relatively  smooth  rocks  over 
one  another  is  stick-slip.  Th.1  purpose  of  this  investigation  is  to  determine  the 
mechanism  of  the  stick-slip  mode  of  frictional  sliding  in  rocks  with  relatively 
light  normal  stresses. 

Scope 

Stick-slip  movement  occurs  when  loading  systems  are  soft;  that  is,  the 
stiffness  k in  the  equation  F * kx  (F  = frictional  force,  x « distance  moved) 
is  small  (Jaeger,  1971).  Gravity  is  such  a loading  system.  A loading  system 
was  therefore  devised  in  the  laboratory  to  allow  the  production  of  stick-slip  under 
controlled  conditions.  The  study  included  variations  in  several  parameters  which 
were  believed  to  influence  a strong  effect  on  stick-slip.  Particular  emphasis  was 
placed  on  temperature,  since  this  variable  had  not  been  included  in  previous 
studies  by  others.  Other  important  variables  included  are  rock  type  and  effect 
of  debris  or  gouge  created  by  sliding.  Other  variables  were  kept  constant  in 
order  to  isolate  those  believed  to  be  critical.  Relying  on  preliminary  tests  and 
the  experience  of  others,  the  nominal  rate  of  movement  was  held  constant  at  about 
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1.5  x 10  * centimeters  per  minute.  Also  held  constant  was  the  surface  rough- 
ness of  the  three  rock  All  were  ground  with  a number-45  grit  aluminum 

oxide  grinding  wheel. 

The  following  parameters  of  stick-slip  by  use  of  force  and  distance  trans- 
ducers and  oscillograph  recorders:  slip  times,  slip  distances,  shear  load  re- 
laxation upon  slip,  time  between  slips,  and  coefficient  of  static  friction. 
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Sliding  Friction 

A review  of  past  and  current  concepts  of  sliding  friction  seems  necessary 
in  order  to  understand  the  various  concepts  which  have  been  advanced  regarding 
stick-slip  of  rocks.  That  an  effort  is  required  to  move  an  object  along  another 
object,  including  the  earth  itself,  has  been  apparent  ever  since  the  first  sub-human 
tried  it.  For  a long  time  we  have  called  the  force  opposing  such  movement  fric- 
tion. In  the  examination  of  sliding  friction,  some  historical  background  will  be 
briefly  examined.  The  modern  concept  of  sliding  friction  will  then  be  presented, 
as  will  some  alterations  and  modifications  to  the  theory  which  have  been  proposed 
to  handle  objections  caused  by  the  nature  of  rock. 

Historical  background 

A review  of  the  historical  background  of  current  concepts  of  sliding  friction 
up  to  about  1830  is  given  in  Bowden  and  Tabor  (1954,  1964).  They  include  an 
extensive  bibliography  of  the  original  sources. 

Leonardo  daVinci  seems  to  have  been  the  first  to  correctly  observe  the  rela- 
tionship between  friction  and  normal  load.  However,  he  confined  his  experiments 
and  observations  to  his  notebooks.  The  first  man  to  publicly  proclaim  the  basic 
laws  of  friction  was  Guillaume  A montons  (1663-1705),  a French  architect.  In 
a paper  presented  to  the  French  Academy  in  1699,  Amontons  described  the  two 
main  laws  of  friction.  These  are: 

1.  Frictional  force  is  proportional  to  normal  load.  Amontons  said  that 
in  most  cases  it  is  about  one-third  of  the  applied  normal  load. 

2.  Friction  is  independent  of  the  size  of  the  bodies  in  contact.  Amontons 
viewed  surface  irregularities  as  the  cause  of  frictional  force. 

Leonhard  Euler  (1707-1773),  the  Swiss  mathematician  and  physicist,  was  the 
first  to  point  out  that  there  is  a difference  between  static  and  dynamic  (or  kinetic) 
friction. 

One  of  the  great  names  in  the  study  of  friction  is  that  of  Charles  Augustine 
Coulomb  (1736-1806),  a French  engineer.  In  between  his  electrical  experiments 
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and  his  earth  pressure  experiments  and  discoveries  he  also,  in  1781  and  1785, 
published  papers  which  revived  Amontons*  work,  which  had  sunk  into  obscurity. 
Coulomb  considered  that  friction  could  arise  from  lifting  over  asperities,  bending 
of  asperities,  and  breaking  of  asperities.  He  described  the  angle  of  friction  as 
the  slope  of  the  steepest  asperity.  He  also  considered  the  possibility  of  cohe- 
sion, but  rejected  it  because  he  thought  that  cohesion  would  have  to  be  propor- 
tional to  the  gross  contact  area.  Coulomb's  explanation  remains  enshrined  in 
basic  engineering  mechanics  books,  and  it  has  some  passionate  defenders  in  the 
ranks  of  adhesion  physicists  (Bickerman,  1970i. 

During  approximately  the  same  period,  scientists  in  England  were  discovering 
cohesion.  Jean  Theophile  Desaguliers  (1863-1744),  an  English  physicist,  dis- 
covered cohesion  between  clean,  like  solids.  Samuel  Vince  (17  9-1821),  a 
theologian  and  physicist,  considered  the  difference  between  static  and  kinetic 
friction  to  be  due  to  adhesion. 

In  the  late  nineteenth  century,  dissatisfaction  with  the  "up-and-over"  concept 
of  the  mechanism  of  friction  grew  with  the  spread  of  thermodynamic  considera- 
tions of  surfaces.  With  many  contacts,  what  slides  up  must  slide  down,  and 
Coulomb's  concept  would  result  in  a zero  energy  loss.  This  is  clearly  contrary 
to  the  obvious  experience  of  the  creation  of  heat  by  energy  loss  during  frictional 
sliding.  As  a consequence,  a molecular  adhesion  theory  of  friction  began  tc, 
draw  attention  with  the  works  of  Ewing  in  1892,  Sir  Will* am  Hardy  in  1919  and 
in  1936,  and  especially  Tomlinson  in  1929.  Karl  Terzaghi  (1915)  came  to  the  con- 
clusions of  the  modern  theory  independently  and  intuitively.  The  real  develop- 
ment of  the  modern  adhesion  theory  of  friction  came  with  the  work  of  Bowden 
and  Tabor  and  their  associates  from  1939  until  the  present.  It  is  the  Bowden  and 
Tabor  theory  which  is  now  considered  the  "classical"  theory  of  friction  (see 
Byerlee,  1967,  and  Lambe  and  Whitman,  1969). 

The  adhesion  theory  of  friction 

Area  of  contact.  At  least  on  a microscopic  scale,  all  solid  surfaces  are 
rough  and  contain  asperities  (Maurer,  1965),  projections  of  material  farther 
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away  from  an  "average"  surface  location  than  other  points  on  the  surface.  Con- 
tact between  even  smooth  surfaces  is  thus  made  only  at  the  contacts  of  the  asperi- 
ties. Obviously,  true  contact  area  bears  very  little  relationship  to  apparent 
contact  area.  How  can  true  contact  area  be  calculated?  Terzaghi  (1925),  sug- 
gested that  the  area  of  true  contact  would  be  related  to  the  normal  load  and  to 
the  yield  stress  of  the  material.  Some  measurements  of  the  true  contact  area 
of  metals  have  been  made  by  electrical  and  optical  means.  The  results  of  these 
measurements  seem  to  confirm  the  following  theoretical  analysis  (Bowden  and 
Tabor,  1954): 

By  the  Hertzian  elastic  analysis  the  area  of  contact  A of  a solid  deforming 
elastically  is: 

A = kW2/3  (1) 

2/3 

i.e. , proportional  to  W ‘ where  W is  the  normal  load.  As  the  load  W is 
increased,  the  mean  pressure  Pm  increases  until  a critical  point  is  reached 
within  the  softer  of  the  materials  in  contact.  When  the  mean  pressure  has  be- 
come high  enough  to  cause  the  regions  in  contact  to  flow  plastically,  the  mean 
pressure  is  directly  proportional  to  the  elastic  limit  Y 


P 


m 


= CY. 


(2) 


Mean  pressure  is  therefore  independent  of  the  size  of  the  asperities,  and  so 
is  independent  of  the  load.  This  makes  pm  the  yield  pressure  p.  At  the 
point  below  the  surface  where  p^  = 1.1Y,  plasticity  is  initiated  at  the  most 
stressed  point.  The  corresponding  load  is: 


W = 


(3) 


for  a spherical  asperity  pressing  on  a flat  plate.  The  radius  of  the  sphere 
is  denoted  by  r,  the  elastic  moduli  of  the  sphere  and  the  plate  by  Ej  and  Eg. 
If  the  asperity  is  of  a conical  or  pyramidal  shape,  the  tip  can  be  considered 
as  a sphere  with  a very  small  radius.  Therefore  a small  load  is  sufficient  to 
deform  the  asperities  beyond  the  elastic  limit.  Thus,  at  least  for  more  ductile 
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materials,  the  tips  easily  deform  plastically. 

Since  the  yield  pressure  is  independent  of  the  amount  of  deformation  which 
has  already  occurred,  the  area  of  real  contact  is  directly  proportional  to  the 
load  W.  That  is: 

A » - W.  (4* 

r p ' > 

Equations  (1)  and  (4)  show  that  viiile  elastic  deformation  is  occurring  on  a 

2/3 

macroscopic  scale,  and  the  apparent  area  of  contact  is  increasing  with  W , 

plastic  deformation  is  occurring  in  the  asperities,  and  the  area  of  real  contact 

increases  with  W.  This  has  been  demonstrated  experimentally  in  an  indirect 

manner  by  the  measurement  of  the  conductivity  of  metals. 

Archard  (1961)  points  out  that  most  asperity-asperity  contacts  seem  to  occur 

with  no  visible  damage.  He  concludes  that,  while  a single  contact  would  deform 
2 /3 

with  A = kVV  , as  the  number  of  contacts  increases  A approaches  kW. 

Thus,  many  asperities  may  deform  elastically  under  some  circumstances. 
Archard  earlier  concluded  (1957)  that,  in  metal,  initial  contacts  are  plastic, 
but  repeated  contacts  are  elastic.  He  deduced  that  the  average  size  of  the 
real  area  of  contact  of  asperities  is  a constant,  but  that  the  number  of  contacts 
increases  with  an  increase  in  normal  load.  This  may  be  true  especially  under 
relatively  low  normal  loads.  Still,  plastic  deformation  of  asperities  seems  to 
offer  the  best  general  explanation  for  frictional  phenomena,  as  will  be  shown 
below. 

Upon  application  of  a tangential  force  in  addition  to  the  normal  force,  the 
area  of  real  contact,  at  least  of  ductile  materials,  increases  (Merchant,  1968). 
By  analogy  with  the  von  Mises  yield  criteria  for  an  ideal  two-dimensional  model 
of  an  intermetallic  junction,  the  yield  criteria  may  have  the  form: 

2 2 2 

P + OS  « PQ  (5) 

where  pQ  is  static  contact  pressure,  s is  tangential  force,  and  a is  a 
constant  with  a value  of  about  10.  The  material  in  the  asperity  is  initially  sub- 
jected to  a normal  stress  equal  to  its  yield  pressure  p = p^.  Thus  as  soon  as 
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the  slightest  tangential  stress  is  applied,  equation  (5)  is  satisfied  only  if 
p diminishes.  This  means  that  the  area  of  contact  must  increase  as  p de- 
creases, following  equation  (4).  This  process  might  go  on  indefinitely  if  it 
were  not  for  the  presence  of  surface  contaminants,  whioh  stop  junction  growth 
and  permit  slip.  Tangential  stress  may  cause  asperities  which  deform  elas- 
tically under  normal  load  to  deform  plastically.  Tangential  stresses  in  brit- 
tle materials  may  generate  fracture  in  or  below  the  sliding  interface  (Tabor, 
1908). 

Once  sliding  has  begun,  the  temperature  of  the  area  in  real  contact  increases 
greatly  (Bowden  and  Tabor,  1954).  This  allows  the  areas  of  real  contact  to  grow 
more  rapidly  than  indicated  by  movement  of  the  gross  block. 

Adhesion  theory  . Two  effects  seem  to  predominate  when  "smooth"  sur- 
faces are  slid  across  one  another.  One,  the  equivalent  of  interlocking  of  the 
Coulomb  theory,  causes  plowing  in  metals  and  may  cause  either  brittle  fracture 
or  dilatancy  in  non-metals.  This  plowing  term  is  estimated  to  cause  a "friction 
coefficient"  of  about  0.1  (Merchant,  1968).  The  other  effect,  which  causes 
what  early  experimenters  called  "true  friction"  though  their  explanations  were 
exactly  opposite  to  today’s,  is  called  "shearing"  by  Bowden  and  Tabor  (1954). 
This  is  the  action  of  one  asperity  lodged  against  another  in  the  direction  of  shear. 

Experiments  and  photomicrographs  by  Bowden  and  Tabor  show  that  the 
shearing  component  of  friction  is  caused  by  actual  "cold  welding"  of  metals  to 
each  other.  This  bonding  or  rather  the  force  required  to  shear  the  bonding 
comprises  sliding  friction. 

It  is  clear  that  when  similar  metals  are  in  contact  the  local  high 
pressure  must  cause  .an  equal  flow  of  both  surfaces  as  that  they  both 
contribute  equally  to  the  formation  of  welded  junctions  and  when  sliding 
takes  place,  both  surfaces  will  be  distorted  and  torn.  Moreover,  since 
the  metal  will  be  hardened  at  the  interface,  the  junction  will  be 
stronger  than  the  underlying  metal,  so  that  the  break  will  occur  in  the 
bulk  of  the  metal  (Bowden  and  Tabor,  1954). 

Similar  effects  are  found  in  non-metals  (Bowden  and  Tabor,  1964).  The 
breaking  point  is  different  in  materials  other  than  metals,  and  at  very  small 
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loads  the  experimental  results  were  similar  to  those  at  high  loads. 

In  metals  the  actual  bonding  is  believed  to  be  of  the  same  nature  as  al- 
loying, although  results  by  Johnson  and  Keller  (1967)  throw  doubt  upon  this  con- 
clusion. When  two  asperities  come  close  together,  the  appvoaching  bodies  have 
equal  effects  on  the  tip  atoms.  Thus,  the  forces  at  the  interface  are  the  same 
as  the  forces  in  bulk.  Because  of  mismatch  between  crystal  lattices  it  would 
be  the  greatest  of  coincidences  if  the  adhesion  were  as  strong  as  in  the  bulk 
metal  (Merchant,  1968),  The  same  general  interaction  occurs  with  dissimi- 
lar meials,  except  that  the  junction  should  be  much  weaker.  Some  non-brittle 
non-metallic  crystalline  substances  such  as  NH^Cl  and  sulphur  seem  to  be- 
have in  the  same  manner  as  metals.  Even  some  brittle  materials  such  as  rock 
salt  and,  under  high  pressure,  sapphire,  seem  to  show  adhesion  between  as- 
perities (Bowden  and  Tabor,  1964).  Contamination  generally  present  on  any 
solid  body  in  air  will  act  to  reduce  bonding  (Merchant,  1968). 

Laminar  solids  apparently  adhere  to  each  other  due  to  Van  der  Waals 
forces.  In  experiments  on  mica  (muse mite)  cleaved  to  present  molecularly 
smooth  surfaces  on  cleavage  faces,  Bowden  and  Tabor  (1964)  found  that  the 
attractive  energy  of  the  flat,  plateau-like  asperities  was  quite  high  over  a re- 
latively small  separation  distance  when  chemically  clean. 

Amontons'  laws  by  the  adhesion  theory.  If  we  consider  only  hard  mater- 
ials that  make  contact  only  at  the  tips  of  the  asperities,  i.e, , a polished  sur- 
face with  no  interlocking,  the  "plowing"  term  can  be  ignored.  In  this  case, 
the  real  area  of  contact  is,  according  to  equation  (4): 


If  s is  the  mean  tangential  stress  required  to  shear  the  junction,  then  the 
shearing  force  (the  friction  force)  is; 

F - As  (6) 

Since  the  real  area  of  contact  A^  depends  only  on  W and  p then: 
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(7) 


F = W-  , 

P 

and  thus  F is  proportional  to  the  normal  load.  This  is  Amontons'  (Coulombs) 
first  law.  Since  frictional  force  F is  directly  proportional  to  load  W,  the 
proportion,  the  coefficient  of  friction  pi , is  independent  of  both  load  W and 
of  apparent  area  of  contact.  In  fact,  the  coefficient  of  friction 

F s 

<8> 

Since  s and  p are  material  properties,  they  do  not  depend  upon  the  size  of  the 
body.  This  is  Amontons’  second  law. 

Bowden  and  Tabor  (1954)  emphasize  that  intimacy  of  contact  and  strength 
of  adhesion  at  points  of  contact  will  be  greatly  influenced  by  the  presence  of 
surface  films.  With  metals,  these  films  are  usually  oxides  of  the  metal;  with 
non-metals,  the  films  are  adsorbed  fluids  and  gases.  The  contaminant  mole- 
cules will  also  touch  and  will  develop  a shear  strength  s0  over  a real  contact 
area  of  A(1  - c),  if  o is  the  fraction  of  direct  non-contaminated  contact  area. 
The  friction  force  will  therefore  be; 

F = A(ps1  + (1  -o)  s2).  (9) 

Since  c and  s2  are  variables,  yu  is  somewhat  variable  even  for  similar 
materials  with  identical  surface  geometry. 

Removal  of  load  . If  true  adhesion  occurs  across  boundaries  of  bodies 
in  contact,  why  do  the  bodies  not  remain  adhered  when  the  normal  load  is  re- 
moved ? Terzaghi  (1925)  hypothesized  that  the  bonds  would  be  broken  by  elastic 
rebound  of  the  macroscopic  portion  of  the  deformation.  Bowden  and  Tabs**  agree 
that  elastic  rebound  of  asperities  and  macroscopic  portions  of  bodies  in  contact 
are  sufficient  to  break  the  bonds  formed  (1954). 

Johnson  and  Keller  (1967)  discovered  no  evidence  of  elastic  rebound  in  soft 
metals  under  ultra-high  vacuum.  With  contaminated  (real)  surfaces  and  with 
harder  materials,  rupture  of  the  always -present  contaminant  layer  is  necessary 
for  good  bonding.  Bonding  is  therefore  weaker,  and  rupture  of  bonds  by  elastic 
relief  seems  possible. 


Summary#  To  summarize  the  Bowden  and  Tabor  theory  of  friction  in  a 
few  words: 

The  experiments  show  that,  under  the  intense  pressure  which 
acts  at  the  summits  of  surface  irregularities,  a localized  adhesion 
and  welding  together  of  the  metal  surfaces  occurs.  When  sliding 
takes  place,  work  is  required  to  shear  the  welded  junctions  and  also 
to  plough  out  the  metal  (Bowden  and  Tabor,  1954). 

This  theory  has  been  extended  to  non-metals  and  is  now  the  generally  acceptor 
"classic"  theory  of  friction.  Under  this  definition  of  friction,  studies  such  as 
those  of  Patton  (1966)  and  of  Ripley  and  Lee  (1961)  which  concentrate  on  macro- 
irregularities along  joint  planes  are  really  studies  of  dilatancy  and  shear  strength 
and  not  of  friction. 

Friction  of  Rock 

Compared  to  the  amount,  of  study  done  on  metals,  the  amount  of  research  on 
friction  of  rocks  and  minerals  is  very  small.  Most  of  the  pertinent  literature 
(Jaeger,  1959,  1971)  concerns  itself  with  the  approach  of  rock  to  Coulomb's 
shear  strength  criteria: 

r " yttcrn  •*  c.  (10) 

It  must  be  kept  in  mind  in  this  section  that  friction  along  separated  planes  and 
joints  in  rock  is  the  topic  of  discussion.  The  "internal  friction"  of  the  intact 
rock  body  is  not  considered. 

Horn  and  Deere  (1962)  experimented  on  the  coefficient  of  friction  of  indivi- 
dual minerals.  They  found  that  the  static  coefficient  of  friction  depended  upon 
whether  the  minerals  were  massive  crystalline  (quartz,  feldspars,  calcite)  or 
laminar  (micas,  talc,  clay  minerals).  They  found  that  up  to  a critical  roughness 
the  coefficient  of  friction  of  massive  minerals  depends  upon  the  amount  of  water 
present,  water  acting  as  an  anti-lubricant.  The  laminar  minerals,  on  the  other 
hand,  were  lubricated  by  water  in  all  cases.  Water  seemingly  separated  the 
plates  and  reduced  Van  der  Waals  forces  from  plate  to  plate.  Coulson  (1970) 
has  done  similar  work  on  rock  systems,  investigating  coefficient  of  friction  or 
rocks  with  different  grinds  and  polishes,  and  moisture  conditions. 
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Rock  is  normally  taken  to  be  a brittle  material.  With  rock;  along  with 
other  brittle  materials,  objection  has  been  taken  to  the  applicability  of  the  ad- 
hesion theory  cf  friction.  Greenwood  and  Williamson  (1966)  concluded  that  the 
behavior  of  a surface  was  absolutely  determined  by  a plasticity  index: 

9 

where  K*  = IjE/fl-vT)  for  similar  materials  and  H = hardness,  <7=  standard 
deviation  ■?'  height  of  asperity,  and  B = radius  of  the  asperity  tip.  Ircrease 
of  load  would  not  make  a material  with'f  loss  than  1 behave  plastically  at  any 
point. 

Jaeger  (1971)  suggested  that  the  mechanism  of  Bowden  and  Tabor  and  of 
Archard  is  not  applicable  to  brittle  materials.  He  suggests  use  of  the  Coulomb 
criteria,  or  a modification  thereof  in  the  form 


, ni 

T - T + k a 
o n 


(12) 


with  k and  m experimentally  determined  constants.  This  criteria  ha*-  been 
advocated  by  Maurer  (1965)  and  others.  Jaeger  does  not  suggest  the  source  of  the 
cohesion  (c)  in  the  basic  Coulomb  equation,  but  he  states  that  it  ranges  from 
0 to  200  pounds  per  square  inch  in  natural  joints. 

Patton  (1966),  Ripley  and  Lee  (1961),  Coulson  (1970),  and  other  investigators 
have  concentrated  on  the  effect  of  macro-irregularities  with  the  resulting  equation 
(in  the  absence  of  cohesion)  of 

T = tan  (0  + 0),  (13) 


where  0 is  the  inclination  of  a slope  in  the  direction  of  r. 

Some  experimenters,  especially  Byerlee  (1967,  1969)  and  Brace  and  Byerlee 
(1970)  have  contended  that  the  mechanism  of  friction  in  brittle  materials  such  as 
rock  differs  from  that  of  the  Bowden  and  Tabor  adhesion  theory.  Byerlee  believes 
that  the  experimental  evidence  showing  cracking  of  quartz  and  granite  in  indenta- 
tion hardness  tests  proves  that  with  such  i j tens  ' strength  materials,  plastic 
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flow  at  the  asperities  is  not  possible.  Byer?ee  agrees  with  the  classic  theory 
in  postulating  asperities  of  basically  conical  or  wedge  shape.  His  conclusion 
is  that  failure  is  more  likely  to  be  by  brittle  fracture  due  to  the  low  tensile 
strength  of  the  material. 

Byerlee  (19(57)  assumes  normal  and  tangential  stresses  at  the  tip  of  a conical 
asperity  of  vanishingly  small  tip  radius  . The  asperities  break  when  the  shear- 
induced  tensile  stress  along  the  outer  edge  of  the  asperity  equals  the  tensile 
strength.  He  can  calculate  a theoretical  initial  coefficient  of  friction  by  cal- 
culating the  strength  of  an  asperity'  on  its  outer  edge  as  it  crushes.  The  calcu- 
lated theoretical  coefficient  of  friction  is  0.10.  The  predicted  value  of  coefficient 
of  friction  for  massive  crystalline  materials  is  quite  similar  to  that  found  ex- 
perimentally by  Horn  and  Deere  (19(52)  for  a clean,  dry  surface.  Experiment 
also  showed  no  variation  of  coefficient  of  friction  with  apparent  area  of 
contact.  The  coefficient  of  friction  does  vary  with  gross  roughness,  and 
Byerlee  assumes  that  with  rough,  surfaces  the  theory  will  not  apply  due  to 
interlocking  of  asperities  predominating  over  asperity  to  asperity  contact. 

The  greatest  objection  to  Byerlee’s  assumption  of  independence  of 
this  theory  from  the  adhesion  theory’  is  the  tangential  shear  forces  assumed 
present  at  the  tips  of  asperities.  No  cause  is  given  for  the  existence  of 
a friction  or  shear  force  between  asperity’  tips.  Even  with  brittle  material, 
the  adhesion  theory  of  friction  cannot  be  excluded.  Tabor  (1968)  considers 
that  the  shearing  of  asperity  to  asperity  contacts  in  brittle  materials  may  occur 
by  tensile  failure  in  the  bulk  material  below  the  contact,  but  that  adhesion 
will  exist  at  the  contact. 

Stick-Slip 

Sti  ck-slip  is  the  movement  of  a sliding  element  across  a stationary  element 
in  abrupt,  rapid  movements  alternating  with  periods  of  no  movement.  The  applied 
tangential  force  seems  to  build  up  to  a critical  value  required  to  cause  movement, 
then  to  suddenly  relax  with  an  accompanying  forward  movement  of  the  sliding 
element.  The  sliding  element  then  stops  moving  (sticks)  and  the  cycle  begins 
anew.  The  variation  of  frictional  force  with  time  is  shown  in  Figure  1,  a tracing 
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of  an  actual  test.  The  velocity  of  the  sliding  element  fluctuates  greatly  during 
movement. 

Stick-slip  is  a phenomenon  which  has  been  observed  in  most  friction  studies, 
particularly  at  light  normal  stresses.  Such  self-induced  vibrations  were  recog- 
nized in  the  nineteenth  century  as  the  source  of  music  from  violin  bows  and 
strings , and  of  the  squeak  of  chalk  on  the  blackboard.  Some  noises  in  natural 
materials  such  as  noises  in  rocks  before  apparent  breakage,  and  the  squeakiness 
when  walking  across  some  sands  may  be  the  result  of  stick-slip.  Wells  (1929) 
was  the  first  experimenter  to  notice  the  phenomenon  scientifically.  Most  of  the 
research  on  stick-slip  has  been  with  metals,  where  stick-slip  at  high  velocity 
of  slider  movement  is  a troublesome  phenomenon  in  the  machine  tool  industry 
(Brockly  et.  al.  19G7,  Derjaguin  et.  al.  1956).  Much  of  this  work  has  been  of  a 
mathematical  nature,  seeking  to  solve  the  general  equation: 

d2x  . dx  dx  . , , . ..  . ... 

mdF  + AdT  * F(dT’  X»  V + k(x‘vt)  = 0 <14> 

m — A— 

where  dt^  is  the  inertial  force  of  an  accelerating  slider  of  mass  m,  dt 

is  the  damping  factor,  F(dx/dt,  x,  tQ)  is  the  frictional  force  which  is  a function 
of  velocity,  time  of  contact  between  surfaces,  and  relative  displacement,  and 
k(x  - vt)  is  the  driving  or  retarding  force,  with  k the  system  stiffness.  At- 
tempts are  usually  made  to  solve  the  equation  in  terms  of  the  natural  frequency 
of  the  sliding  system  or  to  determine  the  critical  velocity  or  the  critical  stiff- 
ness or  other  conditions  under  which  stick-slip  can  be  eliminated.  An  emphasis 
has  been  placed  on  the  discovery  that  stick-slip  can  be  eliminated  by  making 
the  system  stiffer.  Gravity-driven  systems  under  light  normal  loads,  such  as 
rock  joints  in  open  cuts,  are  a soft  system  (Jaeger,  1971). 

Less  work  has  been  done  in  attempting  to  define  mechanisms  of  stick-slip. 
Outstanding  has  been  the  work  of  Rabinowicz  (1965).  He  states  "All  stick-slip 
processes  are  caused  by  the  fact  that  the  junction  force  does  not  remain  con- 
stant as  a function  of  3ome  other  variable, " Irregular  and  harmonic  stick- 
slip  vibrations  can  occur  in  machine  systems.  In  natural  systems,  self-induced 
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vibrations,  or  regular  stick-slip,  is  the  common  type  found. 

Bowden  and  Tabor  (1954)  note  that  if  static  friction  is  higher  than  kinetic 
friction,  and  one  of  the  sliding  surfaces  has  a certain  degree  of  elastic  free- 
dom, the  motion  will  be  intermittant  and  characterized  by  stick-slip.  An  ana- 
lysis of  the  equations  of  motion  show  that  stick-slip  under  these  conditions  must 
occur  when  the  static  coefficient  of  friction  yug  is  greater  than  the  kinetic  co- 
efficient of  friction  A’  and  that  the  magnitade  of  stick-slip  reduces  as  yu^ 
approaches  yug.  Under  stick -slip  conditions,  yu^  is  determined  by  the  shear 
load  at  one-half  the  slip  distance.  Tolstoi  (1967)  showed  that  stick-slip  is  depen- 
dent on  the  degree  of  freedom  normal  to  the  surface  of  sliding.  When  dilatancy 
is  not  possible,  stick-slip  disappears.  This  might  be  significant  under  geologic 
conditions. 

Why  should  the  static  coefficient  of  friction  be  larger  than  the  kinetic  coef- 
ficient? Rabinowicz  (1951),  following  the  initial  suggestion  of  Ishlinski,  proposed 
that  the  static  coefficient  of  friction  varies  with  the  time  of  contact  between  the 
sliding  elements,  building  up  exponentially  within  a short  time  of  contact  to  be 
higher  than  kinetic  friction.  He  suggested  that  junction  growth  was  the  producer 
of  high  values  of  static  friction.  He  equated  junction  growth  with  creep  in  a ver- 
tical direction,  and  suggested  that  yug  is  controlled  by  (1)  the  creep  rate  in 
compression  of  asperities  resulting  in  increase  of  junction  area  and  (2)  shear 
strength  of  junctions  formed.  As  suggested  previously,  contaminants  cause 
junction  growth  to  cease  and  tangential  stress  to  build  up  faster  than  shear 
strength.  The  concept  of  time  dependency  of  static  friction  has  been  supported 
by  the  research  of  Kosterin  and  Kraghelsky  (1952)  who  derived  an  equation  relating 
duration  of  contact  to  stress  and  strain  in  the  contacted  surfaces  and  thus  to  fric- 
tion. Static  friction  was  experimentally  determined  to  follow  the  equation  in  metals 
ami  plastics.  Brockley  and  Davis  (1968)  used  controlled  stick-slip  to  investigate 
postulated  time  dependence  of  static  friction  in  metals.  They  found  that  the 
static  friction  coefficient  does  increase  with  time  of  stick  if  there  is  a difference 
between  and  yu^.  Simkins  (1967)  and  Tolstoi  (1967)  deny  that  there  is  a 
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difference  between  static  and  kinetic  coefficients  of  friction.  They  claim  that 
the  apparent  differences  are  caused  by  observing  the  statistical  maximum  and 
minimum  of  the  same  physical  property. 

A result  of  the  equations  of  motion  is  that  with  a small  average  velocity 
of  sliding  (total  distance  moved/total  time),  the  maximum  velocity  achieved 
during  slips  may  be  quite  large  (Derjaguin,  1956V  This  may  influence  the  sur- 
face properties,  especially  temperature.  The  temperature  rise  upon  slip  is 
more  marked  with  poor  thermal  conductors  (Bowden  and  Tabor,  1954). 

Most  of  the  above-cited  work  has  concentrated  on  the  elimination  of  stick- 
slip  as  a nuisance.  Brockley,  Cameron,  and  Potter  (1967)  determined  that 
stick-slip  will  be  eliminated  by: 

1.  A small  load  and  a large  stiffness. 

2.  A small  difference  in  values  of  ^ and 

3.  A large  damping  ratio. 

4.  A small  value  of  static  friction  growth  constant  (junction  growth). 

Stick-Slip  of  Rock 

Little  previous  work  has  been  done  on  stick-slip  as  applied  to  rocks.  Jaeger, 
in  his  pioneering  study  of  the  frictional  properties  of  rocks  (1959),  encountered 
violent  stick-slip  on  bare,  flat  surfaces  of  plaster.  Small  plaster  fragments, 
frequently  arranged  in  slickensides,  were  present  after  slip  at  high  stresses. 

Jaeger  speculated  from  the  high  calculated  energy  release, that  high  temperature 
at  points  of  contact  had  occurred  during  stick-slip. 

Korn  and  Deere  (1962)  in  their  study  of  friction  of  specimens  of  pure  minerals 
encountered  stick-slip  in  many  of  their  tests  on  quartz.  To  them,  as  to  most 
other  investigators  of  the  friction  of  rocks,  the  phenomenon  was  a nuisance  in 
attempts  to  obtain  static  friction  coefficients. 

Coulson  (1970)  in  a study  of  friction  of  rocks  with  varying  surface  roughness 
found  that  stick-slip  occurred  within  a wide  range  of  normal  stresses  in  several 
rock  types  in  a direct  shear  type  test.  Stick-slip  was  most  common  in  smooth- 
ground  specimens  (number-600  grit)  at  high  normal  stresses.  Coulson  experienced 
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stick-slip  principally  in  rocks  containing  quartz,  but  also  found  it  in  a basalt. 
Interestingly,  he  did  not  obtain  stick-slip  in  tests  on  the  same  Indiana  limestone 
used  in  part  of  this  study.  He  did  encounter  conditions  reflecting  that  the  static 
coefficient  of  friction  is  a function  of  duration  of  stick.  Some  effort  was  made 
by  Coulson  to  eliminate  stick-slip  by  increasing  the  stiffness  of  his  loading 
system. 

Dieterich  (1970)  conducted  experiments  to  determine  the  time  dependence 
of  s.atic  friction  in  stick-slip  sliding.  In  a series  of  experiments  on  sandstone, 
graywacke,  quartzite,  and  granite  with  normal  stress  of  from  0.02  to  0.85 
kilobars  (1  bar  = 14. 5 psi  = 1.  02  kg/cm2),  he  found  that  sliding  on  clean  rough- 
ground  surfaces  was  initially  stable.  However,  as  rock  debris  accumulated  on 
the  surfaces,  stick -slip  became  the  dominant  mode  of  sliding.  Dieterich  found 
that  change  of  the  interval  between  slips  from  15  seconds  to  24  hours  had  about 
a 10%  influence  on  yu.g.  He  related  stick-slip  to  periodicity  of  shallow  earth- 
quakes. 

An  investigation  into  friction  of  rocks  by  Hoskins  et.  al.  (1968)  included 
some  specific  study  of  stick -slip.  These  investigators  determined  that  stick- 
slip  occurred  with  smooth  surfaces  tested  (35  + 5 microinches  roughness),  and 
that  rougher  blocks  demonstrated  a smooth  increase  of  frictional  force  up  to  a 
value  controlled  by  the  normal  stress.  Hoskins  et.  al.  observed  a small  dis- 
placement occurring  before  the  main  slip.  The  time  length  of  slip  was  very 
abrupt,  usually  less  than  2 milliseconds.  They  concluded  that  the  stick-slip 
observed  fit  the  mathematical  pattern  of  even  oscillation  described  by  Derjaguin 
et.  al.  (1956)  and  other  mathematical  investigators  quite  well.  They  made  no 
comments  relative  to  the  mechanism  - f '■tick-slip. 

In  experiments  on  artificial  joint  systems  made  of  blocks  of  precisely 
smooth-ground  plaster,  Brown  and  Trollope  (1970)  discovered  that  sliding  along 
joints  normally  occurs  in  a stick-slip  mode.  With  joints  oriented  at  45°  to 
principal  stress  directions,  stick-slip  occurred  regardless  of  confining  pressure. 
At  ether  orientations,  stick-slip  occurred  at  high  pressures  (over  200  psi). 
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The  only  investigation  to  emphasize  mechanisms  of  stick-slip  has  been  that 
of  Byerlee  (1970),  who  concentrated  on  stick-slip  which  can  occur  at  very  high 
pressure  deep  in  the  earth.  His  investigations  led  him  to  propose  a theory  of 
stick-slip  of  rock  which  is  dependent  on  the  mechanism  of  brittle  fracture. 
Basically,  stick-slip  is  seen  as  a process  whereby  the  irregularities  on  the 
surface  become  locked,  fail  by  brittle  fracture,  and  then  lock  again.  When  suffi- 
cient fault  gouge  builds  up,  the  sliding  becomes  stable  at  low  confining  pressures, 
but  remains  unstable  at  high  pressure  due  to  the  dense  packing  of  the  debris 
grains.  Failure  of  one  grain  then  leads  to  a high  stress  in  adjoining  grains,  and 
the  system  becomes  unstable  again.  This  system  has  been  seen  as  the  explana- 
tion for  shallow  focus  earthquakes  (Brace  and  Byerlee,  1970). 

In  his  recent  survey  of  friction  of  rock,  Jaeger  (1971)  st  ■'  js  that  stick-slip 
is  a complex  phenomenon  related  to  surface  roughness,  average  rate  of  sliding 
of  the  rocks,  time  between  slips,  wetness  of  the  surface,  and  the  amount  of 
gouge  present.  The  present  study  adds  temperature  as  an  additional  variable 
pertinent  to  the  mechanism  for  stick-slip  of  rock. 

NORMAL  LOAD 


Fig.  2.  Schematic  of  shear  apparatus:  (1)  stationary  block  and  holder; 

(2)  sliding  block  and  holder. 
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TEST  METHODS  AND  PROCEDURES 
Apparatus 

Direct  shear  testing  was  performed  on  prepared,  separated  rock  specimens 
to  determine  the  mechanism  of  stick-slip.  Hoskins  et.  al.  (1968)  presented 
objections  to  the  use  of  pre-cut  triaxial  specimens  for  low  normal-stress  fric- 
tion investigations.  Lajtai  (1969)  has  shown  that  the  direct  shear  mode  of  test 
most  closely  resembles  the  actual  condition  of  state  of  stress  along  a fault  or 
joint  after  the  faulting  has  taken  place.  Gravity  is  a low-stiffness  mechanism 
(Jaeger,  1971)  and  the  direct  shear  mode  most  closely  resembles  the  gravity- 
induced  case  of  movement  along  a joint. 

The  equipment  used  to  conduct  the  stick-slip  tests  is  a temperature  controlled, 
constant-rate-of-strain  direct  shear  device.  Figure  2 is  a schematic  of  the  device, 
omitting  the  temperature  control  features.  The  direct  shear  mechanism  is  com- 
prised of  a 15-cm  by  10-cm  stationary  rectangular  specimen  holder,  with  a smaller 
rectangular  moveable  holder  measuring  internally  8. 3 cm  by  7. 5 cm  placed  on 
top.  An  insert  into  the  moveable  holder  will  hold  smaller,  2 cm  by  2 cm,  rock  speci 
mens.  The  rock  specimens  are  placed  in  these  holders  and  held  in  by  set  screws 
against  metal  plates.  Approximately  one  inch  of  travel  is  possible  with  no  change 
in  apparent  contact  area.  A normal  load  is  transferred  to  the  upper  block  by  a 
ball  and  socket,  which  allows  a free -hanging  load  to  travel  with  the  moving  block. 

The  upper  block  is  pulled  across  the  lower  block  at  a constant  nominal  rate  of 

strain  by  a screw  and  wedge  mechanism  driven  by  an  electric  motor.  The  screw 

advances  the  wedge  at  a very  slow  rate.  A pulley  impinging  on  the  wedge  is  pulled 

downwards  along  the  wedge.  This  causes  the  upper  block  to  move  by  means  of  an 

-3 

attached  steel  wire.  The  upper  block  moves  at  a nominal  rate  of  1. 5 x 10 
centimeters  per  minute. 

The  shear  box  is  enclosed  in  a double-walled  asbestos  oven.  Temperature  is 
controllable  within  the  oven  from  room  temperature  to  275°  Centigrade  by  an 
Athena  Model  63  temperature  controller.  The  controller  is  connected  to  a 
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thermocouple  placed  through  the  lid  of  the  oven,  and  to  electric  heating  coils 
within  the  oven  which  activate  when  the  temperature  falls  below  the  set  value. 

Air  temperature  in  the  oven  is  controllable  to  within  about  2°C.  Closer  con- 
trol was  difficult  due  to  slots  cut  in  the  base  of  the  oven  to  allow  the  normal 
load  to  move  with  the  moving  specimen  holder.  A fan  was  also  provided  inside 
the  oven,  but  seldom  used. 

Instrumentation 

The  shear  load  required  to  move  the  block  is  measured  by  a strain  gage 
direct-load  transducer  mounted  as  part  of  the  pulling  mechanism.  Two  arms 
of  a wheatscone  bridge  are  mounted  so  as  to  be  additive  on  a thin  strip  of 
metal.  As  the  metal  strains,  the  load  is  measured  by  previous  direct  cali- 
bration on  an  oscillograph.  Three  different  transducers  were  used  during  the 
tests,  but  all  limestone  tests  with  normal  stress  of  3 kg/cm2  and  over  and  all 
tests  on  the  other  rocks  were  performed  using  the  same  copper  transducer. 

This  transducer  is  accurate  to  within  100  grams.  The  transducers  used  on  the 
lighter  normal  stress  tests  of  the  limestone  are  accurate  to  within  300  grams. 
Temperature  compensation  is  provided  by  mounting  dummy  gages  on  material 
identical  to  that  of  the  transducer. 

Horizontal  movement  of  the  upper  block  was  measured  by  a linear  variable 
differential  transformer  {LVDT)  Schaevitz  Model  1000HR.  The  core  of  the 
transformer  is  attached  by  a metal  rod  directly  to  the  moveable  upper  shear  box. 
The  LVDT  body  is  Teflon  lined  to  insure  smooth  sliding,  and  the  LVDT  is  accurate 
to  within  0. 0005  inches  (0. 0013  cm).  The  LVDT  had  to  be  reset  by  hand  in  its 
holder  whenever  the  total  movement  exceeded  0.  02  inches  (0. 0508  cm)  at  the 
normal  attenuation  used.  Calibration  was  by  a micrometer  accurate  to  0. 0001 
inches.  A second  LVDT  was  arranged  to  measure  vertical  movement  of  the 
upper  block.  Though  this  LVDT  was  accurate  to  within  0. 0002  inches,  the 
blocks  were  ground  too  smooth  for  any  vertical  movement  to  be  recorded,  and 
use  of  this  LVDT  was  discontinued  later  in  the  program. 
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The  horizontal  LVDT  signal  was  conditioned  by  a Schaevitz  Model  SCM 
signal  conditioner  and  further  amplified  by  a Brush  Model  RD  4215-10  DC 
amplifier.  The  load  transducer  was  activated  and  amplified  by  a Brush 
Model  4212  carrier  amplifier.  The  information  gathered  was  recorded  on 
a Brush  Mark  200  ink-writing  oscillograph.  The  char'  was  usually  run  at 
speeds  of  20  or  50  seconds  per  centimeter,  though  faster  rates  were  some- 
times used.  Times,  movements,  and  loads  were  read  directly  from  the 
oscillograph  record.  Temperature  was  read  from  the  set-dial  and  deviation 
dial  of  the  temperature  controller. 

Specimen  Preparation 

Each  block  was  saw-cut  to  the  appropriate  size,  then  ground  "smooth"  with 
a number-45  grit  aluminum  oxide  grinding  wheel  on  an  automatic  grinder. 

This  degree  of  smoothness  is  much  rougher  than  the  "rough"  quartz  block  of 
Horn  and  Deere  (1962),  which  was  ground  to  240  grit;  slightly  rougher  than 
Coulson's  (1970)  rough  blocks  (number-80);  but  smoother  than  the  rough  blocks 
of  Hoskins  et.  al.  (1968)  and  the  natural  fracture  surfaces  tested  by  Byerlee 
(1970).  The  grinding  was  believed  necessary  to  minimize  dilatancy  (inter- 
locking) and  to  give  some  degree  of  uniformity  to  the  block  surface  structure 
so  that  tests  run  on  different  blocks  could  be  compared.  The  number-45  grit 
surface  is  smooth  enough  to  be  controlled  experimentally,  but  rough  enough 
to  be  out  of  the  realm  of  purely  ideal  jointing. 

After  grinding,  the  blocks  were  carefully  washed  in  water  to  remove  all 
grinding  debris.  The  blocks  of  rock  were  then  oven  dried  at  110°C,  removed 
from  the  oven,  and  allowed  to  cool  and  air-equilibrate  with  moisture  at  mom 
temperature.  For  subsequent  test  runs  on  blocks  already  run  once  or  more, 
the  blocks  were  again  carefully  washed,  oven  dried,  and  air-moisture  equi- 
librated. 

Conduct  of  the  Tests 

Test  series  were  performed  on  limestone  with  normal  stresses  of  0. 75, 
1.52,  1.95,  3.0,  5.0,  10. 0,  and  20. 0 kilograms  per  square  centimeter. 
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On  the  basalt,  tests  were  run  with  normal  stresses  of  0.  75,  1. 50,  3. 0,  5. 0, 

10.0,  and  20.0  kilograms  per  square  centimeter,  and  cm  the  quartzite,  1.5, 

3.0,  5. 0,  and  10. 0 kilograms  per  square  centimeter.  Two  series  were  run 
on  limestone  at  0. 75  and  3. 0 kg/cm2,  and  on  basalt  at  10. 0 kg/cms.  Gener- 
ally, each  test  series  is  a group  of  tests  performed  on  one  set  of  rock  blocks. 

This  series  is  given  a number,  prefixed  I,  for  limestone,  Q for  quartzite, 
or  B for  basalt.  Each  run  across  a block  uninterrupted  by  unload  of  the 
normal  stress  is  regarded  as  a test,  and  is  identified  by  a letter  suffixeu  to  the 
series  number.  Each  day's  run  on  a bltxik  is  identified  by  a number  added 
behind  the  test  letter.  Thus  the  third  test  series  on  basalt,  first  test,  second 
day's  run  is  Test  B3A2. 

The  series  attempted  to  establish  temperature  ranges  for  the  initiation  and 
cessation  of  stick-slip,  and  to  determine  the  effect  of  temperature  change  on 
the  stick,  slip,  and  shear  load  relaxation.  Additionally,  an  attempt  was  made 

s 

to  ascertain  the  effect  of  the  introduction  of  preground  rock  debris  between  the 
two  test  blocks. 

In  order  to  examine  the  stick-slip  properties  of  the  rock  with  respect  to 
temperature,  movement  of  the  blocks  was  begun  at  30°C.  After  observation 
of  the  mode  of  movement  at  the  lower  temperature,  the  temperature  was  raised 
incrementally.  For  example  stick-slip  had  initiated  at  30°C,  the  temperature 
was  raised  to  125°,  and  subsequently  to  200°  to  study  the  effect  of  increased  temp- 
erature. 

Observations  were  made  of  time  between  slips  (i.e. , time  of  stick),  amount 
of  slip,  shear  load  relaxation  during  slip,  and  total  shear  load  required  to  ini- 
tiate each  slip.  From  these  measured  parameters  the  coefficient  of  static  fric- 
tion at  the  instant  prior  to  slip,  and  the  cumulative  movement  were  calculated. 

The  character  of  the  slip  in  terms  of  time  length  of  slip,  abruptness,  shape  of  the 
slip  curve,  and  existence  of  creep  between  slips  was  also  observed.  The  tempera- 
ture was  then  reduced,  and  observations  were  made  of  movement  quantities  at 
various  temperatures  during  unloading.  The  temperature  of  cessation  of  stick- 
slip  was  noted. 
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In  some  light  load  tests  on  limestone,  the  most  accurate  means  of  quan- 
tifying the  stick-slip  results  was  time  between  slips.  While  an  absolute  one- 
to-one  correlation  of  length  of  slip  and  time  between  slips  was  not  observed, 
for  the  limestone  tests  the  correlation  coefficient  between  the  means  is  0.85 
to  0.98,  and  means  were  used  in  all  quantifications  in  order  to  observe  patterns. 
In  limestone  tests  with  normal  stress  of  over  3. 0 kg/cm2 , and  in  all  tests  on 
basalt  and  quartzite,  test  results  will  usually  be  presented  using  values  of  the 
length  of  slip.  This  parameter  is  preferred  because  there  are  fewer  gaps  in 
the  record,  and  because  it  presents  a better  picture  of  the  pattern  in  basalt. 
Shear  load  relaxation  data  parallels  that  for  slip,  but  is  not  as  complete  in 
the  limestone  tests. 
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DESCRIPTION  OF  ROCKS  TESTED 
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Three  types  of  rocks  were  used  in  the  testing  program,  all  being  from 
commercial  quarry  sources. 

Salem  Limestone 

The  limestone  tested  was  the  familiar  building  limestone  from  the  vicinity 
of  Bedford,  Indiana.  This  rock  is  from  the  Salem  Limestone  formation  of 
Mississippian  age.  The  Salem  Limestone  is  a very  light  gray,  slightly  porous, 
bioclastfc  limestone,  the  clastic  particles  being  composed  of  rounded  fossils 
shells  of  foraminifera,  with  measured  sizes  up  to  0.82  mm,  occasional  frag- 
ments of  bryozoans,  and  possibly  fragments  of  other  organisms.  The  fossil 
fragments  lie  in  a cement  of  optically  continuous  calcite.  Quarry  information 
stated  that  the  formation  contains  up  to  two  percent  quartz,  but  the  specimens 
examined  petrographic  ally  were  100%  calcite.  The  measured  mean  fragment 
diameter  is  0.  Gl  mm  with  a standard  deviation  of  0.20  mm.  Unconfined  com- 
pressive strength  was  evaluated  as  400  kilograms  per  square  centimeter. 


Basalt 

This  material  is  a basalt  of  Tertiary  age  from  Sonoma  County,  California, 
supplied  by  Smith’s  Granite  and  Marble  Works,  Napa,  California.  It  is  a 
dark  gray,  massive  basalt  with  interlocking  crystalline  texture.  Altered  zones 
which  probably  represent  flow  boundaries  are  spaced  about  one  centimeter  apart. 
Glass,  corroded  olivine  crystals,  microcrystalline  augite,  and  occasional  garnet 
and  magnetite  crystals  are  contained  in  a matrix  of  plagioclase  feldspar  (labra- 
corite)  laths.  Mean  length  of  plagioclase  laths  is  0.13  mm,  standard  deviation 
0.03  mm.  Some  twinned  plagioclase  crystals  are  as  large  as  0.6  mm.  Larger 
rectangular  crystals  of  epidote,  up  to  0.6G  x 0.23  mm,  are  scattered  throughout 
the  rock.  The  feldspar  laths  are  arranged  subparallel  to  the  flow  direction  in- 
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dicated  by  the  flow  boundaries.  Labradorite  and  augite  were  identified  by  X-ray 
diffraction;  the  other  minerals  did  not  show  up  in  the  X-ray.  Mineral  composi- 
tion as  determined  by  the  integrating  stage  is:  plagioclase  87%,  glass  4%,  epi- 
dote  1.5%,  garnet  0.5%,  olivine  3%,  augite  4%,  ....^gnetite,  trace. 

The  average  unconfined  compressive  strength  perpendicular  to  the  flow 
boundaries  was  1369  kg/cm2,  with  a standard  deviation  of  138  kg/cm2.  Two 
specimens  were  tested  with  the  major  principal  stress  axis  parallel  to  the  slow 
boundaries.  One  broke  along  the  boundary  with  an  unconfined  compressive 
strength  of  1181  kg/cm2,  the  other  broke  within  the  unaltered  rock  with  a 
strength  of  1785  kg/cm2.  The  altered  zones  are  not  believed  to  be  extensive 
enough  to  affect  the  surface  properties  of  the  rock.  Tests  were  run  parallel 
and  perpendicular  to  the  flow  boundary  zones,  and  also  on  blocks  with  very 
weakly  developed  altered  zones. 

Sioux  Quartzite 

The  Sioux  Qiartzite  is  a dense,  fine  grained,  very  hard  pre-Cambrian 
quartzite  which  outcrops  in  extreme  northwestern  Iowa  and  southwestern 
Minnesota,  where  it  is  quarried  for  aggregate  and  for  monument  stone.  The 
pink  to  red  color  is  attributed  to  iron  oxide  coatings  on  the  clastic  grains,  which 
in  turn  are  overgrown  by  authigenic  quartz  in  crystallographic  continuity  with 
the  adjacent  sand  grains.  The  result  is  a characteristic  quartzite  interlocking 
structure  with  minimum  wids.  The  mineralogical  content  is  primarily  quartz 
with  accessory  zircon,  anatase,  rutile,  hematite,  and  orthoclase  (Beyer,  1897). 
The  rock  was  not  examined  petrographically  or  tested  for  strength. 
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TKST  RESULTS 


Results  will  be  discussed  in  several  categories:  temperature  for  initiation 
of  stick-slip  (from  smooth  slip),  character  of  stick-slip,  changes  in  stick  time 
or  slip  distance  with  time,  changes  in  surface  character,  temperature  as  a con- 
trolling factor,  and  changes  in  the  coefficient  of  sliding  friction.  Detailed  test 
results  are  given  in  Drennon  (1972)  and  Drennon  and  Handy  (1972). 

Initiation  of  stick-slip 

Initiation  of  stick-slip  was  found  to  depend  upon  four  factors:  (1)  material, 

(2)  previous  frictional  history  of  the  specimen,  (3)  temperature,  and  (4)  normal 
stress.  The  most  extensive  series  comprising  sixty  tests  was  performed  on  the 
limestone,  in  which  temperature  and  normal  stress  were  found  to  have  a pronounced 

2 

effect  on  the  initiation  of  stick-slip.  In  thirteen  tests  with  a normal  load  of  0.75  kg/cm  , 
no  stick-slip  occurred  below  a temperature  of  100°C  (Table  1).  Under  higher  applied 
normal  stresses  stick-slip  became  more  likely  below  100°C,  but  still  occurred  at  this 
low  temperature  in  only  about  one-fourth  of  the  tests.  On  the  other  hand,  stick-slip 
always  was  initiated  when  the  temperature  was  increased,  usually  (58%  of  the  tests) 
beginning  in  the  range  100-179°C.  Extended  movement  or  debris  artificially  added 
between  the  limestone  blades  usually  inci  eased  the  likelihood  of  smooth  slip  ex- 
tending to  higher  temperatures. 


Table  1.  Temperature  of  initiation  of  stick-slip  in  limestone 


Normal  stress 

Temperature 

range 

Number 
of  tests 

2 

0. 75  kg/cm 

30  - 99°C 

0 

100  - 179 

9 

180  - 205 

4 

over  205 

0 

2 

1.52  kg/cm 

30  - 99°C 

2 

100  - 179 

5 

180  - 205 

0 

over  205 

0 

. 2 
1.95  kg/cm 

30  - 99°C 

1 

100  - 179 

0 

180  - 205 

6 

over  205 

2 

2 

3.00  kg/cm 

30  - 99°C 

5 

100  - 179 

7 

180  - 205 

0 

over  205 

0 

2 

5. 00  kg/cm 

30  - 99°C 

1 

100  - 179 

6 

180  - 205 

0 

over  205 

0 

2 

10. 00  kg/cm 

30  - 99°C 

1 

100  - 179 

5 

180  - 205 

2 

over  205 

0 

2 

20. 00  kg/cm 

30  - 99°C 

1 

100  - 179 

3 

180  - 205 

C 

over  205 

0 

33 


Initiation  of  stick-slip  in  basalt  is  more  irregular  in  pattern  as  well  as  more 
frequent  in  occurrence  than  in  limestone.  Most  tests  on  basalt  produced  stick- 
slip  immediately  upon  initiation  of  movement  regardless  of  normal  stress  or 
previous  frictional  history.  These  factors  do  seem  to  have  some  influence  on 
the  initiation  of  stick-slip  in  basalt,  but  certainly  less  than  in  limestone. 

Smooth  slip  occurred  at  the  initiation  of  movement  mostly  at  higher  normal 

2 

stresses.  Three  tests  with  a normal  stress  of  10  kg/cnr  began  smooth  slip 

o 2 

on  washed  blocks  at  30  , whereas  the  initial  basalt  test  at  0. 75  kg/cm  went  into 

stick-slip  after  a short  period  of  irregular  smooth  slip  on  initiation  of  movement. 

These  two  test  runs  continued  and  ended  in  smooth  slip. 

Basalt  tests  with  debris  introduced  between  the  blocks  tended  to  begin  sliding 
with  a chaotic  sticking,  slipping,  and  sliding  in  an  extremely  irregula.-  manner. 
Sometimes  long  periods  of  smooth  slip  would  alternate  with  stick-slips  which 
were  more  likely  to  be  all  pervasive  in  basalt  in  the  middle  range  of  the  normal 
stresses  tested.  However,  with  the  exception  of  Test  Series  B5,  the  high  normal 
stress  groups  which  b*  gan  in  smooth  slip  moved  into  stick-slip  withoiK  any  al- 
teration of  temperature  or  other  test  condition.  Stick-slip  seems  to  be  the  normal 
mode  of  frictional  sliding  for  the  basalt  in  the  range  of  normal  stressi.  and  with 
the  testing  machine  used. 

Quartzite  began  movement  in  the  stick-slip  mode  immediately  upon  'nitiation 
of  movement  at  30°C,  and  the  smooth  slip  mode  of  movement  occurred  n the 
three  more  lightly  loaded  tests  after  the  blocks  had  remained  overnight  ( approxi- 
mately fifteen  hours)  at  a temperature  of  200°.  In  those  three  cases,  movement 
began  with  smooth  slip  and  retained  that  mode  of  movement  with  ai.  increase  of  tem- 
perature. A reduction  in  temperature,  on  the  other  hand,  resulted  in  resumption 
of  stick-slip  at  temperatures  close  to  100°C,  and  stick-slip  continued  upon  re- 
elevation of  temperature  to  200°.  This  series  of  occurrences  was  the  principal 
difference  between  results  from  quartzite  and  the  other  rocks  tested,  though 
it  resembles  some  of  the  basalt  tests  wherein  stick-slip  reduction  occurred  with 
length  of  travel. 
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No  debris  tests  were  conducted  on  quartzite,  since  practically  no  debris  was 
generated  by  sliding. 

Character  of  stick-slip 

At  temperatures  above  100°  most  limestone  stick-slips  were  abrupt,  with 
measurable  time  lengths  from  the  beginning  to  the  end  of  slip  movement  being 
relatively  rare  at  normal  oscillograph  rates  of  20  or  50  seconds  per  centimeter. 
Some  random  measurements  at  125°  ranged  from  0.60  to  8.0  seconds,  but  long 
time  lengths  were  most  unusual.  At  30°,  limestone  slips  sometimes  had  a time 
length  of  up  to  20  seconds,  but  such  lengths  were  not  common  during  regular 
stick-slip.  They  were  normal  when  transition  to  smooth  slip  was  imminent. 

In  limestone,  slips  tend  to  break  off  abruptly  at  the  peak  of  load  on  the  load  - 
time  record  (Figure  3),  and  to  gather  time  length  by  slowing  down  as  the  slip 
comes  to  a stop.  The  time  length  of  limestone  slips  stretches  out  as  the  point 
of  smooth  slip  approaches  during  temperature  decrease. 

Stick -slips  in  the  basalt  tests  presented  a different  appearance.  Figure  3 
represents  a common  phenomena,  called  "30°-basalt  stick-slip"  because  it  was 
most  often,  though  not  exclusively,  observed  at  that  temperature.  This  is  the 
tendency  to  alternate  relatively  large  slips  and  relatively  large  shear  load  relaxa- 
tions with  small  slips  and  small  shear  load  relaxations.  Creep  frequently  appears 
between  slips.  This  phenomenon  would  often,  but  not  invariably,  cease  when  the 
temperature  was  raised. 

The  time  length  of  slip  of  basalt  slips  tended  to  be  longer  than  those  of  lime- 
stone. 

All  tests  of  quartzite  were  very  abrupt  (Figure  3),  with  a small  amount  of 
creep  between  slips.  The  time  length  of  slip,  where  measured  during  Tests  Ql 
and  Q3,  was  less  than  0.05  second.  This  compares  with  normal  limestone  times 
of  0.6  to  8.0  seconds,  and  basalt  time  lengths  which  ranged  from  1 to  22  seconds 
during  a single  test.  As  a rule,  the  magnitude  of  slip  of  the  quartzite  was  larger 
than  that  of  the  basalt  or  of  the  limestone. 
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A comparison  of  the  load  - time  plots  of  Figure  3 shows  clearly  the 
abrupt  break  of  the  limestone  and  quartzite  and  the  "nose  over"  of  load  of 
the  basalt.  Though  this  figure  is  one  of  the  more  exaggerated  examples  of 
"nose-over"  of  load  and  small  creep  movement  before  the  major  slip,  a similar 
effect  is  characteristic  of  all  basalt  tests.  The  shape  of  the  stick-slip  curve 
for  basalt  is  similar  to  that  for  silver,  a very  ductile  material,  described  by 
Brockley  and  Davis  (1968). 

Stick-slip  during  long  runs  at  constant  temperature 

During  most  of  the  tests  the  upper  block  was  allowed  to  move  for  considera- 
ble periods  of  time  over  the  lower  block  at  a constant  temperature.  Measure- 
ments were  made  of  each  time  between  slips,  slip,  and  shear  load  relaxation 
upon  slip.  Groups  of  ten  successive  slips  were  then  averaged  and  studied  to 
determine  whether  any  trends  occurred  with  time,  or  with  accumulated  total 
slip. 

The  most  extensive  constant  temperature  test  at  the  lowest  normal  stress  of 
0.75  kg/cm2  was  the  190°  run  during  limestone  Test  LIB,  which  showed  a steady 

increase  in  slip  distance  throughout  the  run.  This  pattern  of  increase  of  slip  with  time 

2 

was  observed  in  tests  at  9. 75  kg/cm  at  all  different  vemperatures.  It  was 

also  the  pattern  of  change  of  slip  with  "'ccumulated  movement  at  constant  tempera- 
2 

ture  at  1.52  kg/cm  , but  this  tendency  for  time  between  slips  and  slip  of  lime- 
stone to  increase  with  distance  moved  is  a phenomenon  restricted  to  small  nor- 
mal stresses. 

2 

In  tests  with  normal  stresses  from  1. 95  kg/cm  and  up,  the  situation  is 
sharply  reversed.  In  tests  beginning  with  clean  blocks,  the  tendency  is  consis- 
tently for  a reduction  in  time  between  slips,  or  slip,  as  movement  progresses, 
approaching  a constant  value  of  about  0. 0040  cm.  This  value  would  be  a func- 
tion of  the  materials,  temperature,  and  test  apparatus  as  well  as  of  normal 
stress.  In  tests  which  were  stopped  at  125°  and  then  re-started  at  the  same 
temperature,  the  slip  value  remains  the  same  as  had  been  attained  before 
stopping  the  movement.  The  scatter  of  the  first  group  upon  resumption  of 
movement  is  usually  rather  large. 
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The  response  of  the  basalt  blocks  to  constant -rate-  of-movemer.t  direct 
shear  tests  at  a single  temperature  is  in  such  great  contrast  to  that  of  the 
limestone  as  to  suggest  differing  mechanisms,  or  at  least  different  responses 
to  the  same  mechanism.  As  with  the  limestone,  response  of  the  basalt  de- 
pends upon  normal  stress  and  temperature,  but  previous  frictional  history 
has  less  effect. 

2 

At  the  light  normal  stress  of  0. 75  kg/cm  , little  change  in  amount  of 

slip  occurred  with  continued  running  at  a constant  temperature.  With  an  increase 

2 

in  normal  stress  1.50  kg/cm  , conditions  begin  to  change  slightly,  and  be- 
2 

ginning  at  3. 0 kg/cm  , the  mode  of  sliding  most  typical  of  basalt  tests  at 
constant  temperature  takes  over.  Within  the  same  test  run,  a steep  and  linear 
(Figure  4)  or  semilogarithmic  (Figure  5)  increase  in  amount  of  stick,  slip, 
and  shear  load  relaxation  occurs. 

In  basalt  tests  at  lower  stresses,  no  stable  value  of  slip  was  reached.  At 
higher  normal  stress  a tendency  was  observed  to  seek  a stable  value. 

Continuous  runs  of  quartzite  at  a single  temperature  tended  towards  constant 
values  of  slip,  either  increasing  to  the  value,  or  decreasing  to  the  value. 

Surface  changes 

All  blocks  were  prepared  for  sliding  in  a like  manner,  being  ground  with 

a number  45  grit  aluminum  oxide  grinding  wheel.  The  surface  of  the  blocks 

was  examined  under  a binocular  microscope.  One  limestone  block  and  one  basalt 

block  were  further  examined  under  a scanning  electron  microscope. 

Under  the  binocular  microscope  the  ground  limestone,  being  bioclastic,  gives 

a rough  appearance.  Powder-sized  particles  are  first  visible  to  the  naked  eye 

2 

after  sliding  at  a normal  stress  of  1. 95  kg/cm  . Examination  under  the  bino- 
cular microscope  showed  some  calcite  chips  with  sharp  edges  apparently  re- 

2 

suiting  from  brittle  fracture  at  normal  stress  as  low  as  1.52  kg/cm  . At 
2 

1,95  kg/cm  normal  stress  all  visible  debris  could  be  washed  off,  and  no  evi- 
dence of  damage  remained. 

2 

Beginning  with  a normal  stress  3. 0 kg/cm  visible  chipping  occurred. 
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Two  of  the  broken  chips  were  examined  under  the  scanning  electron  micro- 
scope. The  surface  of  one  chip  was  washed  and  dried  according  to  the  normal 
procedure  before  examination,  and  the  other  fragment  was  handled  carefully  to 
avoid  disturbing  the  debris  present  on  the  face.  The  washed  chip,  examined  at 
300x  showed  little  evidence  of  sliding.  Some  cleavage  faces  seemed  pulled  up 
from  high  points.  These  were  usually  on  the  rear  side  looking  in  the  direction 
of  travel.  Little  evidence  was  found  in  any  picture  for  conical  or  otherwise 
sharp-peaked  asperities,  but,  as  would  be  expected  from  a limestone,  the  high 
poi  nts  were  relatively  flat  cleavage  faces.  These  faces  built  up  to  flat  contact 
areas  which  were  a small  percentage  of  the  total  surface  of  the  block.  The 
fragment  examined  with  debris  in  place  showed  many  fractured  crystal  frag- 
ments. Some  debris  is  present  on  the  higher,  flat  areas,  but  most  has  fallen 
into  the  lower  void  areas  which  are  quite  extensive  in  this  limestone.  Figure  6 
and  other  pictures  show  that  the  debris  is  not  created  by  shearing  off  the  tops 
of  high  points  which  seem  little  affected  by  sliding,  but  by  plucking  of  small 
cleavage  fragments  from  the  edges  of  the  higher  points  and  carrying  these 
f ragments  both  into  voids  and  onto  the  surface  of  the  "asperities”. 

The  basalt,  being  a much  harder  material,  showed  little  damage  due  to 

2 

sliding  in  tests  with  normal  stresses  of  less  than  3. 0 kg/cm  . In  the  3.  0 and 

5. 0 normal  stress  tests  minor  striations  appear  on  both  the  top  and  bottom 

2 

blocks.  In  a test  at  20. 0 kg/cm  , a heavy  concentration  of  rock  dust  appears 
in  parallel  lineations.  The  striae,  examined  under  the  binocular  microscope, 
appear  much  more  evenly  distributed  and  less  strongly  evident  than  with  the 
naked  eye.  The  fragments  are  visibly  due  to  the  brittle  fracture  of  small 
crystals  of  feldspar,  as  they  tend  to  be  white  in  color,  as  against  the  black 
or  gray  augite. 

Basalt  fragments  were  also  examined  under  the  scanning  electron  micro- 
scope. In  the  clean  basalt,  ground  with  numben-45  grit  but  never  subjected  to 
sliding,  asperities  appear  as  large,  flat  ground-off  crystals  of  feldspar, 

(Figure  7a).  The  planeness  of  some  of  the  flat  crystal  faces  and  their  rela- 
tively large  extent  are  remarkable.  The  fragment  broken  from  an  upper  block 
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after  sliding  at  10  kg/cm“  was  also  examined  under  the  scanning  electron 
microscope.  Compared  to  the  limestone,  remarkably  little  debris  was 
noted  (Figure  7b).  Most  of  the  debris  had  fallen  into  lower  areas,  but  some 
is  present  on  higher  areas.  As  with  the  limestone,  most  of  the  debris  appears 
to  have  been  removed  by  brittle  fracture  from  the  edges  of  the  higher,  flat 
crystal  faces. 

Because  of  the  hardness  of  the  quartzite,  very  little  surface  damage  was 

2 

noted.  Two  striations  were  observed  after  a test  at  3.0  kg/cm  , but  the 
amount  of  damage  was  not  sufficient  to  be  noticeable  after  the  dust  along  the 
marks  was  washed  away. 

Effect  of  added  debris 

Sliding  obviously  produced  debris,  and  study  of  the  stick-slip  patterns 
in  the  limestone  showed  obvious  changes  in  the  character  of  slip  and  time 
between  slip  patterns  with  continued  forward  movement.  It  was  thought  that 
the  debris  produced  upon  sliding  might  be  very  influential  in  determining  the 
character  of  stick-slip.  Hoskins  et.  al.  (1968)  concluded  that  debris  produced 
by  the  sliding  would  cause  stick-slip  to  disappear,  a conclusion  also  reached 
by  Dieterich  (1970). 

It  was  decided,  therefore, to  artificially  introduce  debris  in  the  form  of 
ground  rock  fragments  between  the  two  sliding  blocks  from  the  beginning  of 
sliding,  in  order  to  observe  the  nature  of  stick-slip  sliding  when  energy  was 

not  being  expended  in  producing  the  debris. 

2 

In  a test  at  1.95  kg/cm“,  0.02  grams  of  powdered  debris  was  added.  Though 
this  was  more  debris  than  was  collected  from  the  surface  after  the  preceding 
test  no  change  in  stick -slip  pattern  was  noted. 

In  subsequent  debris  tests  on  limestone  at  higher  normal  stresses,  0.05 
grams  of  debris  were  added.  The  added  debris  was  not  all  sub-microscopic  in 
size,  as  was  that  created  by  sliding,  but  ranged  in  diameter  from  0. 05  mm 
downwards.  This  debris  had  a significant  effect  on  the  nature  of  stick-slip, 
especially  at  the  initiation  of  movement.  In  the  tests  with  added  debris,  the 
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initial  group  of  slips  is  the  smallest  in  length  of  the  entire  test,  and  the  amount 
of  slip  then  increases  rapidly  to  a value  similar  to  that  finally  attained  in  the 
tests  begun  on  clean  blocks.  Thus,  addition  of  a finite  amount  of  gouge  between 
the  sliding  blocks  appears  to  speed  up  what  occurs  on  an  initially  clean  joint  in 
a longer  time. 

The  result  of  adding  debris  between  blocks  of  basalt  is  very  similar  to  the 

result  on  limestone.  In  the  lightest  normal  stress  test  with  debris  added 
2 

(3. 0 kg/cm  ),  a chaotic  alternation  of  stick-slip,  smooth  slip,  and  creep  oc- 
curred at  30°.  Upon  raising  the  temperature  to  125°,  stick-slip  became  esta- 
blished at  a constant  mean  value,  but  with  a standard  deviation  which  was  large 

2 

and  which  increased  with  further  movement.  In  a test  B4  at  5. 0 kg/cm  , the 
debris  test  also  seemed  to  reach  a constant  mean  value  of  slip  rapidly,  but  this 
value  was  accompanied  by  a large  standard  deviation. 

Debris  tests  were  not  performed  on  quartzite,  as  already  indicated. 

The  creation  of  addition  of  quantities  of  debris  between  sliding  blocks  of 
a joint  system  thus  has  a pronounced  effect  on  stick-slip  between  the  blocks, 
namely  a more  rapid  approach  to  equilibrium  conditions:  In  the  amounts 
added,  debris  did  not  cause  a cessation  of  stick-slip. 

Effect  of  temperature 

Tbmperature  proved  to  be  a major  factor  in  controlling  slip,  time  between 
slips,  and  shear  load  relaxation.  In  all  cases,  increasing  the  temperature  of 
the  rock  while  it  was  sliding  in  a stick-slip  mode  caused  an  increase  in  the  time 
between  slips,  an  increase  in  the  length  of  individual  slips,  and  an  increase  in 
the  amount  of  relaxation  of  shearing  load  upon  slip  (Figures  8 and  9). 

Temperature  could  be  raised  rapidly  by  the  heaters  in  the  temperature  control 
enclosure  of  the  shear  box,  but  had  to  be  reduced  more  slowly.  When  temperature 
decreased,  the  amount  of  stick,  slip,  and  shear  load  relaxation  all  decrease 
with  temperature  until,  in  the  limestone  and  sometimes  in  basalt,  smooth  slip 
finally  became  the  mode  of  frictional  movement. 

The  temperature  of  cessation  of  stick-slip  (Table  2)  varied  from  test  to  test 
for  limestone  and  was  lower  than  the  temperature  of  initiation  of  stick-slip  in  all 
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Figure  8.  Change  of  average  slip  with  temperature;  Basalt 
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cases  where  some  temperature  rise  from  a lower  temperature  was  required 
to  induce  the  beginning  of  stick-slip. 


Table  2.  Temperature  of  cessation  of  stick-slip  in  limestone 


Temperature  of  cessation 
Degrees  Centigrade 

Number  of  Events 
Limestone  Basalt 

30-40 

0 

3 

40-50 

7 

0 

51-00 

3 

1 

01-70 

1 

0 

71-80 

4 

0 

81-90 

2 

0 

91-100 

9 

0 

101-200 

2 

2 

200- 

4 

Would  not  stop 

0 

0 

Generally  the  temperature  of  cessation  varied  with  normal  stress,  with 

most  of  the  lower  temperatures  representing  tests  with  higher  normal  stress. 

In  the  basalt,  elimination  of  stick-slip  by  a reduction  in  temperature  was 

more  irregular  than  in  the  limestone,  anti  in  a number  of  tests,  stick-slip 

could  not  be  eliminated  at  all  by  temperature  reduction.  The  few  tests  in 

which  smooth  slip  was  induced  by  temperatures  reduction  were  in  the  higher 

normal  stress  tests,  with  smooth  slip  initiating  at  low  temperatures  (Table  2). 

A peculiarity  of  the  basalt  was  the  frequency  of  the  mode  of  sliding  moving 

into  smooth  slip  during  continuous  running  at  temperatures  above  30°.  This 

2 

occurred  twice  at  0.  75  kg/cm“  normal  stress,  and  also  occurred  twice  at 
2 2 

20. 0 kg/cm“.  At  20. 0 kg/cm“  the  block  had  moved  the  fairly  extensive  dis- 
tance of  0.2077  cm  when  smooth  sliding  occurred,  but  the  same  block  in  the 
same  test  (BOA)  subsequently  continued  stick-slip  down  to  37°  without  going 
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into  a smooth  slip  mode  of  sliding. 


In  ihe  quartzite  tests,  smooth  slip  occurred  only  after  a period  of  rest  at 
200°C,  and  lowering  the  temperature  below  100°C  caused  a return  to  stick-slip. 
This  is  just  the  opposite  of  the  reaction  of  the  limestone  and  the  basalt,  probably 
reflective  of  a different  surface  property. 

The  coefficient  of  static  friction  was  evaluated  from  the  maximum  shearing 
load  immediately  before  the  beginning  of  slip,  divided  by  the  normal  load.  No 
evaluation  was  made  of  the  dynamic  coefficient  of  friction,  which  in  stick-slip 
studies  is  usually  taken  as  the  load  at  one  half  of  slip  distance  divided  by  normsil 
load.  The  behavior  of  the  coefficient  of  static  friction  was  seen  to  vary  with  vock 
type,  frictional  history,  normal  stress,  and  temperature.  The  evaluation  of 
coefficient  of  friction  was  carried  through  temperature  changes  after  comparison 
of  measured  values  of  coefficient  of  friction  upon  unloading  with  the  calculated 
coefficient  of  friction  of  the  last  slip  immediately  before  unloading  showed  the 
two  values  to  be  the  same.  Gaps  in  the  solid  lines  on  the  plots  of  coefficient  of 
friction  versus  cumulative  distance  moved  represent  the  end  of  a day's  run  and, 
usually,  unload  of  the  shear  load  before  start  of  the  next  day’s  run. 

Under  light  normal  loads  the  coefficient  of  friction  of  limestone  tended  to 
decrease  to  a constant  value  (Figure  10).  This  was  the  exception,  howex'er, 
because  under  higher  normal  loads  (Figure  11)  and  in  all  tests  on  basalt  and  on 
quartzite,  as  the  test  progressed  friction  increased  to  approach  a constant  value. 
This  is  illustrated  in  Figures  12  and  13. 

The  final  coefficient  of  friction  for  limestone  tended  to  increase  with  higher 
normal  stress  (Figure  14)  but  not  for  the  basalt  (Figure  15)  perhaps  because 
the  sliding  distance  was  insufficient  for  friction  to  reach  equilibrium  value.  The 
data  for  quartzite  show  a similar  trend  to  limestone  but  are  rather  sparse. 

The  effect  of  temperature  on  friction  is  shown  in  Figures  11  and  13,  an  in- 
crease in  temperature  increasing  the  coefficient  of  sliding  friction. 

Debris  tests  showed  similar  characteristics  of  a very  rapid  rise  from  a low 
initial  coefficient  of  friction  to  a stable  value.  This  often  is  approximately  main- 
tained despite  any  alteration  of  temperature  up  or  down  or  any  alteration  of  the 
mode  of  sliding  to  smooth  slip  from  stick-slip. 
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COEFFICIENT  OF  STAT 


The  values  of  debris-aaded  coefficient  of  friction  for  the  limestone  con- 
trasts sharply  with  the  values  for  blocks  which  were  slid  with  initially  clean 

2 

surfaces,  being  much  higher  and  resembling  that  attained  in  the  20.0  kg/cm 
normal  stress  test  (Figure  14). 

In  the  basalt  tests,  addition  of  debris  docs  not  change  the  pattern  of  in- 
crease of  coefficient  of  friction  with  accumulated  slip,  but  it  does  seem  to  speed 
is  up  considerably,  being  reached  after  0.1-0. 3 cm  sliding.  In  all  cases,  the 
coefficient  of  static  friction  of  debris  tests  in  basalt  is  lo  ver  than  the  coef- 
ficient of  tests  begun  with  clean  blocks,  debris  tests  approaching  a value  of 
0.6  while  clean  tests  reach  initial  stability  at  0.9  to  1.0.  Temperature  changes 
cause  a greater  reaction  in  clean  tests  than  in  debris  tests. 


CUMULATIVE  SLIP,  cm 


Figure  10.  Pattern  of  coefficient  of  static 

friction  for  limestone.  Test  L7A 

a =3.0  kg/cm3, 
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Figure  15.  Basalt  and  quartzite  final  coefficients  of 
friction , 


DISCUSSION 


A great  deal  of  information  on  the  stick-slip  patterns  of  rock,  some  of  it 
seemingly  trending  in  opposing  directions,  has  been  gathered.  It  is  believed , 
however,  that  most  of  the  patterns  of  change  in  type  and  amount  of  slip  and 
coefficient  of  friction  with  temperature,  cumulative  slip,  and  time  can  be  in- 
tegrated into  a single  model  for  the  mechanism  of  stick-slip  in  rock,  which, 
while  perhaps  not  correct,  will  at  least  aid  in  organizing  results. 

It  is  obvious  from  these  tests  that  stick-slip  will  occur  in  any  rock,  given 
a soft  loading  system  such  as  gravity,  that  used  in  the  experiments,  or  almost 
any  other  system  with  physical  discontinuities  allowing  a degree  of  freedom 
perpendicular  to  the  surface  of  sliding  (Tolstoi,  1967).  Stick-slip  was  the 
prevalent  mode  of  sliding  in  a limestone,  a basalt,  and  a quartzite.  Only  the 
latter  contains  a significant  amount  of  quartz.  Stick-slip  can,  then,  be 
regarded  as  a universal  phenomenon  in  rocks,  given  the  proper  conditions 
of  movement  and  force. 

The  results  of  the  tests  on  the  limestone  suggest  that  smooth  slip,  when- 
ever it  can  be  induced,  may  be  caused  by  a film  of  adsorbed  water.  On  the 
majority  of  low-normal-stress  tests,  stick-slip  did  not  commence  until  the 
temperature  was  raised  above  100°.  Kven  on  the  higher  normal  stress  tests, 
reduction  of  temperature  after  a period  of  stick-slip  at  a higher  temperature 
to  a level  permitting  the  adsorption  of  water  resulted  in  the  cessation  of  stick- 
slip.  Temperature  usually  had  to  be  raised  above  100°  for  stick-slip  to  re- 
initiate. A similar  relationship  between  friction  and  adsorbed  water  was  found 

by  Peterson  and  Murray  (1967)  for  ceramic  materials. 

2 

The  tests  with  normal  stress  of  3 kg/cm  and  oven  and  the  lighter  loaded 
blocks  which  began  stick-slip  below  100°  were  clean  blocks.  The  absence 
of  debris  would  permit  some  major  asperities  to  come  into  direct  contact, 
penetrating,  as  it  were,  the  intervening  film  of  adsorbed  water  and  thus  ini- 
tiating stick-slip  at  low  temperature.  At  the  higher  normal  stresses  the 
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adsorbed  water  may  be  reduced  to  a monolayer  or  removed  entirely  on  the 
higher  flat  asperities.  The  presence  of  this  water  may  also  cause  anomalous 
creep  (Westbrook  and  Jorganson,  1905)  of  the  crystals  in  contact,  increasing 
the  surface  area  and  the  bonding.  This  is  known  as  the  Rebinder  effect 
(Westwood  et.  al. , 1907).  Once  initiated,  stick-slip  might  continue  as  the 
asperities,  and  then  the  debris,  penetrated  and  disturbed  the  film  of  adsorbed 
water. 

Further  support  to  the  concept  of  smooth  sliding  at  low  temperature  on  a 
boundary  layer  of  adsorbed  water  is  given  by  activation  energies.  In  one 
phase  of  these  experiments,  the  measurement  of  activation  energy’  required 
to  initiate  and  maintain  creep  of  limestone  blocks  separated  by  a joint  was 
attempted.  These  experiments  are  described  in  Appendix  A.  The  results 
were  found  to  have  a bearing  on  the  current  problem.  The  activation  energies 
obtained  were  found  to  be  temperature  dependent.  F 'll  creep  begun  below 
100°,  activation  energies  ranged  from  -5.12  to  -10.73  kilocalories  per  mole, 
averaging  -7.  54  ki  localories  per  mole.  This  value  compares  favorably  to 
the  value  for  activation  energy  of  water  of  -4  to  -5  kilocalories  per  mole 
determined  by  Glasstone,  Laidler,  and  Kyring  and  reported  by  Mitchell,  Singh, 
and  Campanella  (1959).  Further,  activation  energies  determined  at  higher 
temperatures  show  a definite  increase  with  temperature,  with  an  average  value 
of  -23. 0 kilocalories  per  mole  for  temperature  between  140°  and  200*\  This 
value  is  similar  to  that  determined  for  dry  sand,  about  -25  kilocalories  per  mole, 
by  Mitchell,  Singh,  and  Campanella  (1959).  According  to  Horn  and  Deere  (1952), 
the  oven  dried,  air-equilibrated  coefficients  of  friction  of  calcite  and  quartz 
are  almost  identical.  Thus  it  is  believed  that  most  adsorbed  water  is  being 
driven  off,  with  time,  when  temperatures  reach  above  100°,  and  that  limestone 
asperity  to  asperity  contact  is  taking  effect.  The  fact  that  the  average  value  of 
the  activation  energy  below  100°  is  above  that  of  liquid  water  is  believed  due  to 
a higher  energy  of  adsorption  (Noble  and  Demirel,  1959). 

In  most  cases,  the  basalt  did  not  begin  movement  with  smooth  slip,  or 
attain  smooth  slip,  even  at  low  temperature.  The  basalt  did,  however. 
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normally  display  "30°  basalt  stick-slip"  at  temperature  below  100°  upon 
initiation  of  movement.  This  irregular  motion,  often  accompanied  by  creep 
between  stick-slips,  is  believed  to  be  affected  by  adsorbed  water.  The 
bonding  between  the  crystals  of  basalt  is  believed  to  be  superior  to  that  be- 
tween those  of  the  limestone.  It  is  entirely  possible  that  the  surface  effect 
caused  by  adsorbed  water  causes  a time-dependent  reduction  of  surface 
hardness.  This  reduction  of  surface  hardness  occurs  in  many  nonmetallic 
crystals  such  as  MgO  (Westwood  et.  al. , 1967)  and  AlgCg  (Westbrook  and 
Jorganson,  1965).  If  it  occurs  more  readily  in  plagioclase  than  in  calcite 
it  would  explain  the  irregular  30°  basalt  movement,  as  caused  by  bonding  al- 
ternating with  water-water  contact. 

The  basalt  did  go  into  smooth  slip  upon  reduction  of  temperature  three 
times.  In  these  cases,  adsorbed  water  could  be  present  in  more  than  usual 
thickness,  and  smooth  slip  could  ensue.  In  all  basalt  temperature  reduc- 
tions, time  lengths  of  slip  stretch  out  as  temperature  reduces  and  bonding 
fails. 

The  rather  unique  behavior  of  the  quartzite  going  into  stick- slip  upon 
a reduction  of  temperature  below  100°C  suggests  a stronger  Rebinder 
effect,  also  indicated  by  the  well-known  property  of  water  as  an  anti -lubri- 
cant on  smooth  surfaces  of  quartz  (Horn  and  Deere,  1962). 

Above  100°  in  the  limestone  and  at  most  temperatures  in  the  basalt  and 
quartzite,  true  asperity  to  asperity  contact  probably  occurs.  These  contacts 
may  be  plastic  and/or  elastic.  Consideration  of  the  characteristics  and  pat- 
terns of  the  tests  leads  to  the  belief  that  plasticity  and  adhesion  by  bonding  of 
the  asperities  as  described  by  Bowden  and  Tabor  is  a critical  factor  in  the 
mechanism  of  stick-slip  in  rock. 

The  asperities  involved  have  been  shown  to  be  small  areas  of  flat  cleavage 
faces  which  are  slightly  higher  than  adjoining  areas  and  thus  can  come  into 
contact  with  the  opposing  flat  cleavage  face.  Above  100°  the  Rebinder  effect 
of  surface  softening  due  to  water  does  not  occur  but  the  increased  temperature 
may  encourage  plastic  deformation  of  the  asperities.  In  limestone,  confining 
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pressure  has  been  shown  to  have  a pronounced  effect  on  ductility.  At  low 
strains  and  confining  pressures  of  around  GOO  bars,  uniform  flow  occurs  at 
room  temperature  (Donath  ct.  al.,  1971).  If  the  asperities  occupy  one  per- 
cent of  the  apparent  area  of  contact,  which  is  not  an  unreasonable  value, 
pressures  of  up  to  2000  bars  were  applied  to  the  limestone  asperit  ies  in  the 
present  tests,  fully  sufficient  for  uniform  flow  and  plastic  deformation  to 
occur.  Upon  close  contact  of  the  asperities,  actual  chemical  bonding  is 
believed  to  take  place,  as  the  atoms  of  the  asperity  are  not  able  to  distin- 
guish between  the  two  sides  of  the  contact. 

When  bond  strength  is  exceeded  by  the  constantly  increasing  shear  force 
of  the  tests,  slip  occurs.  As  the  slip  rate  decreases,  bonding  can  again  take 
over  at  whatever  asperities  are  in  contact.  With  the  time  required  for  the 
shear  load  relaxation  to  recover  to  the  former  load  level  in  a soft  system 
where  the  load  cannot  follow  the  movement,  bonding  will  increase  with  time 
of  asperity  contact.  The  same  general  effect  is  believed  to  occur  in  the 
basalt  and  the  quartzite. 

Increase  of  temperature  allows  further  plastic  deformation  and  larger 
areas  of  contact  of  the  asperities  and  a longer  time  is  required  to  build  up 
the  required  shearing  force  to  break  the  developed  bonds.  The  strength 
of  bonding  thus  increases,  as  shown  by  Kosterin  and  Kraghelsky  (1962). 

The  excess  of  static  coefficient  of  friction  over  kinetic  coefficient  of  friction 
then  means  a larger  slip,  a larger  stick,  and  a larger  shear  load  relaxation 
with  increase  of  temperature. 

All  of  the  test  results  on  limestone,  basalt,  and  quartzite  show  the  de- 
pendence of  stick-slip  quantities  on  temperature,  yet  Heard  (1960)  shows 
that  the  brittle  shear  strength  of  a limestone  similar  to  that  tested  is  not 
affected  within  the  temperature  and  nominal  normal  stresses  used  in  this 
experiment.  Similarly,  a typical  basalt  varied  only  2%  in  its  bulk  modulus 
between  25°  and  200°  (Nafe  et.  al. , 1968), bulk  modulus  being  simply  related 
to  the  elastic  (Youngs)  modulus.  Since  strength  and  elastic  properties  do  not 
vary  appreciably  with  temperature  in  the  range  tested,  the  internal  properties 
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of  the  material  do  not  seem  to  be  involved  in  stick -slip  so  much  as  the  pro- 
perties and  reactions  in  the  surface  and  near  surface;  that  is,  the  asperities 
and  their  roots  and  surroundings.  Plastic  deformation,  which  is  affected 
strongly  by  temperature,  is  regarded  as  the  probable  mode  of  deformation 
in  the  surface. 

In  this  connection  basalt,  unconfined  at  room  temperature,  is  seemingly 
a very  brittle  material.  The  pattern  of  slip,  however,  is  almost  precisely 

that  of  a veiy  ductile  material.  How  can  this  be  explained?  One  factor  is 
alignment  of  the  basalt  crystals  because  of  its  origin  as  a flow.  The  flow  boun- 
daries gave  an  external  lineation  to  the  basalt  which  is  reflected  by  an  internal 
lineation,  at  least  with  the  plagioclase  feldspar  crystals.  Test  series  where 
the  upper  and  lower  blocks  were  arranged  so  the  lineation  was  parallel  gave 
good  stick-slip,  whereas  one  basalt  test  series  had  the  flow  directions  of  top 
and  bottom  blocks  perpendicular  and  gave  poor  stick-slip.  Apparently,  good 
bonding  in  the  basalt  only  occurs  if  the  crystals  are  properly  oriented. 

Secondly,  poor  conductors  of  heat  create  very  high  temperature  rises  upon 
sliding.  Even  at  sliding  speeds  as  low  as  1 or  2 feet  per  second,  tempera- 
tures tend  toward  the  melting  points  at  surfaces  in  contact  (Bowden  and  Tabor, 
1954).  Basalt,  with  a thermal  conductivity  of  4.1  to  6.0  x 10  3 cal/cm-sec  °C, 
is  a poor  thermal  conductor  (Nafe  and  Drake,  1968).  In  metals  most  of  the  heat 
is  diffused  into  the  bulk  of  the  metal,  but  rock  heat  created  upon  sliding  is  slow 
to  diffuse,  and  the  contact  surface  may  reach  the  melting  point  of  plagioclase 
feldspar,  1090°C  at  one  bar  pressure,  but  as  low  as  750°C  at  5000  bars  (Turner 
and  Vcrhoogen,  i960).  /»t  these  temperatures  of  the  asperity,  plastic  flow  easily 
occurs  and  the  bonding  becomes  stronger.  With  the  ambient  temperature  aiding 
in  the  temperature  rises,  intimate  contact  becomes  even  easier,  and  bonding 
even  stronger.  A partial  melting  and  recrystallization  of  the  plagioclase  feld- 
spar with  each  slip  would  shortly  bring  slightly  misaligned  crystals  into  good 
alignment  with  the  direction  of  slip. 


The  combination  of  a softening  crystal  and  shear  stress  would  result  in 
the  same  phenomenon  in  basalt  as  occurs  in  metals  - junction  growth.  The 
tangential  stress  combined  with  the  plasticity  and  bonding  of  the  crystal  per- 
mits a growth  of  the  area  in  good  bonding  contact  of  each  asperity  in  the  di- 
rection of  the  shearing  stress.  If  the  increase  of  shear  force  is  faster  than 
the  junction  growth,  slip  occurs.  The  presence  of  contaminants  will  also 
inhibit  junction  growth.  As  each  slip  occurs,  more  of  the  asperity'  areas  are 
increased.  Thus  as  the  sliding  continues,  asperity  area  increases,  bonding 
gets  better,  and  the  stick-slip  parameters  become  larger. 
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CONCLUSIONS 


The  tests  conducted  on  a limestone,  a basalt,  and  a quartzite  suggest 
the  following  general  conclusions: 

1.  Stick-slip  may  be  a universal  phenomenon  in  rocks,  having  been  demon- 
strated in  three  widely  divergent  rock  compositions. 

2.  Stick-slip  of  lightly  loaded  rocks  is  due  to  asperity  to  asperity  bonding. 
This  bonding  is  affected  by  temperature  and  is  believed  to  be  due  to  plastic 
deformation  of  asperities.  Junction  growth  of  asperities  may  occur  in  some 
rocks. 

3.  The  amount  and  character  of  stick-slip  relate  to  the  loading  system 
involved  and  the  mineralogical  composition  of  the  rock. 

4.  Temperatures  above  those  at  the  initiation  of  stick-slip  result  in 
larger  slips,  longer  sticks,  and  higher  shear  load  relaxations.  The  rates  of 
change  with  increasing  temperature  depends  on  the  material. 

5.  When  smooth  slip  occurs,  it  will  usually  occur  at  temperatures  below 
100°C.  Smooth  slip  at  low  temperature  is  due  to  a boundary  layer  of  adsorbed 
water.  As  the  adsorbed  water  is  driven  off,  stick-slip  is  induced.  Alternatively 
at  low  temperatures  the  Rebinder  effect  may  cause  a softening  of  the  surface  by 
water,  leading  to  creep  of  asperities  and  increased  bonding.  This  results  in 
stick-slip  at  low  temperatures. 

6.  Brittle  fracture  and  the  ^cumulation  of  debris  both  act  to  reduce  the 
amount  of  stick-slip  at  a given  temperature,  by  reducing  the  area  of  bonding  of 
asperit  ies  and  by  producing  weaker  bonds  in  the  system.  Debris  separates  the 
crystals  which  must  be  aligned  for  good  bonding.  Smooth  slip  may  occur  below 
100°  in  blocks  covered  with  debris  in  the  same  manner  as  in  initially  clean 
blocks. 
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APPENDIX  A. 

ACTIVATION  ENERGY  FOR  CREEP  ALONG  JOINTS  OF  LIMESTONE 

Introduction 

In  order  to  allow  movement  of  one  material  past  another,  the  bonds  con- 
necting the  materials  must  be  broken.  If  the  Bowden  and  Tabor  theory  of 
adhesion  friction  {1964)  is  correct,  the  friction  of  rock  joints  should  be  con- 
trolled by  bonding  of  asperities  across  the  joint.  If  the  materials  then  slide 
slowly  past  each  other,  or  creep,  the  rate  of  such  movement  should  depend 
on  the  rate  of  making  and  breaking  of  bonds.  This  rate  can  be  investigated 
using  the  postulates  of  the  rate  process  theory  first  set  forth  by  Eyring  in 
Glasstone  et.  al.  (1941)  for  chemical  solutions  and  later  extended  to  soils  by 
Noble  and  Demirel  (1969)  and  Mitchell  et.  al.  (1968,  1969).  An  attempt  to 
confirm  this  hypothesis  was  made  using  temperature-controlled  direct-shear 
of  limestone.  The  results  are  only  preliminary  because  of  the  overwhelming 
domination  of  tests  by  stick-slip. 

Theory 

The  investigation  of  creep  of  limestone  along  joints  centered  on  use  of  the 
rate  process  assumption  to  determine  the  applicability  of  the  rate  process 
equation  (Noble  and  Demirel,  1969): 

ln6  = In  A - ^-+  p(T  - ^vpn)  (A-l) 

where  5 is  creep  rate,  $1  is  activation  enthalpy,  k is  Boltzmann’s  con- 
stant, T is  absolute  temperature  during  shear,  r = p'AT  where  p'  is  flow 
unit  volume  and  p'r  represents  the  energy  applied  by  the  shearing  stress  r» 
yu.  is  a coefficient  representing  the  effect  of  dilatancy,  pR  is  normal  stress,  and 
A is  a constant  including  the  effects  of  structure  (arrangement  of  crystals  and 
grains  of  the  rock)  and  of  entropy.  Written  in  exponential  form: 

-/U} 

6 = A exp(-gr)  exp  (p  (r- y«pn)  ).  (A-2) 
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The  dilatancy  of  the  smoothly  ground  limestone  was  not  large  enough  to  be 
measured  with  the  measuring  device  available.  Thus,  under  constant  shear 
stress  conditions  the  second  exponential  term  may  be  included  in  the  constant 
coefficient  giving  the  following  equation: 


ft  = x exp( 


(A-3) 


where  R = N k;  N is  Avogadro's  number  and  jH  is  considered  as  activation 
a a 

energy  per  mole. 

Assuming  that  X is  independent  of  temperature  and  that  material  struc- 
ture is  a constant  included  in  X,  activation  energy  can  be  evaluated  by 
determination  of  the  rate  of  deformation  at  different  temperatures. 


ln6l-ln52  -/»! 
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The  analysis  of  the  relative  size  of  flow  unit3  and  thus  the  relative  density 
of  bonds  is  possible  in  tests  conducted  at  a constant  temperature.  Bond  density 
is  analyzed  by  a jump-creep  test  which  resembles  a Dorn  test  with  shear  load 
suddenly  increased  instead  of  temperature.  The  rate  of  creep  decreases  with 
time  in  an  approximately  logarithmic  manner  (Mitchell  et.  al. , 1969).  Even 
relatively  long-term  "steady-state”  creep  is  actually  transient  in  rate.  These 
decreasing  creep  rates  can  be  analyzed  according  to  methods  developed  by 
Singh  and  Mitchell  (1968)  and  described  as  a stress-strain-time  function  for 
soils.  The  equation  described  by  these  authors  is: 

m , \ m 
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in  which  e is  strain  rate,  for  which,  as  dimensions  do  not  change,  defor- 
♦ 

mation  rate  5 may  be  substituted,  t.  is  any  time  during  application  of  the 
initial  shear  stress  rx,  t'  is  the  time  after  commencement  of  the  test  of 
application  of  the  stress  rP , and  t . is  the  time  at  which  the  strain  rate 

i 1 

e(T-,  t.)  is  recorded  under  stress  ra . The  strain  rates  thus  obtained  can 
J 2 

be  related  by 


-/tfl 

c - A exp(-JJr-)  exp(2aT) 


(A -8) 


For  comparison  with  the  rate  process  equation  used  in  this  investigation: 
a - s/2  (A -9) 


~ = 4a  kT  - 26’ 


(A-lO) 


where  the  > component  is  the  distance  between  equilibrium  positions  of 
bonds,  and  S is  the  number  of  flow  units  per  unit  area.  The  size  of  c 
should  decrease  with  increasing  number  of  bonds  in  the  same  material. 
Therefore,  s',  which  has  volume  units,  should  be  a relative  measure  of  the 
number  of  bonds. 


Testing  Procedure 

The  basic  test  apparatus  is  the  same  as  that  described  in  the  body  of  this 
report.  The  exception  is  the  mechanism  for  inducing  the  shear  load.  Since 
the  creep  tests  were  performed  as  constant  stress  type  tests,  the  shear  load 
mechanism  was  simply  a wire  running  over  a low-friction  pulley  with  weights 
attached  to  the  wire. 

The  specimens  were  prepared  in  the  same  manner  as  described  for  the 
stick-slip  tests.  After  preparation  all  blocks  were  raised  to  200°C  in  the 
testing  apparatus  with  the  normal  load  applied  for  four  hours.  The  temperature 
was  then  decreased  to  the  initial  testing  temperature.  This  was  done  so  that 
all  samples  would  have  the  opportunity  to  initially  form  the  same  number  of 
bonds. 
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In  order  to  check  for  full-shear  load  and  prevent  rapid  movement  upon  ini- 
tial creep  loading,  the  upper  block  was  initially  loaded  by  slowly  adding  lead 
shot  to  the  load  device  until  full  shear  occurred.  Full  shear  was  defined  as 
movement  of  the  upper  block  across  the  lower  block  with  continuous  accelera- 
tion. After  the  full-shear  load  was  recorded,  any  debris  (usually  none,  in 
these  cases  with  light  normal  load ) was  wiped  off  with  a clean  cloth.  The  rock 
was  then  repositioned,  re-heated  to  200°,  and  cooled  to  test  temperature.  A 
shear  load  less  than  that  required  for  full  shear  was  applied.  Lead  shot  was 
then  slowly  added  until  creep  was  induced. 

Early  in  the  testing  program  it  was  recognized  that  dilatancy  due  to 
structural  differences  between  blocks  made  it  impossible  to  validly  compare 
deformation  rates  on  one  block  with  those  on  another  block.  It  was  decided 
that  valid  data  could  be  obtained  by  use  of  the  Dorn  temperature  jump  method, 
originally  developed  for  creep  studies  of  metals  (Mitchell  and  Singh,  1969). 

In  the  Dorn  method,  the  specimen  is  allowed  to  come  to  a steady  state  of 
deformation  rate  at  a lower  temperature  T^.  The  temperature  is  then  raised 
as  rapidly  as  possible  to  a higher  temperature  T whereupon  the  strain  rate 
increases  to  a second  steady  state  deformation  rate  5g.  When  these  values 
are  substituted  into  equation  (A -4),  the  activation  energy  is  obtained. 

Results  and  Interpretation 

Calculation  of  activation  energies  by  the  Dorn  method  in  individual  blocks 
rapidly  showed  that  activation  energies  were  not  temperature  independent.  Re- 
sults are  shown  in  Table  A-l.  In  the  table  the  temperature  ranges  represent 
zones  within  which  the  lower  and  higher  temperature  of  the  jumps  are  included. 
Individual  jumps  might  have  different  beginning  and  ending  temperatures  within 
the  ranges.  The  standard  deviation  as  a percentage  of  the  mean  should  be 
noted.  This  percentage  is  low  in  the  temperature  areas  below*  100°C  quite 
irregular  in  the  100  - 150°  range,  and  decreases  again  at  higher  temperature. 
To  insure  that  activation  energies  do  rise  with  temperature,  and  do  not  merely 
reflect  the  size  of  the  temperature  jump  employed,  the  activation  energy  per 
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Deformation 


Figure  A-l.  Dorn  temperature  jump  method  for  determination 
of  activation  energy. 
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Table  A-l 


Activation  energy  values  for  creep  along  limestone  joints 


Test  Normal  stress  Activation  energy  - kilocalories/mole 

kg/cm“  Temperature  ranges 

35-100°  100-150°  140-200° 


29 

0.375 

-9.226 

33B 

0.751 

-10.  73 

34 

0.751 

-10.05 

-36.43 

37 

0.751 

-5.12 

40 

0.751 

-26. 99 

42 

1.220 

-18. 098 

43 

1.220 

-5.369 

-6.487 

-13.917 

46 

1.220 

-8.146 

49 

1.877 

-8.779 

-5. 590 

50 

1.877 

-10.599 

-15. 983 
-16.408 

52 

0.188 

-8. 709 

-27.  503 

54 

1.877 

-5.279 

-21.474 

55 

0.188 

-12.525 

56 

0.188 

-37.045 

64 

0.751 

-6.466 

69D 

0.751 

-5.  742 

-28.651 

69E 

0.751 

-24.673 

69F 

0. 751 

-11.108 

Mean 

-7.540 

-11.664 

-23. 014 

Standard  deviation 

1.991 

7.129 

9.849 

Standard  deviation  as 

26.4 

61.1 

42.8 

percentage  of  mean 
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degree  of  temperature  rise  during  each  jump  was  calculated,  averaged  according 
to  the  range  of  the  jump,  and  is  presented  in  Table  A-2.  The  meaning  of  the 
temperature  ranges  if  as  noted  above. 


Table  A-2.  Change  in  activation  energy  per  degree  of  temperature  jump 


Temperature 

range 

Mean  activation  energy 
(kilocalories/mole/0  C) 

Standard  deviation 

35  - 100° 

0.184 

0.  067 

100  - 150° 

0.281 

0.188 

140  - 200° 

0.  66-1 

0.206 

This  demonstrates  conclusively  the  dependence  of  activation  energy  upon 
temperature. 

Ar.  attempt  was  made  to  determine  if  the  temperature  dependence  could  be 
expressed  by  a pre-exponential  term,  expressing  deformation  rate  as: 

* - BTm  exp  (-^fi/HT)  (A -11) 

Constant  H and  m values  were  not  obtained.  This  is  probably  due  to  the  im- 
portance of  the  structure  factor  in  It,  and  the  inability  to  separate  structure 
factor  from  the  items  strictly  dependent  upon  temperature. 

The  values  of  activation  energies  obtained  seem  to  be  indicative  of  different 
natures  of  sliding  surfaces  at  different  temperatures.  The  activation  energy 
\nlue  at  temperature  below  100°,  -7.51  kilocalories  per  mole,  is  very  similar 
to  values  for  water,  -4  to  -5  kilocalories  per  mole,  determined  by  Glasstone 
et.  al.  (1941).  Values  of  activation  energy’  for  temperature  between  140°  and 
200°  is  similar  to  the  values  of  Mitchell  et.  al.  (i960)  for  dry  quartz  sand. 
These  observations  lead  to  a conclusion  that  below  100°C  the  creep  of  limestone 
is  actually  proceeding  on  a film  of  adsorbed  water.  Above  140°  true  calcite-to- 
calcite  bonds  are  obtained  as  all  adsorbed  water  is  driven  off.  An  intermediate 
condition  may  prevail  at  intermediate  temperatures  until  all  water  has  had  time 
to  be  removed.  Chemisorbed  or  tightly  physically  adsorbed  water  may,  indeed, 
remain  in  places,  perhaps  where  asperities  are  not  actually  in  contact.  That  the 
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activation  energy  value  below  i00°  is  above  that  of  water  may  be  due  to  higher 
asperities  physically  squeezing  out  the  water  and  permitting  actual  calcite- 
calcite  bonds  in  limited  places.  It  has  also  been  speculated  that  the  activation 
energy  higher  than  that  of  liquid  water  may  be  due  to  a structural  arrangement 
of  the  adsorbed  water  molecules  which  is  different  from  that  of  liquid  water. 

Five  tests  conducted  at  100°  proved  amenable  to  analysis  of  the  relative 
size  of  flow  units  by  the  method  described  in  Equations  (A-5)  through  (A-10). 
average  values  of  c obtained  for  the  limestone  are  shown  in  Table  A-3. 


Table  A-3.  Flow  unit  volumes  of  limestone 


Normal  load 
(kilograms) 

Normal  ^ress 
(kg/cm  ) 

Flow  unit  volume 
(cm  ) 

100 

1.88 

5.493  x 10"19 

40 

0.75 

7.  506  x 10"19 

10 

0.188 

592. 0 x 10~19 

These  values  do  tend  to  increase  with  decreasing  normal  load.  Thus, 
the  volume  occupied  by  a single  bond  must  decrease  with  increasing  load. 

This  indicates  a relative  increase  in  the  number  of  bonds.  The  limited  data 
base  for  this  observation  must  be  emphasized. 

Conclusions 

The  following  conclusions  can  be  drawn  from  the  creep  tests  on  limestone: 

1)  Activation  energy  experiments  show  that  $1  is  not  independent  of 
temperature  but  increases  with  increasing  temperature.  This  does  not  refute 
the  theoretical  equations  for  activation  energy,  but  indicates  a change  in  the 
strength  or  nature  of  the  chemical  bonding  with  temperature. 

2)  The  values  of  these  strengths  related  to  the  test  conditions  and  the 
temperatures  involved  lead  to  the  conclusion  that  activation  energies  below  100°C 
reflect  creep  occurring  in  a thin  film  of  adsorbed  water.  Average  activation 
energies  are  above  that  of  water  by  a small  amount.  This  is  believed  to  reflect 
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either  occasional  ealcite-calcite  bonds,  or  to  be  a reflection  of  possible  higher 
activation  energy  of  adsorbed  water  compared  to  that  of  liquid  water.  Above  140°C, 
the  large  activation  energy  is  believed  to  reflect  true  calcite-calcite  bonding. 

3)  Limited-data  permitting  calculation  of  flow  volume  units,  which  are  propor- 
tional to  the  number  of  bonds,  show  that  the  number  of  bonds  increases  with  in- 
creased normal  load. 
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INTRODUCTION 


Statement  of  the  Problem 

The  purpose  of  this  research  is  to  study  and  formulate  a failure  mechanism 
which  could  apply  to  movement  of  rock  bodies  along  their  joints  and  fracture 
planes.  Rock  bodies  can  acquire  these  features  of  layering,  jointing,  and  frac- 
turing either  as  a part  of  their  geological  origin  or  from  external  causes  such  as 
stress  release,  earthquakes,  blasting,  or  other  engineering  activities.  Such 
features  in  rock  bodies  are  centers  of  weakness,  and  control  strength. 

Observed  failures  involving  shear  within  rock  or  soil  masses  often  occur 
in  plane  strain  , i.  e.  with  a relatively  constant  cross  section  normal  to  the 
failure  surface.  Such  a two-dimensional  failure  can  be  modelled  by  plane- 
strain  shear  tests,  whereby  the  rock  or  soil  material  is  confined  to  prevent 
strain  in  the  third  dimension.  In  the  present  study,  such  strain  was  prevented 
not  by  confinement  but  by  the  geometry  of  the  individual  particles,  which  are 
rods  lying  parallel  to  one  another  in  a shear  box. 

Scope 

This  research  involves  subjecting  ideal  assemblages  of  rods  to  a biaxial 
stress  field,  varying  the  principal  stress  ratio,  and  recording  the  stresses, 
strains,  volume  changes,  and  translocations  and  rotations  of  individual  members 
of  the  assemblage.  Simultaneously,  a theory  was  developed  to  predict  the  be- 
havior, such  that  the  theory  could  be  tested  and  revised  as  necessary. 

The  biaxial  testing  method  used  has  the  following  advantages  over  the  tests 
of  natural  materials. 

1.  A continuous  photographically  recorded  inventory  of  individual  particle 
movements  in  relation  to  assemblage  stress,  strain  and  volume  changes. 

2.  Ready  prediction  of  an  ideal  behavior,  such  that  departures  from  the 
ideal  can  be  identified  and  related  to  causes. 
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3.  Controlled  test  conditions  which  minimize  the  sampling  variable 
by  reusing  the  same  sample. 

4.  Ability  of  testing  different  regular  geometrical  array,  affording 
a precise  variation  of  void  ratio  and  packing  density. 

5.  An  accurate  visualization  of  failure  mechanism,  in  that  failure 
by  sliding  and/or  rotation  can  be  differentiated. 
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BACKGROUND  AND  LITERATURE  REVIEW 


In  order  to  formulate  a failure  theory  for  a granular  mass  subjected  to 
a biaxial  stress  field,  it  is  necessary  to  understand  the  past  and  present 
concepts  of  dilatancy  and  sliding  friction. 


D5‘„-  ■'anoy 

Dilatancy  — "The  property  of  granular  masses  of  expanding  in  bulk  with 
change  of  shape.  It  is  due  to  the  increase  of  space  between  the  individually 
rigid  particles  as  they  change  their  relative  positions"  (Century  Dictionary). 

Dilatancy  was  first  described  and  named  by  Professor  Osborne  Reynolds 
(1885),  who  snowed  that  a dense  sand  mass  expands  upon  shearing.  He 
observed: 

...  I would  point  out  the  existence  of  a singular  fundamental 
property  of  such  granular  media  which  is  not  possessed  by 
known  fluids  and  solids.  On  perceiving  some  thing  which  resem  - 
bles nothing  within  the  limits  of  one's  knowledge,  a name  is  a 
matter  of  great  difficulty.  I have  called  this  unique  property  of 
granular  masses  "dilatancy",  because  the  property  consists  in  a 
definite  change  of  bulk,  consequent  on  a definite  change  of  shape 
or  distortional  strain,  any  disturbance  whatever  causing  a change 
in  volume  and  generally  dilation. 

Reynolds  observed  that  with  granular  media,  so  long  as  the  grains  are 
held  in  mutual  equilibrium  by  stresses  transmitted  through  the  mass,  every 
change  of  relative  position  of  the  grains  is  attended  by  a consequent  change  in 
volume;  and,  if  in  ar.>  way  the  volume  is  fixed,  than  all  change  of  shape  in 
prevented.  Professor  Reynolds  made  one  assumption,  that  the  position  of 
any  internal  particle  becomes  fixed  if  the  positions  of  surrounding  particles 
are  fixed.  Without  frictional  bridging  of  particles  over  voids,  this  condition 
would  always  be  fulfilled.  It  follows  that  no  grain  in  the  interior  can  change 
its  position  in  the  mass  by  passing  between  the  contiguous  (touching)  grains 
without  disturbing  them;  hence,  whatever  alteration  the  medium  may  undergo, 
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the  same  particle  will  always  be  in  the  same  neighborhood.  If,  then,  sueh 
a medium  is  subjected  to  internal  strain,  the  shape  of  each  elementary  group  is 
determined  by  shape  and  arrangement  of  surrounding  particles.  Any  distor- 
tion of  the  boundaries  of  such  a medium  will  cause  distortion  of  the  arrange- 
ments of  its  particles,  necessitating  a change  in  volume  and  hence,  the  mean 
density.  If  particles  are  rigid,  the  relations  between  distortion  and  dila- 
tancy  are  independent  of  friction;  that  is  to  say,  the  same  distortion  of  any 
bounding  surface  must  mean  the  same  internal  distortion  whatever  the  friction 
may  be.  The  only  possible  effect  of  friction  is  to  render  the  grains  stable 
under  circumstances  in  which  they  would  not  otherwise  be  stable. 

Mead  (1925),  while  applying  Reynolds'  theory  of  dilatancv  to  solid  rocks, 
concluded  that  incoherent,  granular  masses  such  as  sand,  in  a condition 
approaching  maximum  density  packing  (rhombic  packing),  are  dilated  by 
deformation;  whereas,  in  a condition  of  open  packing  (cubic  packing),  they 
deform  without  di'ation.  Prevention  of  free  dilation  by  enclosing  pressures 
induces  failure  by  fracture  or  shear  when  the  mass  is  deformed,  and  with 
development  of  joints  and  faults  along  thin  zones  of  dilation.  This  manner 
of  failure  icquires  a minimum  increase  in  volume  and  involves  dilation  only 
in  the  sheai  zone.  When  free  dilation  i.°  not  prevented,  the  granular  mass 
deforms  by  flow.  Such  a deformation  of  closely  packed  grains  involves  the 
entire  mass  and  causes  a much  greater  volume  increase  than  that  required  by 
failure  along  definite  shear  planes,  and  can  be  called  plastic  deformation 
or  plastic  flow. 

Performing  an  experiment  somewhat  similar  to  Reynolds'  experiments 
with  a rubber  balloon,  sand  and  water,  Mead  (1925)  observed  that  if  the 
amount  of  fluid  in  the  granular  aggregate  is  only  sufficient  to  fill  the  voids  in 
the  condition  of  maximum  density  packing,  deformation  of  the  mass  requires 
an  increase  in  volume,  and  the  mass  fails  largely  by  fracture  and  not  by 
plastic  deformation.  If  the  fluid  available  is  slightly  in  excess  of  this 
amount,  the  aggregate  is  easily  deformed  plastically  up  to  a point  where  the 
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increased  voids  absorb  the  available  fluid,  at  which  point  the  mass  becomes 
rigid.  On  the  other  hand,  when  the  available  fluid  phase  is  sufficient  to  fill 
the  voids  with  grains  arranged  in  a cubic  or  minimum  density  packing,  the 
m..ss  may  be  deformed  to  any  extent  without  an  increase  in  volume,  and 
failure  may  be  due  to  plastic  flow  offering  very  little  resistance  to  defor- 
mation. 

The  mechanics  of  response  to  deformation  of  incoherent  granular  masses 
can  be  applied  to  solid  rocks  by  conceiving  of  them  as  having  a solid  and, 
potentially,  an  ideally  incompressible  fluid  phase.  The  latter  may  cause 
the  rock  to  yield  to  deformation  by  flowing  or  by  fracture.  In  simplest  terms, 
the  rocks  may  be  regarded  as  granular  aggregate,  the  hard  grains  the  solid 
phase.  To  the  extent  that  the  rocks  are  porous,  the  pores  represent  the  volume 
of  fluid  phase  and  the  material  that  occupies  the  pores  is  the  fluid  phase. 
Alternately,  the  solid  phase  can  be  represented  by  harder,  more  resistant 
minerals  of  rock  mass,  whereas  the  potentially  fluid  phase  can  be  represented 
by  those  constituents  of  rock  which  are  relatively  mobile,  as  evidenced  by 
their  rearrangement  to  schistose  (foliated,  i.e.  capable  of  splitting  up  in 
thin  irregular  plates)  structures  through  processes  of  crystallization  and 
and  recrystallization.  Thus,  potential  fluid  phases  occur  under  certain  con- 
ditions of  composition,  pressure,  temperature,  and  rat  e of  deformation. 

Brown  and  Hawksley  (1947),  while  experimenting  with  regular  two- 
dimensional  packings  of  equal  spheres,  observed  a marked  tendency  for 
the  spheres  to  move  together  in  groups;  thus,  they  found  a third  process 
(besides  the  appearance  of  slip  lines  and  dilatant  expansion)  by  which  the 
tight  regions  in  the  array  break  down.  They  pointed  out  that  relative  move- 
ment between  groups  of  tight  packing  leads  to  the  loss  of  the  peripheral 
spheres  by  a sort  of  "abrasion"  between  the  groups  which  gave  loose 
irregular  packing  between  the  groups. 

In  effect,  their  observation  showed  that  there  is  a tendency  for  local 
regions  in  .packings  to  become  or  to  remain  tightly  packed.  In  general, 
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the  tight  regions  did  not  dilate  uniformly,  but  either  failed  locally  in  a line 
of  slip  or  moved  as  a group.  An  initially  uniform  loose  packing  collapsed  lo- 
cally, while  an  initially  tight  packing  failed  irregularly  along  lines  of  slips, 
in  either  case  giving  a nonuniform  distribution  of  voids.  Groups  of  tight 
packing  were  not  found  to  mesh  with  each  other. 

They  concluded  that  in  tight  regular  arrangements,  dilatancy  is  a geome- 
trical necessity  if  deformation  is  to  occur;  whereas,  in  irregular  arrangements, 
dilatancy  does  occur,  but  the  explanation  lies  in  the  "stability"  of  the  packing. 
Moreover  since  the  groups  of  tight  packing  have  been  observed  to  move  as  a 
whole  without  dilation,  or  when  they  do  dilate  to  fail  along  the  line  of  slip 
(in  preference  to  a uniform  dilation),  it  would  seem  that  these  groups  possess 
some  rigidity.  Thus  an  assemblage,  regarded  as  an  arrangement  of  funda- 
mental units  which  are  semirigid,  may  dilate  through  the  interplay  of  these 
units.  Since  there  cun  be  an  interchange  of  particles  between  adjacent  units, 
such  an  interplay  is  more  "flexible"  than  the  interplay  between  individual 
particles. 

Andrade  and  Fox  (1949),  while  working  with  a two-dimensional  regular 
(hexagonal)  array  of  ebonite  and  polyethylene  rods,  pointed  out  that  the  dilation 
of  a regular  array  is  intimately  connected  with  the  mechanism  of  deformation. 
They  detected  two  classes  of  dilational  deformations,  both  of  them  conse- 
quences of  slip  on  well-defined  planes.  In  the  first,  the  dilation  was  localized 
at  the  boundaries  of  regular  areas  of  rods  which  preserved  their  original  packing; 
in  the  second,  it  was  more  or  less  irregular  throughout  some  of  these  areas. 

The  first  type  is,  in  a sense,  the  more  fundamental  since  the  deformation 
necessarily  involves  slip  which,  in  turn,  involves  the  primary  dilation.  The 
occurrence  of  secondary  dilation  was  ascribed  to  some  extent,  at  least,  to 
surface  cohesion  of  the  array.  In  the  absence  of  secondary  dilation,  the 
primary  dilation  showed  a fairly  regular  alternation  corresponding  to  the 
slipping  and  healing  of  the  array.  The  regularity  of  this  alternation  was  proba- 
bly disturbed  by  the  secondary  dilation  since  secondary  dilation  could  not  be 
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expected  to  show  any  regular  variations.  Secondary  dilation  was  also  found  to 
he  dependent  on  internal  friction  - the  smaller  the  internal  friction,  the  smaller 
the  role  played  by  secondary  dilation.  The  result  was  that  greater  friction 
gave  dilation  which  was  not  as  regular  as  it  was  in  the  case  of  less  friction. 
More  surfaces  of  misfit,  some  of  which  became  healed  in  the  progressive 
stages,  were  observed  with  greater  friction.  This  fact  could  be  held  as  the 
direct  consequence  of  secondary  dilation  and  increased  internal  friction. 

Apparently  internal  friction  stabilizes  the  array  against  the  onset  of  slip; 
but,  when  frictional  forces  arc  eventually  overcome,  the  deformation  becomes 
catastrophic  with,  consequentially,  a greater  tendency  for  the  areas  of  rods 
to  break  up  than  when  the  deformation  proceeds  gradually,  as  in  the  case  of 
low  friction. 

Hills  (1903)  notes  that  the  relationship  of  dilalaney  to  the  spacing  of  shear 
planes  developed  in  a deformed  rock  may  be  very  important.  Since  a con- 
siderable amount  of  work  is  done  in  developing  a single  shear  plane  under  dila- 
tant  conditions,  the  principle  of  least  work  suggests  that  the  number  of  planes 
will  be  small  in  coarse-grained  materials;  whereas,  with  a fine-grained 
aggregate,  each  shear  plane  involves  less  work  and  the  effort  may  be  expended 
with  equal  facility  on  a greater  number  of  planes. 

Friction  of  a Dilatant  Mass 

The  internal  friction  of  a soil  or  of  any  dilatant  mass,  may  be  considered 
to  be  comprised  of  two  components;  sliding  friction,  and  dilatancy  or  "inter- 
locking". Interlocking  has  also  been  recognized  bv  later  workers  in  metals. 

At  the  present  stage  of  knowledge,  it  can  be  said  that  Reynolds'  dilatancy 
theory  superimposed  on  the  adhesion  theory  of  friction  can  roughly  formulate 
the  friction  theory  of  dilatant  masses.  If  F(]  is  the  force  of  friction  due  to 
dilatancy  (Fig.  1)  and  Fa  is  that  due  to  adhesion,  then  F,  the  force  of  friction 
of  the  dilatant  mass,  is  given  by 
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F = Fd  + Fa 

- VV  tano  + K VV 

= W(K  + tano) 

F - (K  + tano) 

VV 


whore  JF  can  be  called  the  coefficient  of  friction  of  dilatant  masses,  K is  a 

w 

constant  reflective  of  adhesion  or  sliding  friction,  and  to  is  the  angle  of  sliding, 


TWO  S URFACES 


W sin  0 


Dilatancy  model:  F{j  - VV  tano 


F cos  <f 


CONTACT  AREA  = A 


Adhesion  Theory.  F0  - KW 

a 


Figure  1.  Two  components  of  internal  friction  of  dilatant  masses. 
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Energy  Theory  of  Dilutant  Masses 
Reynolds'  energy  concept 

Taylor  (1948)  is  usually  credited  with  being  the  first  to  attempt  separation 
of  two  components  of  friction  of  dilatant  masses  by  evaluating  the  work  done  by 
dilatant  expansion  in  direct  shear  against  the  normal  pressure,  but  it  seems 
that  Reynolds  (1885)  had  a similar  energy  concept  in  mind  and  which  he  applied 
to  granular  masses,  with  the  difference  that  he  neglected  sliding  friction.  He 
said. 

If  the  particles  were  rigid  the  medium  would  be  absolutely  without 
resilience  and  hence  the  only  energy  of  which  it  would  be  susceptible 
would  be  kinetic  energy,  so  that,  supposing  the  motion  slow,  the 
work  done  upon  any  group  in  distorting  it  would  be  zero.  Thus, 
supposing  contraction  in  one  direction  and  expansion  at  right  angles, 
then  if  px  be  the  stress  in  the  direction  of  contraction,  and  py, 
pz,  the  stress  at  right  angles,  a being  contraction,  b and  c 
expansion, 

px  a 4 py  b + pz  c - 0 
or  supposing  b - c,  py  = pz, 

px  a 4 py(a  + c)  - 0. 

With  friction  the  relation  will  be  different;  the  friction  always  opposes 
strain,  i.e.  tends  to  give  stability. 

Taylor  energy  theory 

Let  be  the  normal  stress,  r the  shear  stress,  r(j  that  portion  of 
shearing  stress  which  acts  to  supply  the  energy  of  expansion  or  shear  stress 
necessary  to  cause  the  sample  to  dilate  against  the  normal  stress,  6 the 
incremental  displacement  in  direct  shear  test  in  the  direction  of  shearing 
torce,  and  by  the  incremental  volume  per  unit  area  which  is  incremental 
displacement  in  the  direction  of  the  normal  stress.  Then 

Work  done  against  normal  stress  - Work  done  by  rd  in  moving  through  6. 

a 

an 

Td  \ 


85 


Shearing  strength  contribution  by  friction  = r - 


T ' rd  r rd 

tan  o 

"n  an  an  ~ 


= tancv  -*•  — 


tano  - tano  + — 
max  1 r 6 


where  is  the  angle  of  friction  of  the  granular  mass  and  g*r  is  the 
internal  coefficient  of  friction.  The  above  energy  equation  was  later  adapted 
to  the  triaxial  test  by  Skempton  and  Bishop  (1950). 

Newland  and  Allely  Theory 

Newland  and  Allely  (1957)  considered  that  the  direction  of  sliding  of 
one  particle  relative  to  another  is,  in  general,  inclined  at  an  angle  to  the 
direction  of  the  applied  shear  stress.  By  resolving  forces  and  displacements 
of  particles  parallel  and  perpendicular  to  their  sliding  surfaces,  they  ar- 
rived at  the  following  expression  for  the  direct  shear  test. 


tan  <p 


Tmax 


= tan(0f  + 0) 


where  tang  = (5  /5  ) 

v //max 


and  tan<^.  = 

where  9 is  the  angle  of  sliding  surface  with  the  direction  of  shear  stress; 
Oj.  is  the  angle  of  internal  friction,  which  accounts  for  the  influence  of  not 
only  the  coefficient  of  sliding  friction,  but  also  of  the  mode  of  failure; 
t'jj  is  the  stress  required  to  overcome  frictional  force,  assuming  sliding 
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planes  are  parallel  to  the  direction  of  the  shear  stress;  ^ is  the  normal 

stress;  and  6/5  if  the  rate  of  volume  expansion, 
v a 

They  suggested  that  using  experimental  values  of  r and  5/5 

J **  max  n v & max 

in  the  above  equations,  t'^  can  lx1  obtained  and  in  turn,  compared  to  experimen- 
tal residual  stress  r 


Newland  and  Allely  also  extended  their  analysis  to  the  triaxial  test. 

Newland  and  Allely  compared  their  equation  with  that  of  Taylor  and 
Bishop,  and  considered  that  the  difference  between  the  energy  approach 
method  and  their  method  lies  in  the  tacit  assumption  that  the  work  done  in  over- 
coming 5Wf  is  the  same  at  both  the  peak  and  the  residual  states;  whereas, 
according  to  their  theory,  the  shear  stress  has  a component  normal  to  the 
plane  of  sliding  which  contributes  to  frictional  strength  but  which  decreases 
in  magnitude  as  expansion  occurs.  In  other  words,  the  work  done  in  over- 
coming frictional  forces  is  greater  at  the  point  of  maximum  shear  stress 
than  at  residual  stress  state. 

Caquot  (1934)  derived  the  following  expression  relating  angle  of  friction 

at  constant  volume  and  true  angle  of  friction  c between  mineral 
cv  u 

surfaces; 

tan  e>cv  -iff  tan  ou 

Bishop  (1954)  also  developed  an  approximate  solution  in  the  form 


15  tan  <nu 
10+3  tano 


While  the  precise  validity  of  these  equations  is  doubtful,  experimental 
data  fill  the  above  equations  closely. 

Rennie  (1959)  studied  the  least  stress  ratio  which  will  cause  failure  in 
a close-packed  face-centered  packing  and  obtained  the  approximate  solution 
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’’JV! 


— '-  - 2 + 2v/TTu  + (higher  terms  in  u) 
a3 

Thurston  and  Deresiewicz  (1959)  also  considered  a face-centered  array 
of  equal  spheres  with  = <7  and  arrived  at  the  following  expression: 

M 4) 


a 

_1  _ /7Tj_8u_ 
<73  /T>  + 4u 

where  u is  a coefficient  of  friction. 


Rowe  Stress-Dilatancy  Theory 

Rowe  considered,  experimentally  and  theoretically,  the  behavior  of 
assemblies  of  cohesionless,  uniform  rods  in  a parallel  stack  and  spherical  par- 
ticles of  uniform  size,  arranged  initially  in  regular  arrays.  The  assemblies 
of  particles  are  subjected  to  axially  symmetrical  state  of  stress.  From  his 
analysis  based  on  a consideration  of  the  forces  between  particles,  he  arrived 
at  the  following  findings  for  regular  packing: 

1.  Whatever  the  geometrical  arrangements  of  solids,  the  stress  ratio 
at  the  peak  strength  and  during  subsequent  states  of  deformation 
follows  the  law 


Oj  - tan  a tan(ou  + S). 

The  energy  ratio  F for  a fixed  orientation  of  particle  movement 
is  given  by  the  expression 


E = 


tan(<pu  + £) 
tan? 


a,(l  + -7) 

e.. 


3.  Slip  occurs  well  past  the  peak  stress  ratio  at  failure,  thus 
establishing  that  the  slip  plane  is  not  the  cause  but  the  result 
of  failure. 
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In  the  above  findings,  a is  the  angle  which  the  imaginary  plane  of  particle 
interlocking  makes  with  the  direction  of  the  minor  principal  stress,  ? is  the 
angle  of  deviation  of  the  tangent  at  the  contact  point  from  the  direction  of  major 
principal  stress,  and  <*u  is  the  angle  of  sliding  friction. 

Howe  extended  his  analysis  to  a random  mass  of  irregular  particles  by 
suggesting  that  the  form  of  law  which  applies  to  any  individual  packing  may  be 
expected  to  apply  to  mixture  of  packings.  He  observed  that  the  angle  a 
characteristic  of  the  particular  packing  arrangements,  disappeared  in  the  expres- 
sion for  E;  thus,  he  concluded  that  the  energy  equation  applies  to  random  as- 
semblages of  particles,  as  well  as  to  regular  arrays.  Since  in  random,  the 
values  of  = vary  initially  throughout  the  mass,  the  angle  of  sliding  must  also 

be  determined.  He  then  shows  that  there  exists  a critical  angle  of  - 

„ n i -n 

E = t an-  ( I + 2 %)  

i dv 
"3  v v . 

*1 

Departure  from  the  stress  dilatancy  behavior  given  by  the  energy  ratio 
equation  is  explained  by  Rowe  in  terms  of  a process  he  refers  to  as  rear- 
ranging. In  a loose  sand,  and  in  an  initially  dense  sand  when  it  reaches  o point 
near  the  maximum  stress  ratio  o1/oP,  the  length  of  *hc  individual  slide  paths  of 
one  particle  over  another  becomes  appreciable  with  respect  to  the  dimensions 
of  the  particles;  hence,  according  to  Howe,  sliding  is  no  longer  restricted 
to  a value  of  = which  gives  the  minimum  rate  of  expenditure  of  energy  in 

internal  friction.  The  effective  value  of  o then  becomes  nr,  where  > o • 

u f f u 

At  the  ultimate  state  of  deformation  when  the  sample  reaches  the  state  at 
zero  rate  of  volume  change,  the  effective  value  of  ou  at  this  ultimate 
constant  volume  condition,  denoted  by  , is  obtained  from  the  equation 

ai 

— - tan*  (l/4ir  + 1/2<0CV), 
by  using  the  observed  stress  ratio. 
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Rowe  (1963)  applied  {he  stress  dilntancy  theory  to  the  stability  of  earth 
masses  behind  retaining  walls,  in  slopes,  and  in  foundations. 

The  stress  dilataney  theory  of  granular  masses  postulated  by  Rowe  was 
discussed  by  Scott  et  aL  t'19G4),  Gibson  and  Morgenstern  (Trollope  et  al. , 
lfM»:t),  Scott  (19(53),  and  Trollope  ct  al.  (1963).  Their  main  criticism 
was  directed  towards-. 

1.  t he  assumed  mechanism  of  deformation, 

2.  the  assumed  absence  of  tolling, 

3.  tne  assumption  that  the  energy  ratio  K is  minimum  in  a random 
assembly,  and 

•1.  the  significance  of  the  o~plane. 

Horne  confirmation 

A more  general  derivation  was  presented  by  Horne  (1965)  who  did  not 
restrict  his  analysis  to  ideal  packing.  He  obtained  exactly  similar  results 
to  that  of  Rowe,  and  thus  substantiated  his  theory.  Horne  built  up  his  ana- 
lysis on  the  basis  of  the  foil  owing  hypotheses; 

1.  The  assembly  consists  of  rotund,  rigid,  cohesionless  particles 
with  a constant  coefficient  of  sliding  friction.  Elastic  and  plastic 
deformation,  crushing,  and  cracking  are  all  ignored. 

2.  Deformation  occurs  as  a result  of  relative  motion  between  groups 
of  particles.  Motion  is  not  facilitated  by  the  presence  of  individual 
particles  acting  as  rollers  between  groups. 

Horne  obtained  the  expression  for  the  energy  ratio  E by  writing  a virtual 
work  equation  for  the  input  - Then  he  minimized  this  ratio  to  obtain  a 
value  of  r ~ 45  -1m  which  then  led  to 


• • 

For  the  triaxial  compression  test  with  c9  - ^ and  c 0 ~ 1 3.  this 
reduces  to  Howe's  equation.  Horne  thus  established  the  limitation  of  the 
stress  dilatancy  theory  and  concluded  that  the  equation  of  energy  ratio  E 
that  provides  a relationship  between  work  quantities  r.y . and 

does  not  provide  a relationship  between  stress  and  strain  rates  separately. 
He  also  concluded  that  the  relationship  may  not  apply  to  a highly  compacted 
assembly  with  a high  degree  of  interlocking. 
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Tinoeo  ;uid  Handy  Theory 

Tinoco  and  Handy  (19f>7)  considered  a random  particle  assemblage  and 
the  possibilities  for  either  sliding  or  rolling.  They  demonstrated  that 
sliding  occurs  at  contacts  where  the  angle  of  inclination  - is  maximized, 
whereupon 

j,  1 - sinoo 

— tan'  <15  * c '2) 

1 - snio  s 

• s 

which  is  identical  to  the  empirical  Mohr-Coulomb  theory.  Next,  considering 
the  work  of  volume  change,  they  derived  an  equation  for  principal  stresses 
due  to  friction.  In  simplified  form. 


_1  (1  - 


/--.  - (1  - fv  '*e, ) isinc 


where  cs  is  the  angle  of  sliding  friction,  ^ ''6r,  is  the  unit  volume  change  per 
unit  axial  strain,  and  r is  a dimensionless  constant. 

For  triaxial  compression  test,  they  arrived  at  the  expression 


1’C 


3- 


- 1)  (1  - sin^  ) 

s 


] 


where  is  the  axial  unit  strain  contributing  to  volume  change.  The 
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equation  was  tested  by  plotting  + (1  + Sv/fic^)  versus  cr1/o?  + 

(1  + 5v/5f,)  which  should  give  a straight  line  of  slope  arctan  sin,- and  an  inter- 
cept of  f-  All  the  graphs  showed  linear  portions  with  slopes  consistent  with 
mineralogicai  composition,  suggesting  establishment  of  an  equilibrium  inter- 
locking parameter  ~ prior  to  dilation.  After  dilatant  expansion,  a new  line 
is  sometimes  established  at  a lower  f.  but  still  with  the  same  slope,  further 
supporting  that  this  gives  an  independent  measurement  of  sliding  friction.  From 
the  plots  it  becomes  evident  that  if  the  coefficient  of  sliding  friction  is  constant, 
then  the  interlocking  function  ''first  increases  with  increasing  strain,  and 
then  decreases  upon  dilation  and  failure. 
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TEST  METHODS  AND  PROCEDURES 


The  Biaxial  Test  Apparatus 

Brief  description 

Briefly,  the  biaxial  load  test  apparatus  which  formed  the  core  of  the  ex- 
perimental program  consists  of  a horizontal  load  frame  with  an  included 
horizontal  removable  Teflon-lined  test  bed  to  contain  an  assemblage  of  ver- 
tical roils.  Controllable  confining  stresses  are  introduced  along  the  sides  of 
the  bed  by  aluminum  plates  and  pressure  cylinders.  An  axial  strain  is 
applied  at  one  end  of  the  test  bed  by  a screw  and  jack  arrangement  powered  by 
a 1/6  h.p.  motor,  the  axial  force  being  monitored  with  a Dillon  load  cell 
(Fig.  2).  Axial  and  lateral  deformations  are  measured  with  0.001-in.  mechani- 
cal dial  gages  (known  as  Ames  dials),  a 35mm  camera  being  used  to  simul- 
taneously record  the  dial  readings  and  positions  of  the  test  rods  after  every 
0. 005-in.  axial  strain  in  the  beginning,  and  0.  010-  and  0. 020-in.  during  the  later 
stages  of  each  experiment.  The  axial  deformation  rate  is  kept  constant  at 
0.  00405-in.  per  minute  with  a i 2.  5%  variation.  This  gives  a strain  rate  of 
0. 001%  per  second  with,  of  course,  a t 2. 5%  variation. 

Collapsible  box 

The  collapsible  box  (Fig.  3)  consists  of  three  aluminum  5-in.  x 4-in.  x 
1/2-in.  rectangular  plates  on  each  side  of  the  box  and  a bottom  plate  of  Teflon 
16-in.  x 8-in.  x 1/4-in.  supported  by  a plate  of  steel  16-in.  x 8-in.  x 1/16-in. 
The  steel  plate  rests  on  a layer  of  steel  ball  bearings  held  separated  by  an 
aluminum  spacer.  These,  in  turn,  rest  on  an  adjustable  horizontal  rest  plate 
supported  by  the  main  frame.  With  this  arrangement,  the  bottom  of  the  box 
has  Ihree  degrees  of  freedom  of  motion,  namely,  axial,  lateral  and  vertical, 
and  may  be  accurately  levelled.  Each  side  plate  of  the  box  is  connected  to 
four  pressure  units  through  ball-and-socket  joints  at  the  ends  of  the  piston 
rods  to  facilitate  free  rotation  of  the  plate  in  two  directions.  One  end  of  the 
box  consists  of  a 5-in.  x 4-in.  x 1-in.  steel  plate  rigidly  connected  to  the  main 
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Figure  2.  Biaxial  test  apparatus,  fluorescent  lights  and  recording  camera 
at  top. 
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BALL  BEARINGS  TEFLON  AND  STEEL  SUPPORT  PLATES 


Figure  3.  Collapsible  box;  (a)  plan,  (b)  cross  section. 
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frame,  and  the  other  consists  of  a similar  steel  plate  which  is  connected  to 
the  load  cell  by  a ball-and-socket  arrangement. 

Lateral  stress  system 

Constant  lateral  stress  is  applied  through  a system  of  twelve  pressure 
units,  four  for  each  plate,  on  each  side  of  the  collapsible  box.  Each  pres- 
sure unit  consists  of  a pressure  cylinder  and  a piston  fitted  with  a 
"Bellofram  rolling  diaphragm”,  used  because  it  is  capable  of: 

1.  providing  a leak-proof  device  to  convert  gas  or  liquid 
pressure  into  a linear  stroke, 

2.  tolerating  minor  eccentricities  and  cocking  of  the  piston 
rod  and  cylinder  without  affecting  the  operation  of  the  unit, 

3.  providing  a long  stroke  within  a relatively  confined  area, 

4.  responding  to  small  pressure  variations  because  of  very  low 
friction  and  hysteresis,  and 

5.  providing  a working  pressure  area  which  is  constant  (within  + 1%) 
through  its  entire  range. 

All  the  pressure  units  on  both  sides  of  the  collapsible  box  are  connected  to 
the  same  pressure  source.  For  a low  lateral  pressure,  upto  115  psi  gage  or 
30  psi  lateral  pressure  on  the  assemblage,  compressed  air  is  used,  whereas 
for  higher  pressures,  liquid  carbon  dioxide  (CC>2)  is  used.  Liquid  CO 2 
is  capable  of  providing  a constant  pressure  of  up  to  900  psi  at  room  temperature. 
A maximum  gage  pressure  of  200  psi  was  used  in  the  present  experimentation, 
dictated  by  the  capacity  of  Bellofram  rolling  diaphragms.  It  has  been  possible 
to  keep  the  pressure  variations  to  a maximum  of  +.1%  by  the  use  of  precision 
pressure  reuglators  with  operating  ranges  selected  for  the  different  lateral 
pressures  used. 

Calibration  of  Bellofram  pressure  units 

Pressure  units  were  calibrated  in  sets  of  two.  Each  set  was  subjected  to 
gage  pressures  from  0 to  280  psi  increasing  at  an  interval  of  20  psi,  and  from 
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280  psi  to  0 psi  decreasing  at  intervals  of  50  or  20  psi.  A previously 
calibrated  proving  ring  was  used  to  establish  the  relation  between  gage 
pressure  and  force.  The  proving  ring  was  calibrated  by  use  of  a direct 
load,  varying  from  0 to  900  lbs.  — increasing  and  decreasing  at  equal  inter- 
vals of  100  lbs,  on  a calibrated  platform  scale.  Calibration  of  the  pressure 
unity  was  repeated  five  times  for  each  set  to  obtain  average  values.  The 
lateral  stress  c*  and  gage  pressure  correlated  through  the  relation 
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Figure  4.  Lateial  and  axial  deformation  measurement  and  the  axial 
Ic id  measurement  system. 


a3  - 9.2552  (gage  pressure)  + 0.7188 


gage  pressure  = 3.9185  .*t3  - 2.817. 

Lateral  deformation  measurement  system 

Lateral  deformation*.  iire  measured  at  four  points  - two  at  the  ends  and 
two  Hi  between  along  the  length  of  the  collapsible  box  (Fig.  4).  Ames  dials 
capable  of  reading  to  0.  001  in.  with  a total  range  of  1 in.  were  used  to  mea- 
sure separation  of  points  on  two  opposite  sides  of  the  collapsible  box.  The 
dial  stems  were  extended  by  brass  rods  of  9-in.  length  and  1/8-in.  diameter 
threaded  at  both  ends.  Coupling  to  the  box  was  by  two  steel  stops,  one  fixed 
to  the  bottom  of  the  dial  and  the  other  to  the  brass  rod.  The  latter  stop  can  be 
fixed  at  nnv  position  of  the  rod  along  its  length,  thus  giving  a capability  fo 
measure  across  the  collapsible  box  with  variable  initial  width. 

In  order  to  read  all  four  Ames  dials  simultaneously,  along  with  axial  de- 
formations and  loads,  a photograph'..':  arrangement  was  made  at  the  top  of  the 
collapsible  box.  The  arrangement  consists  of  two  concentric  neon  light 
tubes  with  a 35mm  Kodak  f2.  S/SOmm  camera  at  their  center.  The  camera 
is  capable  of  rotation  is  a vertical  plane.  The  assembly  of  light  and  the 
camera  (Fig.  2)  is  mounted  on  a horizontal  leg  and  can  slide  to  any  posi- 
tion along  the  leg.  The  horizontal  leg  can  revolve  around  a column  which, 
in  turn,  is  fixed  to  the  mo.in  frame.  The  horizontal  leg  also  can  revolve  in  a 
horizontal  pl.ne  and  can  move  up  anc:  down  against  the  column.  With  this  ar- 
rangement, the  camera  can  be  easily  focused  and  the  light  inten^'ty  can  be  ad- 
justed fo.*  good  photography.  Photographs  are  mad;  at  regular  time  inter- 
vals, allowing  ail  the  Ames  dials  to  lie  recorded,  w*  ile  the  plan  deformation 
of  the  assemblage  is  photographed,  enabling  location  of  the  plane  or  zone  of 
failure  and  its  inclination.  This  also  facilitates  the  study  of  translational 
and  rotational  movements  of  each  individual  rod,  made  possible  by  marking 
lines  on  the  exposed  rod  ends  prior  to  starting  each  experiment. 
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Axial  deformation  system 

The  axial  deformation  system  consists  of  a Norton-Duff  ten-ton  jaes. 
driven  by  an  electric  motor  of  1 t\  b.p.  through  a complex  speed  reduction 
arrangement.  Application  of  the  axial  load  is  thus  strain-controlled. 

The  electric  motor  is  fitted  with  a precision  motor  speed  control  which 
is  ‘heoretieally  capable  of  regulating  the  speed  in  100  steps.  The  speed 
reduction  system  consists  of  two  Boston  reducers,  five  sprockets  of  30, 
no,  15,  9 and  112  teeth,  and  a worm  gear  jack.  The  Boston  gears  reduce 
the  speed  by  2,000  times,  and  the  sprockets  can  further  effect  a reduction 
of  12  times.  Thus,  the  whole  speed  reduction  system,  including  the  worm 
gear  jack,  is  capable  of  providing  an  axial  deformation  range  of  0 to  2.5- 
thousanc'ths  of  an  inch  per  minute.  An  Ames  dial  is  connected  to  the  worm 
shaft  of  the  jack  such  that  it  directly  measures  the  ial  deformations  in  thou- 
sandths of  an  inch. 

Axi.".l  load  measurement 

For  measurement  of  axial  load,  a Dillon  load  cell  with  its  readout 
system  is  used.  The  load  cell  is  of  a 10, 000-lb.  capacity  and  is  mounted  on 
the  end  of  the  worm  shaft  of  the  lo-ton  jack.  The  other  end  of  the  load  cell 
i»  connected  to  the  collapsible  box  through  a ball  and  socket  arrangement.  The 
Dillon  readout  system  is  in  two  stages,  namely  0 to  5,000  lbs.  and  5,000  to 

10.000  lbs. , which  gives  better  precision  in  reading  the  axial  load. 

Main  frame 

The  main  frnnv  (Fig.  2)  was  designed  as  a rigid  closed  box  section  for 
an  axial  load  of  15,000  lbs.  and  a lateral  load  of  5,000  lbs.  The  maximum 
deflection  in  the  be  section  was  kept  equal  to  one  thousandth  of  an  inch  under 
the  above  load  system.  In  practice,  the  maximum  axial  load  used  is  less  than 

1.000  lbs. , width  gives  a calculated  maximum  structural  axial  deformation 
equal  to  one-quarter  of  a thousandth  of  an  inch,  compared  to  a total  axial 
deformation  of  1 '2  in.  in  the  assemblage.  The  box  section  was  rigidly  welded 
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to  a four -legged  skeleton  supported  on  rollers,  which  gives  further  rigidity 
to  the  box  section. 


Materials  Tested 

The  granular  assemblage  is  modeled  by  cylindrical  rods  of  different 
diameters,  materials,  and  shapes.  Two  different  materials,  namely  steel 
and  Teflon,  were  available  in  geometrical  shapes  of  cylinders  and  square  bars. 
The  following  fiv''  sizes  of  steel  rods  of  circular  cross  section  were  experi- 
mented with: 

1"  diameter  cold  rolled  steel  rods 
3/4*'  diameter  cold  rolled  steel  rods 
1/2"  diameter  cold  rolled  steel  rods 
1/4"  diameter  eteel  welding  rods 
1/8"  diameter  steel  welding  rods 
3/4"  diameter  Teflon  rods 

All  the  sizes,  except  the  1/8-in.  diameter  welding  rods,  were  machine  cut  and 
lightly  polished  on  a lathe  with  a fine  emery  cloth.  The  rods  then  were  covered 
by  a thin  layer  of  light  machine  oil  to  prevent  rusting.  The  lengths  of  the  rods 
'*ve  accurately  cut  so  as  to  give  a reasonably  smooth  surface  while  standing 
their  axes  vertical.  The  1/8-in.  steel  welding  rods  were  saw-cut,  so  the 
len^.h  could  not  be  controlled  as  precisely,  with  the  result  that  a smooth 
planar  surface  was  not  attaino  v.-hen  they  wpre  marie  to  stan.:  with  their  axes 
vertical  in  the  collapsible  boy. 
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Sample  Preparation 

Before  testing,  the  rods  were  washed  with  acetone  to  remove  any  loose 
material  and  lubricant  from  the  surface,  and  were  then  cleaned  with  a clean 
dry  cloth.  A second  acetone  bath  was  given  to  the  rods  before  they  were 
used  for  experimentation. 

Countless  possibilities  existed  for  an  initial  arrangement  of  regular  and 
random  rod  assemblages.  To  simulate-  a densest  packing,  a rhombic  geometri- 
cal array  was  used  in  all  tests,  the  long  axis  of  the  rhombus  coinciding  with  the 
major  principal  stress  Erection  (o^).  Care  was  taken  in  arranging  the  rods 
in  a regular  geometrical  Array  with  its  center  line  coincident  with  the  line  of 
application  of  the  load,  a string  line  being  used  for  this  purpose.  In  a rhombic 
array,  the  rods  in  alternate  rows  number  r and  r-1,  and  it  was  always 
seen  that  the  first  and  the  last  rows  of  the  assemblage  had  the  same  number 
(i.e. , r)  of  l'ods  so  as  to  insure  a symmetrical  assemblage.  Before  applica- 
tion of  axial  strain,  the  assemblage  was  subjected  to  the  required  predeter- 
mined constant  lateral  stress  (a3)  and  the  loaded  end  platen  moved  axially 
with  the  help  of  the  disengaged  jack,  so  that  contact  with  the  load  corresponding 
to  crj  ~ was  shown  by  the  Dillon  load  readout  system.  At  this  point,  the 
axial  deformation  dial  was  set  to  zero.  Lines  making  an  angle  of  60°  with  the 
a„  direction  were  then  drawn  with  a felt-tipped  pen  through  the  centers  of  the 
exposed  rod  ends.  These  lines  helped  in  establishing  movements  of  the  rods 
during  deformation.  After  the  linear  grid  was  marked  with  black  ink,  the 
four  lateral-deformation  modified  Ames  Dials  were  fitted  on  to  rest  on  the 
top  of  the  side  platens.  The  first  lateral  deformation  dial  was  placed  on  the 
center  line  of  the  cylinders  in  the  first  row  from  the  load  side,  and  the  last 
dial  was  placed  on  the  center  line  of  the  last  row;  whereas,  the  second  and 
third  dials  were  exactly  5 inches  apart  on  the  edges  of  the  central  side  platens. 
The  camera  was  adjusted  and  focussed  in  such  a manner  that  the  lateral  defor- 
mation dials  could  be  read  comfortably,  along  with  a reasonably  sharp  picture 
of  the  deformation  pattern.  Special  care  was  taken  in  aperture  and  shutter 
adjustment  for  the  sake  of  neat  photography. 
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Testing  Technique 

Five  values  of  constant  lateral  stresses,  namely  10,  20,  30,  40,  and 
50  psi,  were  used  for  the  experiments.  After  all  preparations  are  completed 
and  the  axial  deformation  dial  is  set  for  an  initio.?  ?ero  deformation  reading, 
a constant  axial  strain  at  the  nominal  rate  of  0. 001%  per  second  is  applied. 
Photographs  are  taken  at  the  interval  of  0. 005  in.  of  axial  deformation  during 
the  initial  stage;  then  this  interval  is  increased  to  0.01  in.  and  finally  to 
0.02  in.  Total  time  for  which  the  axial  deformation  is  applied  to  assemblage 
along  with  the  total  axial  deformation  itself,  is  noted  to  find  the  true  axial 
strain  rate.  After  the  assemblage  has  been  subjected  to  a required  axial 
deformation,  the  experiment  is  stopped  and  the  film  is  processed.  Processed 
film  is  projected  through  an  enlarger  to  read  the  lateral  deformation  dials  and 
observe  the  movement  of  individual  rod  and  failure  planes.  Dial  readings  are 
converted  to  volume  change  by  the  formulations  given  below: 

Effective  width  of  the  assemblage 

= distance  between  the  centers  of  the  first  and  .ast  rods  in  the 
first  row.  With  symbols  as  shown  in  Fig.  5, 

= 2 (r  - 1)  d sin> 

Effective  length  ot  the  assemblage 

= distance  between  the  center  lines  of  the  first  and  last  rows 
= (n  - 1)  d sin  (90  - X) 


Axial  strain  ( r, ) = — 

' (n  - l)d  sin(90  - X) 

Initial  volume  (vQ)  = 2(r  - 1)  n - 1)  d2L  sin(90-X)sinX 


Change  in  volume  (dv)  = ^ 


-A+iLa  + B±Cb  + C + D 


j — e - 2<r  - 1)  dE  sinxll. 


Volumetric  strain  (dv/v^) 


i r A + B_  B + C,  C + D „ ~\ 

2(r  - 1)  (n  - l)d2  sin(90  - X)  sinX  L 2 a + 2 b + 2 e ” 2(r  " 1)dE  8in*  J 
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Axial  stress  (a,)  “ ~ 

2(r-l)d  sinX  L 

where  X = 30  in  the  present  experimentation. 


LATERAL 
DEFORMATION 
DIAL  READINGS 


E AXIAL  DEFORMATION  DIAL 

r = Number  of  rods  in  the  first  cross  row 
n * Number  of  cross  rows 
d = Diameter  of  rods 
L = Length  of  rods 
G=  Load  cell  reodout 
X = Distribution  angle  (Fig.  6). 


Figure  5.  Mathematics  of  collapsible  box. 
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Reproducibility 

The  rate  of  application  of  axial  strain  and  axial  load  was  reproducible  with 
a high  degree  of  accuracy,  and  the  use  of  a two-stage  pressure  regulator  made 
it  possible  to  attain  an  almost  constant  lateral  pressure.  The  reproducibility  of 
the  regular  geometrical  array  with  its  center  line  coincident  with  the  line  of 
application  of  load  was  somewhat  less  precise,  in  spite  of  the  extreme  care 
taken  in  arranging  the  rods.  In  addition,  there  was  unavoidable  variation  in 
the  surface  finish  of  the  rods. 

Another  important  factor  which  may  affect  reproducibility  is  the  initial 
zero  setting  of  the  axial  deformation  dial.  It  was  observed  that  after  applica- 
tion of  constant  stress,  the  axial  deformation  dial  could  not  be  uniquely  set 
for  load,  corresponding  to  a.  = c„  on  the  Dillon  load  cell  readout  system. 

i «> 

For  example,  after  setting  the  axial  dial  at  zero  with  the  readout  load  cor- 
responding to  q.  ~ a_,  if  the  assemblage  was  subjected  to  some  axial  load 

* iJ 

and  left  for  some  time,  then  after  bringing  the  axial  load  to  its  original  value, 
the  axial  deformation  dial  would  not  come  to  the  initial  zero  setting.  This 
effect  was  particularly  pronounced  in  smaller  diameter  rods.  Some  of  the 
possible  reasons,  in  the  case  of  the  smaller  diameter  assemblage,  can  be.* 

(1)  it  is  extremely  difficult,  if  not  impossible,  to  align  the  center  line  of  the 
assemblage  along  the  Cj-axis;  and  (2)  the  smaller  diameter  rods,  being 
relatively  more  slender  than  the  larger  diameter  rods,  can  bend  more  easily 
along  their  length,  thus  rendering  the  perfect  fit  of  the  assemblage  impos- 
sible. 

This  could  be  a reason  for  a floating  respect  to  axial  strain, 

or  a shifting  dv/vQ -versus curve  along  the  direction  of  the  axial  strain  on 
plots.  These  effects  were  apparent  from  repeated  tests.  Flotation  of  the 
(ct,  /"*  ) point  and  the  dv/v  -versus  - * curves  with  different  lateral  stresses 

can  also  be  due  to  different  initial  elastic  deformations,  since  increasing  lateral 
stresses  will  induce  greater  initial  elastic  deformations. 
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THEORETICAL  INVESTIGATIONS 


General  Considerations 

When  an  assemblage  of  rods  is  subjected  to  a stress  field,  then  the 
forces  between  the  cylinders  are  at  first  indeterminate  because  each 
cylinder  is  in  contact  with  six  neighboring  cylinders.  Any  small  change 
in  the  geometry  of  i's  boundary,  e.  g.  a change  in  the  shape  of  the  sur- 
rounding frame  which  is  brought  about  by  moving  the  load  side  platen, 
results  in  a geometrical  change  of  pattern  which  always  includes  a series 
of  gaps  or  slack  contacts  (this  follows  from  Osborne  Reynolds'  theory  of 
dilatancy).  A very  small  change  of  shape  in  the  boundary  is  sufficient 
to  produce  this  effect;  it  needs  to  be  only  of  the  same  order  of  magni- 
tude as  the  elastic  strain  in  the  cylinders.  The  gaps  are  definite  limits 
to  the  lines  of  action  of  the  forces.  Thus,  creation  of  two  gaps  around 
a cylinder  will  reduce  the  neighboring  contacts  from  six  to  four,  which 
will  make  the  forces  between  the  cylinders  determinate.  Since  dila- 
tancy may  be  looked  upon  as  a cause  of  gaps  and  also  of  slack  contacts 
(no  visible  gaps),  its  intimate  connection  with  the  pressure  distribution 
is  apparent. 


Forces  Acting 

Let  us  consider  an  assemblage  of  cylinders  having  a distribution 
angle  equal  to  X (angle  for  some  skew  arrangements  such  that  X = 30° 
gives  us  a hexagonal  array  and  X = 45°  gives  us  a cubic  array,  shape 
factor!  = 1).  When  it  is  subjected  to  a biaxial  stress  field,  a sort  of 
mass  transmission  of  forces  takes  place  in  the  axial  and  lateral  direc- 
tions. If  a cylinder  K (Fig.  7)  in  an  upper  layer  is  supported  by  two 
cylinders  P and  L in  the  lower  layer,  cylinder  K tends  to  push  cylin- 
ders P and  L apart,  tending  to  break  the  contact  between  them.  Now 
consider  a array  (distribution  angle  = x)  of  cylinders  in  which 
cylinder  P is  surrounded  by  six  cylinders,  L,  M,  N,  O,  J,  and  K 
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Figure  fi.  Regular  geometrical  arrays  ;i  is  the  shape  factor  and  > is  the 
distribution  angle;  (a)  hexagonal  array  of  oval  shapes,  (b)  hex 
agonal  array  of  circular  shapes,  > - 150  , (c)  cubic  array  of 
circular  shapes,  1 -If)  , (cl)  hexagonal  array  shown  with 


(Fig.  7).  This  is  the  geometric  arrangement  which  gives  the  minimum 
void  ratio  when  X = 30°.  In  the  axial  direction,  rods  J and  K are 
supported  by  rods  0 and  P,  and  P and  L,  respectively.  Similarly, 
rod  P is  axially  supported  by  rods  N and  M.  In  the  lateral  direc- 
tion, rod  O is  supported  by  rods  J and  N,  and  rod  P by  K and  M, 
assuming  no  contacts  between  O and  P,  and  P and  L.  This  gives 
rise  to  the  mass  transmission  of  forces  from  cylinder  to  cylinder.  If 
f is  the  axial  force  on  one  cylinder  and  f3  is  the  lateral  force,  then 

f,  = 2d  ^ cos(90  - » (1) 

f3  = 2d  c3  cos  X (2) 

where  d is  the  diameter  of  cylinders,  ^ and  a3  are  axial  and  lateral 
stresses  on  a gross  area  basis,  respectively,  and  X is  the  distribution 
angle. 
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Next,  consider  a cylinder  P of  the  array  and  resolve  the  axial  and 
lateral  forces  f^/2  and  fs/2  to  give  a normal  force  F acting  along  the 

(90  - X)  -direction  from  the  CT3  -axis,  and  a tangential  force  T at  a 
right  angle  to  the  (90  - » -direction.  Then 

fa  U 

F =— cos>  4 — sinX 


Substituting  the  values  of  f and  f3  from  Eqs.  (1)  and  (2),  respectively, 
then 


F = da1cos(90  - X)cos>  4 da3cosX  sinX 
= d cosX  sinX  + cr3) 


,sin2x 

= d—  (Oj  + ob 


(3) 


T =2"  sinX  --5- cosX 


= dcijCos^O  - X)sinX  - da3  cosXcos* 

= d^si^X  - da3cos2X 

= d(a1sins>  - o2cossX)  (4) 

The  remainder  of  the  theoretical  development  is  rather  long,  and  is 
presented  in  Appendices  to  this  report. 

The  results  may  be  summarized  as  follows: 

Elastic  Strain 
*dv^e  = 4(o;  4 6tanX) 

vq  d3sin2x  (10) 


where  (dv)e/vQ  is  the  elastic  volumetric  strain  of  the  assemblage  and 
d and  6 are  radial  and  tangential  deformations  at  cylinder  contacts. 
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Analytical  solutions  for  6 and  ot  for  cylindrical  rods,  to  the  best 
knowledge  of  the  authors,  are  not  available.  As  previously  shown, 

F - d^p-  H + n3)  (3) 

T = d(at  sin  ?X  - a<,cossA)  (4> 


By  definition  within  the  elastic  regon,  ^/n3  vai  .^-s  linearly  with  ^ ; this 
also  is  seen  in  the  experiments.  If  c 3 is  kept  constant  throughout  the 
experiment, 


/T 

"1 


E f 


(Hi 


where  E^  is  a constant  and  can  be  called  an  elasticity  modulus  of  the 
assemblage. 

From  the  above,  it  can  be  seen  that  for  a particular  diameter  d 
and  distribution  angle  x of  the  assemblage, 

F ■ some  constant  x * 

T = another  constant  x 

Thus,  once  a relation  is  found  between  a,  6,  F,  and  T from  the  theory 
of  elasticity,  one  could  easily  find  a relation  between  (dv)p/v  and  axial 
strain  ^ . 

Dilatancy  Stress  Ratio 


At  failure 


1 

— - (f  A sinicosX) 

03 


cosX 

sin3X 


(26) 


where  f is  the  coefficient  of  sliding  friction. 
Post  failure 


rCOl>[Mn>  sin()  4 + «*(>  + *] 


(29) 


no 


^>^j,i^ji^UlK_4Vy«i<  ^ 5^- g jUt* ^-WV f ?? S4"  i>_  


Thus,  the  stress  ratio  should  change  according  to  the  above  relation  with 
increasing  deformation  after  initial  failure.  The  relation  between  u> 
and  the  axial  unit  strain  €1  may  be  shown  to  be 


where 


t - arc  tan 


E sin* 


- cot  2) 


(39a) 


E = (n  - l)d  fj  cosX 


(42) 


Example  solutions  of  Eq.  29  for  a hexagonal  array  of  cylinders,  i.  e., 
for  A = 30°,  are  given  in  Table  1.  For  a maximum  stress  ratio,  uj  = 0 
and  Eq.  29  reduces  to  Eq.  26. 

Complete  derivations  and  intermediate  equations  are  given  in 


Appendix  A. 


Tabic  1.  Theoretical  maximum  stress  ratios  for  hexagonal  array  of 
cylinders 


Coefficient  of  sliding  friction, 
f 

Maximum  stress  ratio, 

(cTj/03) 

0. 00 

3.00 

0.01 

3.07 

0.  05 

3.35 

0.10 

3.69 

0. 12 

3.33 

0. 14 

3.97 

0.16 

4.11 

0.18 

4.24 

0. 20 

4.39 

0. 22 

4.52 

0.24 

4.66 

0.26 

4.80 

0.28 

4.94 

0. 30 

5.08 
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Volumetric  Strain  as  a Function  of  A rial  Strain 


On  the  basis  of  observed  failures  along  discrete  failure  planes,  dv/v0 


-1 

\ _il 

(43) 


— — , { r sin(2X  + arctanf- r— — ~ - cot2>l  ) -ll 

n - 1 \sin2X  ' L(n  - 1)  f1sin2X  J ’ ) 

Where  n is  the  number  of  cross-ro  v.3.  The  derivation  of  this  equation  is 

given  in  the  Appendix. 

Effects  of  multiple  failure  planes  and  particle  (rod)  size 

The  coefficient  of  friction  on  the  surfaces  of  the  particles  (rods)  is 
not  constant,  but  is  randomly  distributed  both  in  space  and,  as  slipping 
occurs,  in  time.  Therefore,  the  resolved  total  friction  will  be  dif- 
ferent on  the  different  possible  failure  planes,  failure  following  a plane 
with  the  lowest  resistance.  As  slipping  occurs  on  a failure  plane,  which 
may  be  called  the  active  shear  plane,  the  coefficient  of  friction  on  that 
plane  will  change,  whereas  the  coefficient  of  friction  on  other  possible 
shear  planes  will  remain  constant  due  to  absence  of  movement  along  them. 
Simultaneously  with  increasing  t,  the  value  of  axial  stress  ^ decreases 
(Eq.  29).  Thus,  in  the  process  of  slipping  along  the  active  shear  plane, 
if  the  friction  along  it  momentarily  exceeds  friction  on  another  possible 
failure  plane,  slip  will  become  arrested  and  shift  to  the  plane  with 
the  lowest  friction,  where  the  process  will  vepeat.  However,  as  the 
axial  strain  progresses  and  axial  stress  ^ reduces  considerably  fol- 
lowing Eq.  29,  variation  in  friction  eventually  will  be  insufficient  to 
trigger  slipping  on  a fresh  plane,  and  slipping  will  continue  on  the  same 
plane.  Hence,  with  increasing  axial  strain,  the  process  of  shifting  to 
new  failure  planes  slows  down  and  ultimately  stops. 

Since  the  tendency  for  arrested  slip  to  occur  depends  on  variability 
of  the  cumulative  friction  along  each  available  slip  plane,  it  may  be 
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FT? 


I 


: 


affected  by  element  size,  in  that  the  smaller  the  elements,  the  more  planes 
are  available  for  slipping  to  occur.  The  number  of  planes  p may  be  esti- 
mated from  the  number  of  rods  r in  a cross  row  and  number  of  cross  rows 
0: 

P ” ■iV^  + r - 1 (-13a) 

This  assumes  no  dead  zones  at  the  ends  of  the  box,  and  no  cross  shears,  which 
would  double  p.  Cross-shears  were  occasionally  observed,  but  were  not 
common  due  to  blocking.  Following  are  the  planes  possible  in  the  experiments: 


Rod 


s.ljp 

Planes,  p» 

contacts 

d,m. 

r 

Knd 

Central 

Total 

per  plane 

1.0 

5 

8 

4 

12 

8 

0.  75 

7 

12 

5 

17 

12 

0.50 

10 

18 

17 

35 

18 

0. 25 

20 

38 

34 

72 

37 

0.125 

40 

78 

or, 

144 

77 

•May  be  multiplied  by  2 for  cross-shears 
Probabilities 

The  probability  of  a particular  number  of  slip  zones  developing  during  a 
♦est  may  be  calculated  from  the  number  of  possible  combinations.  For  exam- 
ple, let  us  assume  that  due  to  end  restraint  only  the  four  central  planes  in 
the  array  of  one  inch  rods  are  eligible  for  slip.  There  is  an  equal  likelii  'od 
that  slip  will  initiate  on  any  one  of  the  four  - 1,  2,  3,  or  1,  whereupon  fric 
will  increase  and  slip  will  most  likely  become  arrested.  The  second  slip 
then  may  involve  several  combinations:  1-2,  1 ‘3,  1 • 4 , 2 ~ 3 ,2  - *1 , 3 s 4 . 
Similarly  a third  slip,  if  it  occurs,  could  involve  1 4 2 • .3,  1 • 2 4 -1,  1 * 3 - 4, 

2 • 3 4 -1,  whereas  a fourth  slip  could  only  involve  1 4 2 » 3 * 1.  In  total  there 
are  -1  1 0 •*  4 * 1 15  possibilities.  If  we  assume  that  failure  involved  an  equal 

likelihood  for  any  of  the  15,  the  probability  of  one  plane  at  failure  is  4 f 15  0.27; 
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affected  by  element  si/.e,  in  that  the  smaller  the  elements,  the  more  planes 
are  available  for  slipping  to  occur.  The  number  of  planes  p may  be  esti- 
mated from  the  number  of  rods  r in  a cross  row  and  number  of  cross  rows 
n: 


This  assumes  no  dead  /.ones  at  the  ends  of  the  box,  and  no  cross  shears,  which 
would  double  p.  Cross-shears  were  occasionally  observed,  but  were  not 
common  due  to  blocking.  Following  are  the  planes  possible  in  the  experiments: 


d,m. 

r 

Slip 

End 

Planes,  p* 
Central 

Total 

Rod 

contacts 
per  plane 

1.0 

5 

8 

a 

12 

8 

0.  75 

7 

12 

5 

17 

12 

0.  50 

10 

IS 

17 

35 

18 

0.25 

20 

ns 

34 

72 

.17 

0. 125 

40 

78 

00 

144 

77 

"May  be  multiplied  by  2 ror  cross-shears 
Probabilities 

The  probability  of  a particular  number  of  slip  zones  developing  during  a 
test  may  be  calculated  from  the  number  of  possible  combinations.  For  exam- 
ple, let  us  assume  that  due  to  end  restraint  only  the  four  central  planes  in 
the  array  of  one  inch  rods  are  eligible  for  slip.  There  is  an  equal  likelihood 
that  slip  will  initiate  on  any  one  of  the  four  -1,2,  .1,  or  1,  whereupon  friction 
will  increase  and  slip  will  most  likely  become  arrested.  The  second  slip 
then  may  involve  several  combinations:  1 • 2,  1 • .1,  1 • -1,  2 - 3,  2 4 1,  3 • -1. 
Similarly  a third  slip,  if  it  occurs,  could  involve  1 * 2 • 3,  1 • 2 • 4,  1 • 3 • 4, 

2 • 3 • 4,  whereas  a fourth  slip  could  only  involve  1 1 2 » :i  • 4.  In  total  there 
are  4 * 0 * 4 ♦ 1 15  possibilities.  If  we  assume  that  failure  involved  an  equal 

likelihood  for  any  of  the  15,  the  probability  of  one  plane  at  failure  is  4 4 15  0.27 
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of  two  planes,  0.40;  of  three  plane?,  0.27;  and  of  four  planes,  0.  07.  Based 
on  these  assumptions  the  most  likely  number  of  slip  planes  is  two,  but  there 
is  still  a 60  percent  chance  that  the  number  involved  will  be  other  than  two. 

The  above  relation  involves  elementary  probability  theory,  whereby  the  n 
number  of  combinations  is 


CP  = Ei 

m m!  (p  - m)! 


(43b) 


where  p is  the  number  of  planes  available  and  m is  the  number  participating. 

In  the  above  example,  for  two  active  planes 

r4  4 • 3 1 2 ' 1 
2 2 * 1 (2  * 1) 


The  total  combinations  possible  is  then 

p _ p!  p!  p!  p! 

p!=  1!  (P-  1)!  4 2!  (p  - 2>!  " 3!  (p-3)!  + +p!  (p-p)! 

Where  0!  is  assigned  a value  of  1.0.  In  the  above  example 

4 4!  Ji  . 4|  4! 

•1 ! = 3!  + 2(2)  r 3!  + 4!  (0)! 

= 4-r6  + 4-fl  = 15. 


(43c) 


The  probability  of  occurrence  of  any  particular  number  of  acting  planes  m 

may  be  calculated  by  dividing  CP  by  CP, . This  has  been  done  for  several 

m p! 

values  of  p in  Table  2.  The  highest  probability  in  each  column  is  underlined. 
In  fhe  case  of  one  inch  rods  with  p = 4,  the  most  likely  number  is  2,  whereas 
for  0.  r.  - inch  rods,  p = 12,  and  the  most  likely  number  is  6,  predicting  the 
typical  results  presented  in  Fig.  9, 

The  data  of  Table  2 are  plotted  on  a cumulative  basis  in  Fig.  18. 
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Table  2.  Theoretical  probabilities  of  multiple  slip  planes 


m 

4 

5 

12 

17 

34 

1 

0,267 

0.161 

0.  003 

* 

* 

2 

0.400 

0.322 

0.  016 

0.001 

* 

3 

0.267 

0.322 

0.054 

0. 005 

* 

4 

0.067 

0.161 

0. 121 

0.  018 

* 

5 

0. 032 

0.193 

0.047 

♦ 

6 

0.226 

0.094 

* 

7 

0.193 

0.148 

* 

8 

0.121 

0.185 

* 

9 

0.054 

0.185 

0.  003 

10 

0.  016 

0.148 

0.007 

11 

0. 003 

0.094 

0.015 

12 

* 

0. 047 

0.028 

13 

0.018 

0.048 

14 

0.005 

0.071 

15 

0.001 

0.  095 

16 

* 

0.113 

17 

* 

0. 120 

18 

0. 120 

L 


117 


Figure  8.  Typical  experimental  slip  zones,  (a)  1-in.  steel  rods,  two  planes 

(b)  0.5-in.  steel  rods,  six  planes 


from  continued  axial 
strain. 


TESTING  OF  THEORY 


Experimental  Measurement  of  Sliding  Friction 
The  coefficient  of  sliding  friction,  f,  appearing  in  the  stress  ratio 
equations  (26)  and  (29)  was  independently  evaluated  by  the  method  of  Bahadur 
(1970),  also  described  in  Logani  (1973). 

The  experimental  coefficients  of  friction  from  a section  of  the  1-inch 
diameter  steel  rods  varied  somewhat  with  normal  load.  Fig.  10.  The  mean 
from  14  determinations  was  f = 0.1922  t 0. 0052,  the  t entry  indicating 
90%  confidence  limits  on  the  mean  from  experimental  error.  Similar  tests 
performed  on  a section  of  Teflon  rod  also  indicated  a tendency  for  the 
coefficient  of  friction  to  decrease  with  increasing  normal  ad  , the  mean 
and  90%  confidence  limits  from  eight  tests  being  0. 0144  t 0. 0018.  Scanning 
electron  micrographs  of  the  steel  rod  surfaces  showed  major  irregularities 
indicating  that  the  actual  coeffieicnt  of  friction  of  different  areas  on  different 
rods  probably  varies  considerably  more  than  indicated  by  the  + entries. 

To  test  Eq.  29,  which  expresses  a relation  between  stress  ratio  (a:  /a3) 
and  axial  strain  (Cj),  the  equation  is  x - 0.  The  equation  is  reduced  for 
maximum  stress  ratio  (a.  /<y3)  by  substituting  a>  = 0.  The  reduced  equa- 
tion  becomes  the  same  as  Eq.  26.  By  substituting  various  values  of  coef- 
ficients of  friction  in  Eq.  26,  the  corresponding  predicted  maximum  stress 
ratios  were  found  and  tabulated  against  the  coefficient  of  friction 

in  Table  3.  The  .experimental  maximum  stress  ratios  {<jy/o3)  are  tabulated 
against  corresponding  lateral  stresses  (a3)  in  Table  3 for  all  sizes  of 
rods  tested.  A very  slight  linear  trend  of  increasing  (aa  / tf3)max  with 
increasing  lateral  stress  (a3)  has  been  observed,  and  is  particularly 
conspicuous  in  the  case  of  l"-diameter  and  l/8"-diameter  rods;  but,  in 
the  case  of  other  sizes  of  rods,  this  trend  is  very  slight.  Data  with  the 
l/8"-diameter  rods  are  less  reliable,  due  to  the  presence  of  Teflon  side 
plates  in  early  tests. 
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COEFFICIENT  OF  SLIDING  FRICTION 


Figure  10.  Experimental  coefficient  of 

sliding  friction,  steel  on  steel. 


The  mean  experimental  maximum  stress  ratios  and  90%  confidence  limits 
on  the  means  are  4.06  t 0. 15  for  <eel  and  3.22  + 0.05  for  Teflon. 

Equation  26  solved  for  f and  with  > = 30°  becomes 

f=  0.1443  0.433  (26a) 

Substituting  the  values  for  ^ /os  into  this  equation  gives  f = 0. 15  + 0.  02 
for  steel  and  f = 0. 032  + 0. 008  for  Teflon,  reasonably  close  to  the  measured 
values  of  0.19  and  0.014,  respectively,  in  the  one  case  the  comparative 
value  being  higher  and  in  the  other  case  lower. 

Predicted  Stress  Ratio  in  the  Post- Failure  Region 

To  test  Eq.  29  in  the  post-failure  region,  curves  of  the  theoretical 
stress  ratio  against  axial  strain  (e^  have  been  drawn  for  various  levels 
of  coefficient  of  friction.  For  drawing  the  above  curves,  the  axial  strain 
at  the  maximum  stress  ratio  was  assumed  to  be  zero,  because  we  are  assuming 
that  dilatancy  begins  after  the  maximum  stress  ratio  is  reached.  These 
theoretical  curves  are  superimposed  on  experimental  stress  ratio  curves  in 
Figs.  11  - 14,  selected  to  illustrate  different  phenomena.*  This  super- 
position was  done  in  such  a way  that  the  point  on  the  theoretical  stress  ratio 
curve  corresponding  to  x = 0 or,  in  e*her  words,  the  point  of  maximum  stress 
ratiq  coincided  with  the  maximum  stress  ratio  on  the  curve.  The  theoretical 
formulation  for  the  stress  ratio  does  not  take  into  account  arrested  slip,  which 
means  that  slip  is  assumed  to  be  taking  place  only  on  a definite  shear  plane. 

The  theoretical  and  experimental  stress  ratio  curves  are  in  good  agree- 
ment in  the  case  of  larger  diameter  rods,  i.e. , 1"  and  3/4",  where  arrested 
slip  is  less  probable.  In  Fig.  11  can  be  seen  an  initial  slip  at  f = 0. 16  falling 
to  0. 14,  and  then  suddenly  increasing  and  stabilizing  at  about  0.20. 

In  the  case  of  0.5-inch  diameter  rods,  Fig.  12,  the  development  of  new 
shear  zones  causes  the  "'xperimental  curve  to  cross  the  theoretical  curves, 
which  should  be  translated  as  a group  to  the  right  to  account  for 


*Other  experimental  curves  are  presented  in  Appendix  B. 
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Figures  11-M.  Selected  experimental  stress-strain  curves  for 
’•arious  rods  and  several  lateral  stresses. 
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initiation  of  strain  in  a new  location.  Other  curves  with  0. 5-inch  rods 
were  almost  horizontal  because  of  increasing  shear  zones.  This  ten- 
dency was  less  pronounced  with  the  0. 25n  and  0.126-  inch  rods.  Fig.  13, 
because  of  the  more  rapid  decrease  in  a:  /a3  as  a function  of  axial 
strain.  Evidence  for  increasing  shear  planes  is  the  unrealistically  high 
f if  the  theoretical  curves  are  not  translated. 

The  experiments  with  Teflon  rods  were  complicated  by  a much 
larger  elastic  strain,  Fig.  14,  and  suggest  multiple  shear  planes. 

Figures  11-14  are  representative  of  many  tests;  other  curves  are 
in  Appendix  B. 

Two  modes  of  failure  appear  to  be  defined  by  the  experimental  and 
theoretical  curves.  These  are  illustrated  in  Fig.  15.  In  the  first, 
elastic  deformation  is  followed  by  slip  along  a single  zone,  with  the 
coefficient  of  sliding  friction  f relatively  constant,  and  in  the  second, 
shear  slip  transfers  from  plane  to  plane  throughout  the  post-maximum 
stress  ratio  dilation.  That  is,  slip  takes  place  on  one  shear  plane, 
becomes  arrested,  and  shifts  to  another,  and  so  on.  If  the  coefficient 
of  friction  is  still  relatively  constant,  the  experimental  stress  ratio 
curve,  instead  of  following  one  theoretical  f-level  curve,  will  stay  more 
of  less  parallel  to  the  axial  strain  axis.  This  can  be  explained  by  con- 
sidering that,  with  each  shear  plane  abondonment  and  formation  of  a 
new  plane,  the  theoretical  stress  ratio  curve  starts  anew,  i.e.,  is  shifted 
to  the  right  along  the  axis  of  axial  strain.  Thus,  the  failure  point  will 
move  parallel  to  the  axial  strain  axis.  In  the  soil  mechanics  literature^ 
such  a stage  in  the  stress  ratio  versus  axial  strain  curve  is  called  a resi- 
dual stage  of  assumed  constant  volume  stage,  whereas  in  the  above 
mechanism,  the  coefficient  of  friction  stays  constant  while  the  volume 
may  go  on  increasing. 

If  the  coefficient  of  friction  does  not  remain  constant  the  point  P 
will  not  move  parallel  to  the  axial  strain  axis,  but  will  shift  (up  or) 
down  to  various  f-level  curves  at  each  particular  axial  strain. 
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(c)  (d) 


Figure  15.  Axial  strain  to  stress  ratio  curves;  (a)  failure  mechanism  1, 
(b)  failure  mechanism  2,  (c  & d)  two  combinations  of  failure 
mechanisms  1 and  2 
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Volumetric  Strain 

On  the  basis  of  experimental  observations,  the  following  postulates 
are  suggested: 

1.  The  rate  of  change  of  the  number  of  failure  planes  is  inversely 
proportional  to  axial  strain.  This  implies  that  the  width  of  the 
failure  zone  increases  sharply  at  first,  and  then  slows  down  with 
increasing  axial  strain. 

2.  For  the  same  strain,  the  smaller  the  diameter  of  the  rods,  the 
larger  is  the  number  of  failure  planes. 

If  in  the  first  postulate  dL  is  the  change  of  number  of  failure  planes, 

L is  the  number  of  failure  plane  at  any  instant,  and  d^  is  the  axial  strain. 
Then  the  first  postulate  may  be  expressed  as  an  equation: 


where  - — is  the  slope  of  the  tangent  to  L -versus  curves  and 
is  a constant.  Rearranging  and  integrating  Eq.  44, 

c dL  = dc* 

oL  = lncj  + In  k 
=lnc:k 

L = clnfik  (45) 

where  k is  a constant  of  integration,  or 


<>  - i ccL 
k 


Constants  c and  k can  be  found  from  experimental  boundary  conditions. 
A verification  of  the  form  of  Eq.  46  will  be  shown  by  a linear  relationship 
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between  L and  in  from  which  k was  evaluated  experimentally,  giving 


% 

■? 

ft 

' 

(• 


L = n 417  Cj  (47) 

, -3 

where  is  the  unit  strain  x 10  . 

To  express  postulate  2 in  a mathematical  form,  a general  tabulation 
of  the  number  of  major  failure  planes,  i.e.  those  contributing  significantly 
.o  volumetric  strain,  was  made  against  axial  strains  for  all  the  experiments. 
On  the  examination  of  this  data,  the  following  generalization  was  made: 

On  the  average,  if  one  failure  plane  is  taking  part  in  l"-diameter  assem- 
blage, then  — , — ; and  —r  failure  planes  will  take  part  in 
and  l/8"-diameter  assemblage,  respectively.  Therefore,  if  L is  the  number 
of  major  failure  planes  and  d (in  inches)  is  the  diameter  of  the  particles 
(rods),  the  above  generalization  can  be  approximately  expressed  as: 


L = (d) 


-0.33 


(48) 


Combining  Eqs.  47  and  48,  we  obtain 

In  417  (j  (49) 

L=  '(d)  °-33 

To  incorporate  the  requirement  of  postulates  1 and  2,  the  final  rela- 
tion for  volumetric  strain  is  obtained  by  multiplying  Eqs.  43  and  49: 


dv 

v 


0 


In  417 


{«iSxs1n(2> 


+ arctanU-l><1sin2X 


oot2>]  )-l} 
(50) 


which  should  hold  outside  of  the  region  of  elastic  deformation.  Elastic 
deformations  have  been  observed  to  be  confined  to  the  initial  region  of 
0-0. 2%  axial  strain. 

Equation  50  was  tested  by  superimposing  calculated  and  experimental 
volumetric  strain  (dv/v^)  versus  axial  strain  fa  carves.  This  super- 
position should  be  done  such  that  the  point  of  zero  volumetric  strain 
(dv/vQ)  on  the  theoretical  curve  coincides  with  the  point  on  the  experimental 
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curve  where  elastic  volumetric  strain  becomes  complete.  This  point 
is  uncertain,  and,  for  the  sake  of  simplicity,  it  is  assumed  to  the  point 
of  minimum  volumetric  strain.  The  difficulty  is  all  the  more  proble- 
matic when  one  considers  that  elastic  and  dilational  deformations  probably 
overlap. 

In  spite  of  the  simplifying  assumptions,  it  can  be  seen  from  the  super- 
imposed curves  that  there  is  fairly  good  agreement  between  the  theoretical 
and  experimental  date  (Fig.  16  and  17).  A slight  deviation  from  the  theo- 
retical curves  has  been  observed  in  the  intermediate  stage  of  axial  strain 
in  that  the  experimental  volumetric  strain  exceeds  the  theoretical  value 
corresponding  to  a particular  axial  strain.  Deviation  from  the  theoretical 
volumetric  strain  (dv/v  ) versus  axial  strain  curve  also  becomes  signi- 
ficant at  very  large  axial  strains  because  of  more  chaotic  and  less  pre- 
dictable movements  of  the  particles  at  large  induced  deformations. 


PROBABILITY  IN  RELATION  TO 
ROCK  CREEP  AND  PROGRESSIVE  FAILURE 


Rock  Creep 

Rock  creep  is  a widely  known  phenomenon,  descriptive  of  progressive 
movements  which  may  lead  to  eventual  catastrophic  failure  of  rockslides, 
cavern  roofs,  etc.  In  this  connection  Tcrzaghi  (1900)  suggests  that  a ran- 
domly jointed  rock  behaves  like  a stiff  clay  or  an  impure  sand  with  considera- 
ble cohesion.  Following  the  behavioral  analogy  to  clay  soils.  Barton  (1972) 
relates  creep  to  a loss  of  cohesion  exemplified  by  a residual  strength  contrasted 
to  a peak  strength.  Brawner  et.  al.  (1972)  and  other  authors  have  recognized 
the  variable  nature  of  friction  in  jointed  rocks,  but  to  our  knowledge  have  not 
attempted  to  relate  this  to  a creep  behavior.  The  zonal  nature  of  rock  creep  is 
suggested  by  Zischinsky  (19(56),  who  differentiates  between  rock  sliding  (Gleitung) 
and  rock  creep  (Sackung)  on  the  basis  that  the  latter  shows  a velocity  dependent 
on  distance  above  stable  rock,  therefore  constituting  a zonal  failure. 

A prevailing  assumption  with  regard  to  rock  creep  appears  to  be  that  clay 
in  the  joints  must  be  a causal  factor.  For  example,  Cording  et.  al.  (1971),  in 
discussingthe  stability  of  caverns  at  the  Nevada  Test  Site  and  other  locations, 
suggests  that  strain  is  a function  of  t exp  m where  t is  lime  and  m is  a 
constant  that  may  relate  to  plasticity  index  of  the  fault  gouge.  Broadbcnt  and 
Ko  (1971)  also  report  gross  semilogarithmic  relationships  between  slope  move- 
ments and  time,  but  analyze  the  data  with  a rheological  model  and  replot  on  a 
linear  scale  to  show  that  movement  actually  is  in  a series  of  jumps  which  may 
or  may  not  show  a progressive  increase  in  magnitude. 

Probability  and  Creep  Rate 

The  calculated  probabilities  for  number  of  shear  planes.  Table  2,  are 
plotted  on  a cumulative  basis  in  Fig.  18.  The  resulting  curves  bear  a striking 
resemblance  to  classical  creep  behavior,  by  considering  that  the  total  shear 
deformation  is  proportional  to  the  number  of  shear  planes,  and  cumulative 
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probability  is  representative  of  elapsed  time.  Both  of  these  assumptions  appear 
reasonable,  in  that  in  a given  assemblage  exhibiting  the  arrested  slip  phenomenon 
the  cumulative  slip  will  equal  the  average  individual  slip  distance  times  the  num- 
ber of  slip  planes,  and  secondly  slip  results  from  minor  fluctuations  in  stress 
or  friction  which  are  random  in  time,  the  cumulative  probability  of  slip  will  bear 
a linear  relation  to  time. 

Let  us  therefore  reexamine  the  factors  and  assumptions  made  to  obtain  data 
in  Table  2 and  Fig.  18: 

1.  Slip  will  initiate  along  the  weakest  plane  or  the  surface  with  the 
lowest  accumulated  friction. 

2.  By  bringing  new  surfaces  into  contact,  slip  will  tend  to  induce  an 
increase  in  friction,  which  usually  will  be  sufficient  to  cause  slip 
to  become  arrested  and  shift  to  the  next  stronger  plane. 

3.  As  a result  of  arrested  and  relocated  slip,  a constant  stress  ratio 
vs.  strain  curve  may  be  developed. 

4.  Each  possible  combination  of  the  potential  failure  planes  either 
singly  or  in  multiple?-  has  an  equal  likelihood  of  occurrence. 

There  can  be  little  argument  that  the  first  assumption  is  valid  for  fractured 
rock  masses,  so  long  as  frictional  resistance  is  taken  to  include  both  sliding  fric- 
tion and  dilatant  effects. 

The  second  assumption  that  sliding  increases  friction,  appears  valid  so  long 
as  the  reduction  is  not  overcome  by  strain  weakening  through  thixotropy  or 
slickensiding.  If  the  second  assumption  is  invalid,  strain  will  continue  along 
the  same  weakest  plane,  with  friction  decreasing  or  remaining  constant.  If 
stress  is  not  somehow  relieved,  the  result  can  only  be  sudden  failure. 

The  third  factor,  development  of  an  interval  of  constant  stress  vs.  strain, 
has  been  demonstrated  for  the  steel  and  Teflon  rods  (Fig.  11,  12,  14)  and  is 
fairly  common  in  triaxial  testing  of  rock. 

The  fourth  assumption  affects  the  shape  of  the  probability  curve.  While  the 
test  data  are  insufficient  in  quantity  to  validate  this  assumption  they  do  comply 
with  >ts  prediction. 
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Perhaps  most  significant  is  that  an  ideal  creep  behavior  is  predicted  for  a 
fractured  rock  assemblage  which  is  devoid  of  cohesion.  The  requisites  include 
a stress  level  close  to  failure,  and  several  planes  available  and  engaged  in 
arrested  slip.  Should  such  an  action  be  involved  in  field  occurrences  of  rock 
creep  along  fractures,  it  may  also  explain  the  subaudible,  clicking  "rock  noises" 
which  characterize  creep  and  warn  of  an  impending  landslide. 

Finally,  if  rock  creep  field  data  can  be  fitted  to  probability  based  curves 
as  in  Fig.  18,  there  should  be  a considerable  advantage  for  prediction  of  ulti- 
mate failure,  compared  to  empirical  semilogarithmic  data  plots  which  predict 
a smoothly  increasing  rate  of  movement  with  no  sudden  changes.  Such  a plot 
of  the  right  half  of  data  in  Fig.  18  shows  two  straight  lines.  Fig.  19,  the  second 
line  being  indicative  of  catastrophic  failure.  The  lines  intersect  at  a cumula- 
tive probability  of  about  0. 99,  where  about  60-70  percent  of  the  available  planes 
are  participating  in  the  failure. 

Thus  it  would  appear  that  the  time  of  failure  may  be  rather  accurately  pre- 
dictable from 'a  changing  thickness  of  the  failure  zone.  This  may  be  monitured 
by  inserting  vertical  pipes  through  the  failure  zone  and  periodically  measuring 
their  curvature  with  depth.  The  zone  thickness  plotted  against  time  may  give 
graphs  as  in  Fig.  18  or  19.  From  Fig.  18  may  be  seen  that  as  a rough  guide,  a 
doubling  in  thickness  of  the  shear  zone  should  be  sufficient  to  bring  about  catas- 
trophic failure. 
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CONCLUSIONS 


P 


Relative  to  Test  Data 

1.  In  a regular  array  representing  a homogeneous,  cohesionless  assem- 
blage, the  stress  ratio  (aj/ca)  is  a predictable  function  of  geometry 
of  the  array  and  the  coefficient  of  sliding  friction  at  contacts. 

2.  Previous  dilatant  work  theories  derived  for  soils  do  not  describe  the 
observed  results  because  they  neglect  elastic  deformations  prior  to 
dilatancy,  and  they  assume  sliding  friction  to  be  a constant. 

3.  After  an  elastic  volumetric  decrease  stage,  dilatancy  initiates  by 
sliding  along  the  least  resistant  chain  of  frictional  contacts.  Friction 
in  this  active  zone  is  thus  less  than  the  average.  Therefore,  as  sliding 
proceeds  and  brings  new  surfaces  into  contact,  friction  in  the  active 
zone  almost  always  increases , tending  to  cause  slip  to  become  arrested 
and  relocated  to  a new  zone. 

4.  A densest  packed  starting  geometry  such  as  the  one  chosen  (a  rhombic 
arrangement  of  cylindrical  particles)  gives  a decreasing  stress  ratio 
(ctj/os)  as  strain  proceeds;  however,  a recurrence  of  arrested  slip 
will  cause  0^/03  to  increase,  decrease  less  rapidly,  or  hold  constant, 
(the  latter  simulating  creep  behavior),  as  the  shear  zone  increases. 
After  the  maximum  shear  zone  becomes  established,  the  stress  ratio 
decreases  in  the  predicted  manner. 

5.  For  volume  change  of  a dilating  array  to  be  predicted  from  geometry 
of  the  array,  allowance  must  be  made  for  the  gradually  enlarging 
shear  zone  caused  by  arrested  slip,  and  the  volume  change  usually  is 
larger  than  would  occur  from  a single  slip.  This  violates  a basic 
premise  of  the  minimum-work  theories. 

6.  The  number  of  slip  planes  appears  to  be  predictable  from  laws  of 
probability,  the  most  likely  number  being  one-half  of  those  available. 
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Relative  to  Rock  Creep 

7.  The  probability  model  for  number  of  shear  planes  reasonably 
simulates  typical  creep  stress-strain  behavior,  the  rate  slowing 
down  at  first  (primary  creep),  then  remaining  constant  (secon- 
dary creep),  then  speeding  up  and  going  into  rapid  failure  (ter- 
tiary creep).  However,  the  mechanism  differs  from  viscous  creep 
in  that  the  particulate  array  may  be  cohesionless,  and  the  active 
zone  is  ever  increasing  in  thickness. 

8.  The  probability  model  for  creep  behavior  requires  that  the  thick- 
ness of  the  shear  zone  be  increasing,  and  that  catastrophic  failure 
becomes  imminent  when  the  thickness  approximately  doubles  its 
initial  value. 
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Appendix  A 

Derivations  of  Equations 
Elastic  Strain 

Consider  cylinder  pairs  J and  P,  and  K and  P,  in  an  array,  Fig.  7. 

A compressive  force  acting  through  the  centers  of  a pair  will  bring  them  closer 
by  a distance  here  designated  a,  known  as  the  distance  of  relative  approach. 

A tangential  force  similarly  will  cause  a tangential  displacement  6.  In  a 
symmetrical  array  the  lateral  components  of  both  ol  and  6 are  equal  and  op- 
posite in  sign  and  add  to  zero  in  the  ?3  direction,  giving  e3C  = 0.  The  sum 
of  axial  components  of  a and  6 is 

_ 2 (geos)-*-  6sin>) 
dcosX 

9 

= T(or+$tan»  (5) 

Since  epQ  - e3C.  - 0,  and  the  elastic  volume  change  per  unit  length  of  one  cy- 
linder is 


(fo)e  - Cj e * f?e  ^ f?e 
2 

= (a  -1-  *an»  (6) 

For  a cylinder  pair, 

<dv)e  - 2(%)e  (7) 

Multiplying  by  the  number  of  cylinders  in  the  assemblage, 

(dv)e  * (n  - 1)  (r  - 1)  - (a  + 6tan»  (8) 

The  orignial  volume  of  the  assemblage  is 


v = (n  - 1)  (r  - 1)  d?  sin2> 

Dividing  (8)  by  (9) 

^ ^e  _ 4 (g  + 6tan>) 
v ' d3sin2> 


(9) 


(10) 
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a,  = a 3 cots>  (16) 

at  which  value  of  > the  force  S is  being  transmitted  along  the  centerline  of 
the  contacting  cylinders.  From  Eq.  (14),  the  four  axial  forces  symetrically 
disposed  in  four  directions  around  P became 


= d / (cr3  co^X^sin2!  + a32  cos2X 

(17) 

= d cotXa3 

(18) 
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=^‘*  U=-" 


r^TTm 


As  the  stress  ratio  further  increases,  this  condition  is  abandoned,  and  the 
normal  component  of  S becomes 


S - S cos  (y  + X - 90) 
n 

The  tangential  component  is  no  longer  zero,  but  is 
St  = S sin  (y - 90) 


(19) 


(20) 


Experimental  data  indicate  that  this  occurs  before  0.1  to  0.2%  axial  strain. 
Eventually  it  is  observed  that  two  contacts,  as  JP  and  PM,  slide,  and  two, 
as  KP  and  PN,  roll,  the  remaining  two  contacts  OP  and  PL  becoming  separated. 
The  limiting  tangential  force  at  JP  and  PM  when  the  contact  stress  is  a maxi  - 
mum  is 


S,  - f(d  cotX)  a, 
i max  3 


(21) 


where  f is  the  coefficient  of  sliding  friction.  The  resulting  couple  is: 


f(d  cotX)  a,d 

which  is  opposed  by  a couple  at  the  other  rolling  contacts: 

S sin(-u+  X - 90)d 

Equating  and  sutstituting  for  S from  Equation  (14)  gives 

f cotX  c,  d2  = d?  sin  (y  + X - 90)  / sin^X  + 03p-  cosyl 

Simplifying  and  utilizing  Equation  (15)  for  X gives 

tanX 

sin  y = 


(22) 


(23) 


and 


cos  y 


/ tan;jA  + 


/■a  * tan*X  + c?- 


which  may  be  substituted  in  the  above  to  give 

— it  | \ rt  rtfl  \ \ _ 

sin3X 


1 ,c  vv  CoSl 

— = (f  +sin  X cosx) 

^3 


(24) 


(25) 


(26) 


(See  alsoLogani,  1973). 
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Post  Failure 

As  the  array  deforms  from  the  initial  angle  X by  an  amount  <x,  the  indi- 
vidual force  s makes  an  angle  (y  + X + -j)  - 90)  with  a line  connecting  centers 
of  pairs.  The  force  given  by  Equation  (18)  will  still  be  acting  and  may  be  re- 
solved as  before.  The  two  opposing  couples  become 

f(d  cot  X)  a3d  (27) 


S sin  (v+  X + a*  - 90)d 
A similar  development  then  gives 


— = COtX  T — ~~r ; + COt(\  + 0))1 

a3  L smx  sin(X  + o ) J 


Volumetric  Strain 

Assuming  the  wedge  ABC  moves  as  one  unit  parallel -to  AB,  the  movement 
of  gf  parallel  to  AB  relates  to  the  axial  movement  by 


’ sin(90-X)  * 

The  lines  connecting  centers  of  a unit  cell  in  the  array,  as  JPNO  in  Figure  7 
form  a parallelogram  with  area  ds  sin2>.  The  area  of  void  involved  is 

d3(sin2X-|).  (31) 

As  X increases  by  the  amount  uj,  the  new  void  area  becomes 

d2[sin(2X  + w)  - -]  (32) 

Subtracting  and  multiplying  the  result  by  the  number  of  distorting  unit  cells  in 
two  shear  planes,  4(r  - 1),  gives 

dv  = 4(r  - l)d3  £sin(2x  + yj)  - sin2xjx  rod  length  (33) 

The  original  volume  is 


vq  = £(n  - l)(r  - l)ds  sin2xjx  rod  length 
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Dividing, 


dv 

v 


4jjsin(2x  + «))  - sin2X  J 
(n  - 1)  sin2X 


(35) 


Since  the  length  of  the  assemblage  is 


Lq  = (n  - 1)  d cosX, 


(36) 


substituting  in  Equation  (34)  gives 


dv  _ 2d  sin(2X  + u?)  - sin  2>. 
v L sin> 


(37) 


The  angle  may  be  solved  as  follows:  Let  the  distance  moved  by  a contact 
be  gg' , the  chord  of  an  arc  which  for  the  usual  two-shear-plane  zone  has  a 
radium  2d.  The  chord  length  is 


gg'  = 2(2d  sin  -) 


(38) 


Since  /gg'  f = (90  - ^)  and  /gfg’  = 2>  + u>,  by  the  law  of  sines 

gg.  Si  sin  gX  + <>') 

sin  (80  -|) 


Combining  with  Equations  (30)  and  (38), 

2(2d  sin|)  = - 31-  ^ +:^ 

sin  (90  - ~)  sin  (90  - » 
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which  simplifies  to  (Logani,  1973) 


tan»>  =[e75T>  - 'o'2']''. 


The  value  for  u>  may  be  substituted  in  Equation  (37)  to  give 
dv  4 ( 1 r d 


- iri  { drr» sin(2>  + aro  la°  [ii5i  - cot  2xT')  -1) 


— -lsin(2l  + arc  tan  f rr~~:  - cot  2)T  )-sin2)\ 

v L sin)  t.  ' L E sin)  J ' J. 


o o 


By  definition  axial  strain  ^ = E 4-  length  of  the  sample  or 


(n  - 1)  d cosX 


which  solved  for  E and  substituted  in  Equation  (40)  gives 


— = f —~r~  sin(2)  + arc  tanf  — r -rr  - cot2)T  )-lj 

v n - 1 i sm2)  L (n  - 1)  sin2X  J J 
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